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ABSTRACT 


' An investigation is conducted into the use of optical data from 
onboard television cameras for the navigation of interplanetary space- 
craft during the planet approach phase* Three optical data types are 
investigated — the planet limb with auxiliary celestial references, the 
satellite-star -and the planet-star two-camera methods. Analysis and 
modelling issues related to the nature and information content of the 
optical methods are examined. Dynamic and measurement system modelling, 
data sequence design, measurement extraction, model estimation and orbit 
determination, as relating to optical navigation, have been discussed. 

The various error sources are analysed. The methodology developed has 
been applied to the Mariner 9 and the Viking Mars Missions. Navigation 
accuracies are evaluated at the control and knowledge points, with parti- 
cular emphasis devoted to the combined use of radio and optical data. 

A parametric probability analysis technique is developed to evaluate 
navigation performance as a function of system reliabilities. 

It has been determined that Optical Navigation can be a very 
effective means of navigating an interplanetary spacecraft during its 
approach phase to the planet, particularly with the combined use of radio 
and optical data. Of the three observation methods examined the 
satellite-star method is found most suited for the knowledge point and 
the planet-star two-camera method for the control point. It has been 
shown that optical and radio data provide complementary navigation 
information and their major error sources are different; their combina- 
tion yields the best results. However, delaying the maneuver timing as 
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much as feasible maximizes the benefit from the strength of optical data. 

A method developed to evaluate consistency between the optical and radio 
solutions is shown to be very effective in the detection of data anomalies. 
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CHAPTER 1 


INTRODUCTION 

It is in the nature of Man to explore his universe, to venture 
into outer space because, like the mountain, ”it is there/* 

Functionally any exploration system must include the following 
three components: 

i> the vehicle or means of transportation along with the 

associated equipment and the investigative instrumentation 
required for the exploration, 

ii) the navigation system — methods for directing the course 
of the vehicle to the desired destination and for pointing 
the instruments towards the desired sources of information,, 
iii) communication method of transmission of commands and 
retrieval of information gathered by the instruments .on 
the craft,- particularly during the navigation of largely * 
non-autqnomous vehicles . 

Each of these three components is indispensable for the success of any 
exploratory project. However, after the vehicle has been built and 
the communication system and hardware set up, the principal job 
remaining is to navigate the craft effectively. At this point, 
assuming that all the hardware performs according to specifications, 
the navigation accuracy constitutes the chief source of uncertainty 
regarding the success of the mission* 
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1.1 Navigation for Space Exploration 

Any navigation system in essence involves the relationship of two 
frames of reference, where one of these must be tied to the vehicle 
being navigated and the other to the destination. Navigation systems 
differ in the method employed to establish this relationship. The 
errors associated with the navigation system are therefore different 
not merely owing to the different characteristic instrumentation errors 
but also through the accuracy of establishing the relationship between 
the two frames of reference. 

The earliest sophistication in the art of navigating a vehicle 
came about on the high seas in the use of the sextant (an optical 
device!) where the measurements essentially determined the orientation 
of the ship relative to the fixed stars in inertial space. This 
information was combined with a knowledge of the earth’s orientation 
relative to the stars to determine the location of the ship on the 
earth, i.e. ,_r_elative to an earth fixed coordinate frame of reference. 
This was possible since marine navigation is a two-dimensional 
problem, with the vehicle always (hopefully!) on the surface of the 
sea, and therefore orientation information is equivalent to location 
information. 

Navigation methods developed along with the vehicles being 
navigated, through gyroscopic devices and acceleration sensors to 
radio navigation systems for aircraft. For space exploration the 
sole mode of navigation has historically been radio navigation 
wherein electromagnetic waves with stable frequency are transmitted 
to the spacecraft from an earth-based tracking station. The 
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spacecraft retransmits the signal back to the tracking station where 
it is received with a change from the original frequency. This 
doppler shift, Af, in the frequency of transmitted electromagnetic 
waves is used to infer p the instantaneous range-rate of the spacecraft 
relative to the earth-based transmitting station, through the 
relationship 



where, 

c = velocity of light 

- transmitted frequency. 

In addition the tracking .station can measure the time the signal 
takes to travel to the spacecraft and back again, which can be used 
as a range measurement from the station to the spacecraft. 

If now the station is related to a geocentric frame of reference 
and this in turn to the target, the relationship between the spacecraft 
frame of reference and the target frame of reference is available. The 
accuracy of this depends upon the cumulative effect of the error 
sources from each of the steps in the process. 

Figure 1-1 illustrates the geometry involved in the acquisition, 
of radio data. The various coordinate frames of reference Involved 
in the process of obtaining target relative spacecraft trajectory are 
indicated. The directions of the axes and the origins of the coordinate 
systems are both relevant; for this simplified interpretation, the 
spacecraft state vector would be referenced to each of these coordinate 
frames successively: 
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S: The spacecraft and its coordinate system 

T: Topocentrlc (station-centered) earth-fixed frame of reference 

E: Geocentric (earth-centered) earth-fixed frame of reference 

G: Geocentric space-fixed frame of reference 

3: Barycentric (Earth-Moon) space-fixed frame of reference 

Heliocentric (sun- centered) space-fixed frame of reference 
P: Target centered frame of reference. 

The purpose of navigation as outlined above is to (1) determine 
the course of the spacecraft with respect to the target along with 
its precision, and (li) guide it to the desired configuration at the 
desired time within the desired accuracy. We restrict ourselves here 
largely to the former task. The problem can be regarded as the 
determination of the location and velocity of the spacecraft at an 
epoch along with a description of the accuracy of this determination. 
Given a model to propagate these through the use of "Newton’s second 
law and the framesof reference involved^ we can determine the target 
relative spacecraft state at any time. 
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1.2 Nature and Information Content of Radio Data 


Space exploration requires very accurate determination and 
control of the location and destination of the spacecraft. Radio 
data has proved to be a powerful method for navigation particularly 
for near-earth spacecraft such as earth-orbiting satellites and 
lunar missions. The problems for interplanetary exploration, 
however, become much more demanding; this is particularly so for 
those space missions where the intent is to place the spacecraft 
into orbit around the target planet. To examine the reasons for 
the additional problems introduced in these situations we first 
discuss a simple model for the information content and strengths of 
radio navigation in this section (originally from References 1, 

2, and 3). The error sources associated with it are discussed in 
Section 1.3, 

Figure 1-2 illustrates the geometry for doppler tracking of 
a spacecraft from a station on the surface of the earth. The 
X y 2 frame of reference is a geocentric space fixed coordinate 
system with the z axis aligned with the spin axis; we disregard 
any errors for the discussion in this section. 

The topocentric position and velocity vectors to the spacecraft, 
and are given by 

£ = £ - 1st (1-2-1) 

and 



( 1 - 2 - 2 ) 
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where 



Letting 

r = ^ = geocentiric range-rate of the S/C, 

a = Earth’a rotation rate, 

A = longitude of tracking station from prime meridian G, 
r^ = distance of tracking -station from -the Earth's spin axis, 

_r^ = unit T/ector in the direction of increasing X, 

•a = geocentric right ascension of the S/C, 

5 = geocentric declination of the S/C, 

i 

and t ■= time elapsed since coincidence of prime meridian with x axiso 
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we may write 


r = ^ ^ (1-2-7) 

and 

ir . • r = (or r, ♦ r • (1-2-8) 

= -air cos6 sin(oit + A - a) . (1-2-9) 

s 

Substituting into equation (1-2-6) , we obtain 

p = r + oir cos6 sin(ait + X - a) , (1-2-.10) 

s 

which is the same result as originally derived by Hamilton and. 
Melbourne in Reference 1. This indicates that the signature of one 
pass of Doppler tracking data (see Figure 1-3) yields the geocentric 
range rate, right-ascension and declination of the spacecraft which 
are implicitly assumed to^ be constant over the pass. This is due to 
the so-called ’Velocity parallax" induced by the motion of the 
rotating tracking, station; attention to this was originally called 
in Reference 2. 


From several days of tracking data it is then reasonable that 
we should be able to determine the time rate of change of these 
quantities — viz the geocentric acceleration and the right ascension 
and declination rates. Reference 3 showed this and pointed 'out that 
the first of these quantities has information in it tO' determine 
the geocentric range of the spacecraft. This was done by extending 
the model to a six parameter model (References 3, 4). 


p - a + b sinwt + c coswt 

+ dt + et sincot + ft coswt , 


( 1 - 2 - 11 ) 


where 
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a = r 

b = o)r cos6 
s 

c ” (or^ cos6(AA - Aa) 

d = a + rCa^cos^S + 6^) (1-2-12) 

s 

m 

e = -wr sin56 
s 

f === -wr {(AA - Aa) sin66 + a cos6} . 

3 

* 

In the expression for the second term is the familiar V^/r 

centripetal acceleration and a is the remaining geocentric acceleration 

S 

due to gravitation, solar pressure etc,, the centripetal acceleration 
term depends on the range, x. This can be supplemented by direct 
measurement of the range p as discussed in Section 1.1. 

Thus all six components of a spacecraft’s geocentric trajectory 
can be determined from a few days of Doppler tracking data. The 
precision, or lack thereof, with which these quantities are determined 
is of course the critical question, and that forms the subject of 
the next section. 
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1.3 Radio Navigation Error Sources and Limitations 

In this section we list the major error sources for radio 
navigation, followed by a simple analysis of the effects of these 
on navigation accuracy. An excellent exposition on this topic 
is given in Reference 5. A slightly different viex^rpoint is taken 
here in presenting the error sources. We categorize the errors 
into three types: 

(i) Measurement imprecision errors 
(ii) Spacecraft acceleration model errors 
(lii) Errors in relating frames of reference. 

The first type of error stems from spacecraft and tracking 
station equipment hardware limitations. The most significant 
contributor of error in this category is the instability of the 
reference frequency in the station master oscillator. Other 
components include cycle slipping, ranging system biases and drifts 
and clock syn chro nization errors. In practice all these errors 
combined turn out to be much less significant than the other 
categories of errors, for the state of the art as it exists ,at 
present. 

Spacecraft Acceleration Model Erro,rs 

This type of error includes any unmodeled or mlsmodeled 
forces that affect the motion of the spacecraft on its trajectory. 

Some of these forces are non’-gravitational, resulting from un- 
certainties in spacecraft associated effects, e.g., solar pressure 
and gas leakage forces. Others are gravitational and arise from 
attractions by massive bodies; if the masses are not known accurately. 
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the computed trajectory will be in error. 

1) Nongravitational: There are two principal nongravitational 

acceleration error sources. The first of these is caused by ■, 
spacecraft generated effects such as gas leaks from the attitude 
control system and spacecraft propulsion system. .These leaks 
arise from imperfections in the valve seats for -these. systems. 

The second kind of nongravitational force is caused by the 
environment such as the radiation pressure due to the solar energy 
falling upon the spacecraft (S/C) structurso 

Solar radiation pressure in the sun-S/C direction is modelled by 

AR = , ‘ (1-3-1) 

MR2 

where 

R IS the probe-Sun distance 

0 

K is a solar radiation constant (= lo031 x 10 ) 

A is the spacecraft effective area normal to the sun-spacecraft 
dir ectio n 

M is the spacecraft mass 

G is the reflectivity coefficient 

R 

Acceleration errors from these sources are very small (“^10 
km/s^) , and are significant not for the amount by which they perturb 
the actual motion of the spacecraft, but rather for the way they 
affect the tracking data — i*e., it is the uncertainty in the 
solar pressure accelerations C^ot the accelerations themselves) 
that impacts navigation accuracy* 

2) Gravitational: In addition to the two types of non- 
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gravitational forces, a third spacecraft acceleration error source 
is gravitational forces stemming from mismodelling (i) the masses 
(as represented through the .gravitational constants) and (ii) the 
harmonic fields of the .perturbing bodies* For the latter, only 
the second harmonic, J 2 (of the Legendre harmonic expansion — see 
Reference 6) is considered to be of any significance for the near- 
planet phases of interplanetary orbit determination. 

Errors in Relating Frames of Reference 

We adopt the same viewpoint here as in Section 1.1, where the 
frames of reference listed were, in order (see Figure 1-1), 
S-X-E-G-B-H-P 
(1) (2) (3) 

to establish the relative vector between the target P and the 
spacecraft S. Errors in this category fall into three subdivisions 
in .accordance with the grouping of transformations indicated ahove^ 

1) Locating the spaceraft in topocentric coordinates: Errors 

in this involve the effects of the intervening medium. The radio 
signal in its path from the station to the spacecraft and back 
again is distorted by the intervening medium through which it 
travels. This occurs due to 

(a) refraction effects in the earth* s txopos,phere which cause re- 
tardation and bending of the electromagnetic beam {Ret. 7); the 
observation appears to have travelled through a longer distance. 

(b) effects due to the radio signal interaction with charged par- 
ticles in the earth^ s ionosphere and the plasma in interplanetary 
space. The variation in particle density causes a time rate of - 
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change in the ray path length (Ref. 8) and thereby introduces 
errors in the radio doppler measurements. 

Representing the modelled spacecraft in topocentric coordinates 
— i.e., relating frames of reference S and T — therefore 
requires calibrating the effects in the transmission media. 

Residual errors after calibration corrupt the information content 
of the radio observations. 

2) Locating the tracking station in geocentric space-fixed 
coordinates: This involves (a) the transformation from a topo- 

centric frame of reference T to a geocentric earth-fixed frame of 
reference E and (b) the rotation from earth-fixed coordinates to 
the space-fixed geocentric frame G. The first consists of estab- 
lishing the actual location of the tracking-station relative to 
the ”earth-crust" by means of latitude, longitude, and radius 
vector, or equivalently in cylindrical coordinates — i.e., distance • 
from spin axis r^, longitude A relative to the prime meridian and 
height above the equator 

Errors in the rotation from E to G consist of three effects 
as follows : 

(i) Polar motion — the motion of the earth’s crust relative to^ 
the spin axis which causes a change of position of the point 
on the earth’s crust through which the spin axis passes 
(Ref . 9) . 

(ii) Timing irregularities caused by non- uniformities in the 
rotation speed of the earth. These errors can cause a 
degradation in the tracking data quality (Ref. 10). 


PXGEJ is 
OK POOR QUaI2^, 
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(ill) Precession and nutation models relating the instantaneous 
spin axis direction to a space-fixed frame of reference. 
Errors in these, however, are small (Ref. 11). 

All of the effects discussed above in 1) and 2) have diurnal 
signatures. At any instant their combined effects can be represented 
by an equivalent set of station locations which would produce the 
same effect. They are therefore referred to as equivalent 
station location errors (ESLE) . 

3) Locating the earth in a target frame of reference: Errors 

in this for an interplanetary target arise' from transformations 
between frames of reference (a) G and B (b) B and H and (c) H and P 
(Figure 1-1) . 


(a) The error in the earth ephemeris relative to the earth- 
moon bary center stems primarily from uncertainties in the 
precise masses of the earth and the moon; the lunar ephemeris 
relative to earth is very accurately known. If r , r are the 
position vectors to the barycenter and the earth respectively, 
in some inertial frame of reference, then by definition of 
the barycenter. 


— B — E y. 


M ’'■e 


-EM 


(1-3-2) 


where 


r^ = earth centered lunar position vector 
= gravitational constant of the moon 

y = gravitational constant of the earth . 

K ' ' 
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Letting 


•'e 


(1-3-3) 


-E -B ^ n -EM ■ 


(1-3-4) 


Taking variations on this equation, we obtain 


r = r “ — T 
-E ^2 -EM 


(1-3-5) 


where r^ has been held fixed since it is very accurately determined. 
Ignoring 6r for the present, since the effect of that is dealt with 


in effects (b) and (c) , this yields 


j: — EM ~ 

6r 6 t] 

Hi 2 

n 


(1-3-6) 


Since r (t) has a 28 day period, the error, <Sr„(t) , in r (t) 
will vary in a sinusoidal fashion. 

Using values 

n 80 , 


o ~ 4 X 10 

n 


= 4 X 10 km 


-we obtain. 


.02 km . 


This is'ao small it would appear this error source is insignificant. 
However, differentiating Equation (1-3-6) we obtain 

^-E 2 -EM ■ 

T1 

Using an r(0^ centripetal acceleration, with 


to - radians /day 
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this yields 


*15 X 10 km/sec^ * 

This is not insignificant relative to uncertainties in the 
unmodelled accelerations acting directly on the spacecraft, e.g*, 
solar pressure. 

(b) and (c) J The second and the third of the sources of error 
arise from uncertainty in the Heliocentric ephemeris of the 
Earth-Moon barycenter and that of the target planet respectively* 
^ In Qur discussion later these are referred to as errors in the 
planetary ephemeris. 

Glassification of Error Sources; Effect on Navigation 

It should be emphasized that the discussion above is a simplified 
picture of the errors. In practice correlations exist between 
various error sources. For instance errors in the heliocentric 
ephemeris of the~target planet would be correlated with errors in 
that of the earth-moon barycenter. This is because the source of 
the ephemeris information is largely, if not entirely, observations 
taken from the earth* However, for ease in analyzing their effects, 
we group all the errors other than data measurement noise and random 
components of error, into the following classes according to their 
effect on navigation accuracy, 

1) Equivalent station location errors (ESLEs) 

2) Ephemeris errors 

3) Unmodelled spacecraft accelerations. 
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These constitute the major error sources for space navigation 
based on earth-based radio data. We now discuss the effects of 
each class of errors keeping the discussion as simple as possible. 
The analysis, thougli greatly simplified shows some of the salient 
features of radio navigation limitations. 


i) Effect of ESLEs 


The effect of equivalent station location errors can be well 

represented through the Hamilton-Melbourne model. Following the 

terminology of Section 1.2, let Ar , AX and AZ be the errors in 

s s 

the equivalent station location coordinates. The variation of 

equation (1-2-10) is given by 

Ap = Af + o)sin(a)t + X - a){Ar cos6 - A5*r sin6} 

s - s 

+ mr cos6 cos(mt + X - a) {AX - Aa} . (1-3-7) 
s 

The last term in Equation (1-3-7) indicates that an error in longitude 

will map directly into an error in spacecraft right-ascension, 

Aa - AX , (1-3-8) 

yielding a position error in this direction given by 

0 ^ = r cos6 0^ . (1-3-9) 

a 

Similarly, the second term of Equation (1-3-7) indicates that an 
error in station-radius produces an error in the declination of the 
probe 


AS = 


1 s 

tan6 r 

s 


( 1 - 3 - 10 ) 


which yields a position error in this direction given by 

‘ _ s 

£ ^ tan r 


(1-3-11) 
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Equations (1-3-9) and (1-3-11) , show the linear dependence of the 
navigation accuracy for these components on the distance of the 
spacecraft from the earth. The latter equation also shox^s that the 
position uncertainty in this direction becomes very large when the 
geocentric declination approaches zero. 

2) Effect of Ephemeris Errors 

These are perhaps the most directly understood error sources, 
since the target relative spacecraft vector is the difference between 
the vectors to the target and to the spacecraft respectively. We 
note that (i) when the target is the Moon, these errors essentially 
vanish since the geocentric lunar ephemeris is known very accurately 
and (ii) the target ephemeris error increases with distance of 
the target from Earth since the ephemeris is established basically 
through observations from Earth. 

3) Effect of Unmodelled Accelerations 

The basic reason for these effects is the nature of the radio data. 
The "d" term of the six parameter model mentioned in Section 1.2 is 
useful in explaining the effects of the spacecraft acceleration 
errors. Since 

d = a + r(d^ cos^S + 5^) (1-3-12) 

g 

and it is the only term containing r, any unmodelled acceleration 

will cause an error in the determination of the range from doppler 

-12 2 

tracking data. Errors in acceleration of the order of 10 km/sec 
can cause a range error of hundreds of kilometers. 

When the range data type, p , is included in the processing^ 
range errors are essentially eliminated. In this case any errors in 
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acceleration are thrown into a and 6, i.e., the velocities perpendicular 
to the line of sight. This does not cause errors in position at 
the time of the range measurement but the trajectory does not propagate 
well. More details on this can be found in Reference 3, 

Summary 

In summary then 

(i) errors in radio navigation increase as the spacecraft 
distance from earth increases, owing to ESLE effects 
and target ephemeris effects. 

(ii) particularly large errors can be induced fox low 
geocentric declination geometries 

(iii) lunar ephemerides are much better known than the planetary 
ephemerides 

(iv) the unmodelled acceleration errors are peculiar to the 
nature of the radio data. 
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2.4 Optical Navigation 

Optical Navigation, as used in this dissertation, is the use 
of observations of celestial bodies viewed by the television 
camera on-board the spacecraft to supplement tlie earth- based 
radio data, in order to determine the orhit of the spacecraft. 

In the previous sections we have seen that there are certain 
limitations in the use of radio data for space navigation. As 
pointed out, these limitations can become significant for trans- 
lunar space exploration. The use of optical data from the TV 
cameras can alleviate some of the difficulties, giving Improved 
navigation accuracies. As we shall see in the next chapter optical 
data and radio data tend to complement each other for interplanetary 
navigation; we therefore stress that the term ’’optical navigation"*’ 
is used in the context of the combination of radio and optical data 
and not merely optical data by itself. 

In our discussion we will restrict ourselves to only the 
interplanetary portion of the mission where the spacecraft is 
essentially in an elliptical path around the sun on a hyperbolic 
trajectory relative to the target planets As we shall see the 
primary benefit in the use of optical data is gained when the 
spacecraft reaches reasonably -close to the planet -- how close 
depends on the characteristics of the instrumentation available and 
on the particular planet in question. Thus the utility of optical 
navigation is primarily during the phase of .the spacecraft’s 
’’approach” towards the planet; navigation accuracies in this phase 
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are very Important and can be quite critical. Since the laws of 
kinematics do not generally distinguish between the direction of 
motion and since the hyperbolic trajectory of one body relative to 
another is symmetric between the approach and departure legs, the 
method is equally applicable for the departure leg if there is 
any for the mission under consideration. 

The idea of using imaging data from the television cameras 
turns out to be quite effective as we shall see in the following 
chapters, and yet is simple to execute. The simplicity arises 
from the fact that no additional instrumentation is required other 
than that which normally exists on most interplanetary spacecraft. 
Instrumentation 

The instrumentation that is normally aboard interplanetary 
spacecraft and relevant to optical navigation consists of (i) the 
television cameras, (ii) a platform, moveable relative to the 
spacecraft, upon which the television cameras are mounted along with 
other science instruments and (iii) an attitude control system to 
orient the spacecraft relative to inertial space. 

The four photovoltaic solar panels in Fig, 1-4 (Ref, 12) which pro- 
vide the energy for spacecraft functions, need to be oriented facing the 
sun. This is accomplished through the use of the sun sensor, the 
line of sight of which is perpendicular to the solar panels and 
establishes one basic orientation axis — the spacecraft-fixed 
roll axis. This sensor, nominally pointed towards the sun, provides 
pitch and yaw axis control signals. The second reference direction 
required for the attitude of the spacecraft to be fixed is provided 
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Figure 1-4. Spacecraft Schematic 



by a roll axis control signal from the star tracker* The star 
used is normally Canopus; any bright star would do in general for 
this purpose so long as its direction is sufficiently non-parallel 
to the sun direction reference axis. The spacecraft orientation 
is maintained by using the pitchy yaw and roll control signals to 
actuate the attitude control reaction gas jets. 

The television cameras along with other science instruments 
are mounted upon the scan platform. This is a support structure 
moveable relative to the spacecraft with two degrees of freedom 
by drive signals to the scan actuators. Platform position is 
detected by potentiometers mounted on the actuator shaft, yielding 
measurements of the two axes gimbal angles* 

Observation Methods 

The celestial bodies that the television cameras can view 
include stars, planets and satellites of planets. Figure 1-5 
shows the three observation methods that we shall investigate. 

In the first method a sequence of pictures of the planet is taken* 
Orientation information for the television is derived using signals 
from the sun sensor, Canopus tracker and potentiometers for scan 
platform gimbal angles. 

The second method uses the fact that the inertial directions 
to the stars are known very accurately — better than 1 arc-second. 
Moreover, since they are extremely distant from the solar system 
the -direction to a star is -unaffected by change of position within 
the solar system. They -can therefore provide accurate reference 
directions , in ;Lnertial space. Thus in this method a sequence of 
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pictures is taken, viewing the natural satellites of the target 
planet against a star background. The star background serves as 
an accurate reference for determining the camera pointing direction 
while the natural satellite images yield the navigation information. 

The second method would not normally be used to image the 
planet against a star background because the planet brightness is 
much larger than that of the stars. Thus camera exposure times 
that would be adequate for imaging a planet would typically be too 
short to simultaneously image stars in the same picture. If however 
two narrow angle television cameras are available on the spacecraft 
we may employ the third method where one of the cameras is used to 
image the planet as in method (i), while the second is used with 
longer exposure time to photograph stars. The star field, as in 
method (ii) , provides significantly more accurate pointing information 
than is obtainable from the scan platform and attitude control 
data as in method (i) . 
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1.5 Coordinate Systems and Transformations 


The major coordinate systems to he used for optical navigation 
measurements are described in this section. The associated coordinate 
transformations are also derived here for use later in Chapter Three* 
Notation: We define a rotation R such that 

% = . (1-5-1) 


i.e., it is a transformation which carries a vector with components 

given in system J by k into a vector In system K with components 

J 

given by Similarly, 

-5 “ ’ (1-5-2) 

where 

^ ’ (1-5-3) 

We note that 


(i) Since these are orthonormal matrices, we have 




(1-5-4) 


and 


(ii) the associative rule holds, 

^JK ®KL °° (1-5-5) 

Inertial xyz Coordinate System 

One of the basic coordinate systems used in celestial mechanics, 
and the primary system for a space- fixed frame of reference in this 


dissertation is the 1950.0 Earth Mean Equator and Equinox coordinate 
system (abbreviated to 1950,0 EME) • The fundamental plane in this 
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system is the plane of the earth mean equator as of Jan* 1 Oh 1950 
and the reference direction . is toward the vernal equinox, which is 
the point of intersection of the plane of the sun^s apparent motion 
about the earth where the sun crosses the equator from south to north. 
The X axis is the direction of the vernal equinox, the z axis is 
normal to the fundamental plane and positive toward the north, and 
the y axis completes a right handed system. 

Celestial Coordinate System, ABC 

As mentioned in Section 1.4, the orientation of the spacecraft 
during its flight is maintained through the use of the attitude 
control sensors which are nominally pointed towards the sun and 
Canopus respectively. The directions to these celestial bodies 
are used to establish the ABC Celestial coordinate system, as shown 
in Figure 1-6. The three axes are defined as follows: 

^ = unit vector to sun 

^ - un it ve ctor to reference star 

A A A.iA At 
b = _c X _s/ X ^1 

AAA 

^ ^ X jc 

p ^ - cone angle of star = cos ^ (c • ^ 

where the "cone angle" of a vector direction p is defined as the 
angle that the vector makes witd the positive £ direction (see Figure 
1-7). The "clock angle" of the vector p is defined as the angle 
which the £ x p vector makes with the positive ^ direction (see Figure 
1-7). 

The rotation matrix to transform a vector from the Inertial 

s 

coordinate system xyz to the Celestial coordinate system ABC is 
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Figure 1-6. 
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therefore given by 


or 



(1-5-6) 


(1-5-7) 


where a^, and are the components of the vectors _a, and ^ 
expressed in the inertial xyz system. 

Spacecraft-Fixed Coordinate Systems 

There are basically two spacecraft structure related coordinate 
systems. They are (i) the sun- star sensor coordinate system, S, and 
(ii) the spacecraft attitude control sensor related coordinate 
system, X. 


(i) Corresponding to the celestially referenced ABC Celestial 
coordinate system, it is useful to define a Spacecraft- 
fixed coordinate system Sg which would coincide with 
the ABC system when the pitch, yaw, and roll attitude 
control system signals are zero. It is emphasized that 
this is a system of coordinates physically related to the 
instruments — the sun sensor and the Canopus tracHer -- in 
the spacecraft. When these sensors are pointed exactly 


towards the sun and Canopus, respectively, the system is 
coincident with the G system and the pitch, yaw, roll angles 
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are zero. 


(ii) XYZ is the spacecraft pitch, yaw, and roll control axis 
coordinate system and is defined relative to the fixed 
spacecraft configuration as shown in Figure 1-4, The 
X and Y axes are aligned along the solar panels and the Z 
axis is antiparallel to the axis of the spacecraft- 
fixed sun- star sensor coordinate system. 

It is convenient at this point to introduce a new celestial 
coordinate system X^Y^Z^ such that if the pitch, yaw, and roll control 
signals are all zero the X Y Z system would be coincident with the 
XYZ system. This implies that the X^Y^Z^ system has the same 
relationship to the ABC system as the XYZ spacecraft-fixed system 


has to the (also spacecraft- fixed) system. 

As shown in Figure 1-8, a rotation R _ through the angle a 

CXq! X 

about the <3 axis (axis # 3) followed by a 180° rotation about the X 


axis (axis # 1), E. , takes the ABC coordinate system to the 

XYZ system, 
o c o 

The rotation from ABC to X Y Z is therefore given by 

o o o o V 




(1-5-8) 
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Figure 1-8. Space craft/Scan Platform Coordinate 
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To rotate from the X Y Z system to the XYZ system requires 

o o o 

rotating about the pitch (X) axis, yaw (Y) axis, and roll (Z) axis 
through the angles 9 , 0 > 9 respectively. The angles 0,0, and 

p y ^ p y 

9^ are determined from the control system signal voltages from the 

sun sensor (pitch and y«w signals) and the star tracker (roll signal) . 

The three rotation matrices are, respectively, 

0 0 

cos 6 sin 0 

P • P 

-sin 0 cos 6 

P P 






The transformation ^ is then given by combining Eqs. (1-5-9), 

XqX 

(1-5-10) and (1-5-11), 
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COS0 cosS fsin8 cos0 +cos© sinQ sin© 
r r p r y p 

-sinQ cos9 * (cos0 cos0 -sin0 sia0 sin0 ) 

r y, '■ r p r y p'' 


sin0 ;-cos0 sin.0 

y I y P 


sin© sin0 -cos0 sin0 cos0 


COS0 sin0 
r p 


+sin0 sin© sin© 

r . y P 


COS0 COS0 

y P 


(1-5-12) 

The rotation transforming from the spacecraft XYZ system to 
the Spacecraft system, Ry^ is just the inverse of the 

transformation from the Celestial ABC system to the XYZ 
system, R obtained in Eq. (1-5-8). Using Eq. (1-5-4) we have 


"CX, 


^SX “ \s 


(1-5-13) 


We observe from Eq. (1-5-8) that is a symmetric matrix for 

the particular definition of the XYZ system here, therefore 


\s " 


cx, 


0 


This would not be true in general. Here, then 




cos o; 


sin O' 


sin O' -cos a 
X X 


0 

0 


(1-5-14) 


L 0 0-1 

The composite rotation from the Celestial ABC system to the 
Spacecraft Sg system is then given by 


R, 


CS 


®XS ®X X “CX 


o o 


(1-5-15) 
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which is equivalent to writing 


^CS ®X X ^SX 

o 

r 

^ ®XS \ X ®xs 
o 


(1-5- 15 a) 


" ®XS ®X X ®XS 

O 

Scan platform Coordinate System 

' This Platform MNL coordinate system is defined relative to the 
spacecraft coordinate system through the^ clock angle oTp and the cone 
angle pp (as shown in Figure 1-8) of the Scan Platform pointing 
vector The L axis is in the direction of the M. axis is in the 
direction of increasing cone angle pp and the N axis is in the 
direction of increasing clock angle 

The transformation from the Spacecraft to Platform coordinates, 
R , can be written as 

^SP " ^ P ^SP (1-5-16) 

O' ‘ or 


where 


R = transformation fronL spacecraft coordinates to P 
Sir ^ 

a 

coordinates, through a clock angle rotation 

E_ _ = transformation from P coordinates to Platform 
Tf P O' 

a 

coordinates through a cone angle rotation Pp. 


These are given by 



’ cos Qfp sin Ofp 0 

-sin ofp cos Qfp 0 

0 0 1 


(1-5-17) 
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cos 0 - sin Pp 

0 10 

siuPp 0 cos'pp 


( 1 - 5 - 18 ) 


The angles and are determined from clock and cone gimbal 
potientiometer signals. 

Television Camera Coordinate System 

The Television Camera MNL coordinate system definition relative 
to the Platform MNL coordinates is also shown in Figure 1-8; 
arriving at the television camera coordinates involves three successive 
rotations through angles tjr, x> these must be taken about the 

positive N axis/ the negative M axis and the positive ^ axis, 
respectively. The rationale for defining the coordinates in this 
fashion is the following: 

positive ^ => increasing p ; 


positive ^ -> increasing or (for ^mall ^ angles) ; 
positive 0 ) positive rotation about X, the look direction; 
-and therefore the angles are referred to as the cone, cross-cone, 
and rotation offsets. 

The transformation from Platform coordinates to Television 
coordinates is given by 

^ ^i2 (1-5-19) 


where 



“ 

cos ijr 0 

0 1 

■ sin i|f 0 


-sin ijf 

0 

cos ^ 


( 1 - 5 - 20 ) 
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^^2 


and 




10 0 
0 cos X -sin X 

0 ' sin X X 


cos 0) sin 0 ) 0 

•sin to cos (0 0 

0 0 1 


(1-5-21) 


(1-5-22) 


Substituting (1-5-20) to (1-5-22) into (1-5-19) we obtain 

r 

cosijrcosa) - sin\Jrsinxsin{JD cos^sinu) -sini|fcosa) - sinxcosiijsina) 

= -cos^sino) - • sinijfsinxcosa) cosycoso) sin\}isinu) - cosijfsinxcosa) 
cosxsin^ siux cosxcos\}i 

(1-5-23) 

Composite Transformations 

The overall transformation from inertial to television coordinates 
is given by 

■ VsAs^ic 


where R is given by 

Lf o 


®cs "" ^sx®x x^sx 


(1-5-25) 


and where R^p and Rp^ are given by Eqs* (1-5-16) and (1-5-19), 
respectively. 

Through the use of the sensor angles 0^, 0^, 6^, and these 
relationships then yield the transformations to go from an inertial 
frame of reference to a television frame of reference. 
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1.6 Historical Perspective and Scope of the Dissertation 


In this section we give an outline of the historical background of 
optical navigation, followed by a synopsis and a brief description of 
the scope of this dissertation. 

As outlined in Section 1.1 space navigation began with the use of 
radio data for the early earth orbiters and for* subsequent lunar 
exploration. Radio navigation methods developed through the years into 
use for interplanetary missions. The area of optical navigation' is a 
relatively new one, pioneered over the last few years. 

A simplified investigation was made in Reference 13 to make a pre- 
liminary assessment of the use of an onboard navigation instrument. 
After the proposal for exploiting the "Grand Tour" mission (Ref* 14) 
opportunities for multiple encounters with the outer planets, there was 
interest (Ref. 15) in an on-board capability for interplanetary naviga- 
tion exploring the use of various instruments such as the sextant 
(Ref. 16) and the planet sensor (Refs. 17 and 18). The use of an 
onboard television camera to view the planet (method (i)) was suggested 
in Ref. 18 and a feasibility demonstration, made in Ref. 19. Refer- 
ence 20 proposed the use of satellite-star data (method (ii)) using a 
television camera for navigation of the Grand Tour missions. Further 
studies on the effects of spacecraft acquired optical measurements upon 
the orbit determination of two Grand Tour trajectories were carried out 
in Ref. 21. The particular camera configuration in the Viking Mission 
made the two-camera planet-star method — method (iii) possible, and 
Ref* 22 made a preliminary relative evaluation of the three methods. 

An in-depth application of the satellite-star and planet limb methods 
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was studied in Refs.. 23 and 24 respectively. With interest aroused in 
optical navigation, the emergence of the Mariner-Jupiter Saturn mission 
led investigators (Ref. 25) to study the requirements for the quality 
of ra4ipi!ietric and .optical -data. 

In the following. Chapter 2 deals with analysis and modelling 
issues related to the nature of the optical observation methods. 
Beginning with the definition of the data type, we continue with an 
examination of light time and stellar aberration effects, and their 
interaction. The errors affecting the optical data are classified 
according to their effects and an analysis is made of the information 
content of the data. 

Chapter 3 describes the overall functions involved in spacecraft 
navigation. These functions include dynamic and- measurement system 
modelling, data sequence design, measurement extraction, model esti- 
mation and orbit determination. The relevant details as applicable to 
optical navigation, for each of these functions, are developed here. 

In Chapter 4 we present the application of the satellite-star and 
the planet-limb (one camera) methods to data obtained from the 
Mariner 9 mission. This includes analysis of optical navigation sensi- 
tivities to significant parameters, with ^n emphasis on sensitivity to 
the amount of radio tracking data used. 

Chapter 5 uses the Viking Mission to examine in detail navigation 
accuracies obtained with the use, respectivity, of radio, plus optical 
data and of optical data only. This is done for the satellite star and . 
the .planet-star two camera methods. The sensitivity to the timing of 
the midcourse maneuver is .examined. In addition, a general method is 
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developed for evaluating the consistency between any two estimates. An 
application of this method, with interesting results, has been made using 
radio and optical estimates from simulated data. Finally we develop a 
parametric probability analysis to evaluate overall navigation system 
performance as a function of optical system reliabilities. 


41 



CHAPTER 2 


OPTICAL DATA ANALYSIS AND MdDELLING 

INTRODUCTION 

The optical data types concept was introduced in the first chapter. 
In this chapter, we shall first -formalize the concept to define the data 
type and observation model with respect to the coordinate system in 
which they are measured in Section 2.1. 

In Section 2.2 we shall examine the basic geometrical relationships 
between the S/C and the celestial bodies involved in the interplanetary 
navigation problem and how this relates to the optical navigation data. 
This is followed by a discussion of the effects of the finite velocity 
of light on the data recorded. These effects consist of an interaction 
of light-time and stellar aberration. After these are factored in, the 
definition of the data type is complete. 

The next section. Section 2.3, deals with the sources of error 
contributing to the inaccuracy of the optical navigation process. The 
magnitude of the errors are presented along with a discussion of the 
reasons for these. 

Section 2.4 examines the basic information content in the data as 
it relates to the determination of the orbit of the spacecraft. This 
information is studied in the context of the desired navigation data we 
are seeking. From the analysis presented here, restricted to two dimen- 
sions so as not to confuse the basic issues involved, some conclusions 
regarding the observability of the system have been drawn. 
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2.1 Definition of Data Types 

The images viewed by the television camera for optical 
navigation can be a planet, one or more of its satellites and stars. 

The optical data types are defined to be the centers of the 
images of the, objects in TV coordinates as viewed in the picture. 
These are defined relative to the scanning process as illustrated 
in Fig. 2- la. 

After the camera is shuttered, the light falling on the photo- 
sensitive vidicon surface causes electron-hole pairs to be created 
due to the photoelectric effect, leaving a net positive charge on 
the surface of the vidicon. The amount of charge is a function of 
the number of photons striking the surface. This charge is then 
erased by a read-out beam which electronically scans the surface of 
the vidicon in a matrix of MxN digitized samples. The sampling is 
accomplished with the use of pulses of electrons falling over a 
certain "sample ar ea" w hich shifts from pulse to pulse due to a 
horizontally imposed ramp bias to scan along a line (see Fig. 2-l(b). 
After the end of the line the value of the ramp bias restarts at zero 
and simultaneously the vertical bias shifts the sample area in the 
vertical direction to begin the next sampling line. Each sample 
area is called a picture element or "pixel." The magnitude of 
the current read by each pulse gives the amount of charge for that 
sample and this can be related in turn to the intensity of light 
falling upon the vidicon surface at^ >that pixel. The resulting 
charge read for various samples is digitized in n binary levels which 
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provide up to 2^ digitized ”gray levels’^ o£ intensity called the 
Data lJumber (DN) ♦ 

Coordinates 

The scanning process for a given picture begins as shown in 
Fig* 2- la near the upper left-hand corner. The sample there is 
designated Pixel number = 1, Line number = 1. The scanning 
directions — increasing pixel numbers along a scan line and 
increasing line number — provide the two reference directions, the 
p and the 1 coordinate directions. These along with the origin 
at the upper left-hand corner provide the (p, 1) coordinate system. 

Fig. 2-la also shows the x, y, z coordinate system the origin 
of which is on the surface of the vidicon at the central line and 
pixel The x and y directions are parallel to the p and 1 

directions, respectively, while z is along the line of sight (LOS) 
of the camera and completes the orthogonal x, y, z right-hand system. 
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2.2 Basic Geometry; Stellar Aberration and Light-Time Effects 
The optical data consists of the image of the object in TV 
coordinates. The scanner provides a window into space enclosing 
a solid angle which is the field-of-view (FOV) of the camera. 

The center of this solid angle (the camera line- of- sight or the camera 
look- direct ion) provides the pointing of the camera relative to some 
reference direction in space. The location of an image relative to the 
center of this solid angle therefore yields the information regarding 
direction to the object in space when transformed through the appropri- 
ate transformation of coordinates. We shall examine these in detail in 
Chapter Three; here we examine the overall geometry of the process of 
imaging an object for a single observation. 

2-2 shows the relative vectors between the sun, spacecraft, 
the object being viewed. If 
A Heliocentric spacecraft vector, 

A Plane tocentric spacecraft vector, 

A Planetocentric object vector, 

A Heliocentric planet vector, 

t = £ - ^ (2-2-1) 

A spacecraf t--centered object vector 

Z = £ + ^ (2-2-2) 

A heliocentric object vector. 


Figure 
planet, and 


X 


and 

then 


J1 

£ 


where 


and 


where 
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If the passage of light were instantaneous, the image would 
yield the direction where 



unit vector in the direction of 


t. 


However the situation is complicated by two distinct but intimately 
related effects* These are (i) light-time correction and (ii) stellar 
aberration. Both of them are due to the fact that the velocity of 
light is not infinite. The first can be considered to , depend on 
the motion of the object and the latter on the motion of the 
observer; their intimate connection stems from the fact that there 
is no absolute frame of reference. 

Light-Time Correction 

Light-time correction as the name implies is the motion of 
the object during, the time taken for the light to traverse the 
distance between the body and the observer^ Let 
vector position, at time t,, of the spacecraft relative to any 

I 

origin in some <non-rotating)^ frame of reference F; let 

be the vector position at time t^ of the object being viewed in 

the same frame of reference F. 


Then the light- time ,equation ignoring relativistic effects can 


he expressed as 
T = t. 

3 



u. (t .) - w. (t 0 

— ‘j j —i 1 


( 2 ^ 2 - 3 ) 


where the light travels from w. to u, in frame F, and t is the 

— 3 

light-time. It should be noted that in general t depends on the 
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frame of reference. The direction from the spacecraft to the object 
before applying the light- time correction, in the reference frame F, 


is 




- u. 



(2-2-4) 


The light-time equation can be solved by iterating the- equation. 

An adequate approximate expression can be derived by assuming uniform 

motion. The component of the object's velocity in the direction of 

increasing range is given by 

V = w. t , • (2-2-5) 

r T — 

where w. is the velocity of the object at time t.. Therefore, 
c(t. - t,) » |w. - u. I - t}(t - t ), , 


or 


|w, - u, I 

t. - t. 

3 X 


c u + 


' W, 

n 


(2-2-6) 




The effect of the factor irr the denominator is of the order of 

V 

sine a the value af x itself is of the order of — . We shall see 

2 "" 

V 

shortly that effects' of the order of can be ignored; in fact that 

c 

is the basis upon which we can afford to neglect relativistic effects. 
Thus the value of the light- time t can be approximated by 


T = t . - » 


[w. - U. i 
Zll 


(2-2-7) 


Stellar Aberration Correction 


Stellar aberration is the difference between the direction 
of the incoming light and the apparent direction from which the 
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observer sees it coming. This effect is due to the velocity of the 
observer not being negligible with respect to the velocity of a 
photon, and was originally discovered in 1725 by James Bradley. 

(He was attempting to find the distance to a star based on observing 
the parallax due to the position of the earth in its heliocentric 
orbit. The effect he observed was rotated approximately 90® in phase 
from the result he expected, i.e., it was in phase with the velocity of 
the earth.) Ignoring relativistic effects again, the effect is easily 


understood using the particle theory of light in terms of the velocity 
of the observer relative to the incoming photon of light. If the S/C- 
object direction is , then referring to the velocity triangle in 
Fig. 2-3 a we have, 

^ > (2-2-8) 


where is the velocity of the spacecraft at time t.. Thus the 
^ 3 

observed direction to the object before applying the stellar aberration 
correction, i.e., the actual observed direction, is given by 



(2-2-9) 


Coordinate Frame 

As noted earlier, the results obtained are independent of the 
coordinate frame used to carry out the computations. To illustrate 
this and to get a physical feeling for the combined magnitude of 
the effects, consider the situation shown in Figure 2-3b. In the 
figure represents the trajectory of the observer- spacecraft 


50 




based TV camera — and represents the trajectory of the object 
being viewed. In general these trajectories will lie in different 
planes , 

The light- time t is computed as 



c 


Now if ^ is the vector angle correction in magnitude and direction 
(follox^ring the right-hand rule) due to light- time, then we can write 

/I 

_t^ X (w. - w.) 

lu. • ( 2 - 2 - 10 ) 


Now 


w. - W. = W.T + OCt^) 

-3 -I -3 

2 

= w. — — * + 0 (t ) 

-j c 


Therefore, 


e 

T 


t * X W. 




+ 0(t ) 


( 2 - 2 - 11 ) 


or,. 

jt* X w. 

0 ^ = — 1 + O(t^) . (2-2-12) 

If is the vector angle correction due to stellar aberration, 
then from Fig. 2-3b 


u. X 

^SA = • ' (2-2-13) 

The combination of the two effects is defined as "planetary 
aberration." Letting 0p denote planetary aberration, adding the 
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two effects we obtain 


£’ X . (w. - u.) 


0p = 0^ +lsA- 


We note that (i) this result is independent of the frame of reference 

that we choose to work in, even though each of 0 and ^ are not. 

This points up the close connection between the two effects, (ii) each 
i- 1 V 

of the corrections is of order — and we are neglecting quantities of 

order “-2- This is justified since values for ^ are typically^ 10 
2 c c 

_8 

so 10 , and the pixel resolution of typical cameras on board 

- 

interplanetary spacecraft ranges from 10 y radians to 50 y radians. 

The situation as outlined above is^ complicated by the fact that 
since the stars are at enormous distances from the solar system, 
catalogs for their ephemeris are not made for proper time, i.e., it 


is not their actual location at time t that we record but rather their 


apparent location in a heliocentric frame of reference. The procedure 
followed when doing this is to correct for the stellar aberration 
caused by the heliocentric motion of the observer on the earth, but no 
account is taken of the light- time. This is to avoid complications in 
the star catalogs pertaining to the distance of the stars from and 
their motion relative to the solar system. 

Thus for images of stars the light- time correction does not apply 
leaving only , the stellar aberration (hence the name). This means then 
that the stellar aberration for the star images must be applied in 
the heliocentric frame of reference. 

Differential Aberration 

A further simplification results in the stellar aberration effect, 
when stars are used as the reference for the determination of camera 
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orientation, as in the satellite-star and planet-star (two camera) 

methods. What is important here is just the relative magnitude of the 

stellar aberration effects between the star and non-star images. This 

is caused solely by the difference in the unit vector directions to the 

two objects and is termed differential aberration. 

To evaluate the approximate value of the differential aberration 

effect, consider tiro objects (one of which is a star) in directions 

_t^j^ and ^2 respectively; the stellar aberration effects will then be 

given by — (u. x ^ ) and “(u. x respectively. Letting a be the angle 
c 3 J- ^ J ^ 

between\t^, and writing 


/\ A ^ ^ 

• -£2 - 

then for small values of a we wij i ) ave 

a | a ^[ 

The differential aberration- effj-.rt is then given by 

= — (u. X At) 

SA c — j — 


(2-2-15) 


(2-2-16) 


(2-2-17) 


The maximum value of this, for a given u^ and a will occur when A^ is 


normal to u . . This maximum will be 
/MAX c 




(2-2-18) 


Since the two objects must be viewed with a fixed- orientation 
camera structure the maximum possible value of a is governed by the 
field of view of the camera for the satellite-star method, or the 
combined fields-of-view of the two cameras for the planet-star method. 
Thus using 


a '^0.04 radians 
max 


(2-2-19J 



op; ppoR 
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( 2 - 2 - 20 ) 


in Eq- (2-2-18), we obtain, 

(“sa) - 

\ /max 

- 4/1 radians 

This effect is therefore small in comparison with the other error 
sources and can be neglected for the current state of the art-, 
Differential Planetary Aberration 

To compute the net effect we define another term^ differential 
planetary aberration, as the combination of the differential aberration 
effect above and the light- time effect for an object other than a star. 
This quantity could also be alternatively obtained as the difference 
between the stellar aberration effect on the star and the planetary 
aberration effect on the other object. 
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2.3 Optical Data Error Modelling 

Figure 2-4 illustrates the geometry involved in the use of 
optical data for spacecraft navigation. Conceptually the following 
frames of xeference are successively involved in the process of 
obtaining the planet relative spacecraft state vector: 

S: The spacecraft and its coordinate system 

TV: Television camera frame of reference 

N: Natural satellite frame of reference 

F: Target centered frame of reference. 

The N coordinate frame is only implicitly involved and that only for 
satellite observations. For planet observations the referencing is 
directly between the television frame of reference and the planet. 

We notice that there are fewer transformations between frames 
of reference in the use of optical data compared with those for 
radio navigation. This fact means little without consideration of * 
(i) the accuracies obtainable at each step in the process and (ii) the 
basic information content inherent in the data i^e., the observability 
of the navigation quantities from the data — and the effect of the 
errors on orbit determination. These two issues are qualitatively 
discussed in this and the following .sections, respectively. The 
problem being highly Jion- linear, detailed quantitative studies are 
only possible numerically in specific applications ^ we do that in 
Chapters 4 and 5. 

In this section we list the major error sources resulting from 
the use' of optical data for navigation. The television camera is 
taken to he the primary sensor for all of .the optical data types, the 
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VECTORS REFERENCED 
TO THE S/C FRAME 

VECTORS REFERENCED 
TO THE TV frame 


S - SPACECRAFT COORDINATE SYSTEM 
T - TV CAMERA COORDINATE SYSTEM 
N - NATURAL SATELLITE COORDINATE FRAME 
. P - planet COORDINATE FRAME 


Figure 2-4. Coordinate Frames of Reference for Optigal Navigation 



other sensors (mentioned in Section 1*1) providing the means to 
relate frames of reference* 

Paralleling the categorization of radio navigation errors, we 
can categorize the errors into the following three types: 

(1) Imaging instrument measurement imprecision errors 
(ii) Spacecraft acceleration model errors 
(iii) Errors in relating frames of reference* 

Imaging Instrument Measurement Imprecision Errors 

These errors arise within the television from effects related* 
to the optical or the electronic part of the imaging system or to 
the photosensitive surface of the vidicon* Some of these are 
invariant characteristics of the television camera that affect the 
entire picture frame* These are independent of the image* Others 
are a function of the specific image being viewed* Unlike instrument 
related errors for radio data, these errors for optical data are not 
small; they can be quite significant in their effect on navigation 
accuracies • 

1) Image- independent errors: -These consist of three types of 
errors listed below: 

(a) Uncertainty in the television system parameter values such 
as, (i) the focal length of the optical system (ii) the scale 
factor and non-orthogonality terms of the transformation from 
the optics system Image plane to^ the line-pixel coordinates on 
the picture frame. 

(b) Corruption of the image due to geometric distortions in 
the telescope of the instrument and the electronics of the 
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vidicon tube. These are referred to as optical distortion and 
electromagnetic distortion respectively. Sources of optical 
distortion include (Reference 26) (i) imperfections in the 

telescope lens (ii) misalignment of the lens optical axis in 
direction (with respect to the normal to the target raster) 
and position (with respect to the center of the target raster) . 

Sources of electromagnetic distortion include (i) non- 
uniformities in the magnetic deflection and electric deceleration 
fields (ii) fringe field outside the deflection region of the 
vidicon tube (iii) interaction between the focusing and deflection 
fields (iv) rotation, non-orthogonality and center displacement 
of the scan line and pixel direction deflection fields. 

(c) The television picture resolution. The image received 
through the optics causes a continuous charge distribution on 
the photosensitive surface of the vidicon. However in the 
electronic read-out process the electron scanning beam has a 
finite size which determines the size of the pixel. This forms 
the lower limit for the accuracy of the imaging system. 

2) Image-dependent errors: Much less is known or can be 

determined from calibration about these errors relative to the 
image independent errors described above. They are related to t 
image fomation and detection processes and may be caused by: 

(a) an interaction of the characteristics of the* imaging 
system with the method used for locating the center of the 
image. In addition to all of those listed above, the relevant 
characteristics include the point spread function of the vidicon 
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. which causes an effect known as "image blooming, "where the 
charge at a pixel is spread over neighboring pixels. This 
effect increases with the brightness of the image and leads to 
an error in locating the image center. 

(b) an interaction between the readout scanning beam with the 
point spread function and with the remaining unerased charge on 
the vidicon- This effect, known as "beam-bending," is small for 
small images such as those of stars but may be quite significant 
for large images such as those of planets. 

(c) smearing effect of the image caused by motion of the camera 
during the exposure time of the optical shutter. The motion is 
due to the limit cycle of the S/C attitude control system. This 
effect would increase with increasing exposure time. It would 
therefore usually be largest when imaging dim objects such as 
stars since a large exposure time would be designed in that 
case to insure a detectable signal on the vidicon*. 

Spacecraft Acceleration Model Errors 

Errors in this category consist of' only the gravitational 
acceleration errors applicable for radio navigation. These affect 
the motion of the spacecraft and hence the propagation of the 
traj;ectory. The non-gravitational acceleration errors do not affect 
optical data; their effect is peculiar to the nature of the radio 
observations as explained in Section 1.3. 

Errors in Relating Frames of Reference 

Referring to Figure 2-U the frames of reference involved in 
relating the target P and the spacecraft S are, in order, 
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S - TV - N - P 
• (13 (2) (3) 

where N would be omitted in the case of planet limb observations* 

The errors in this process then fall into the three .subdivisions 
indicated- The first is the error in establishing the orientation 
of the instrument relative to inertial space be" it through use of 
star images or of the attitude control system and scan platform 
sensors. The second results from errors in the model of the planet 
or satellite; the last stems from inaccuracies in the satellite 
ephemeris * 

1) Instrument Orientation Errors: Depending on whether or 

not stars are used to establish the inertial pointing direction of 
the camera, these errors can be very different in magnitude - 

(a) Pointing information from engineering data: In this mode 

the camera pointing direction is reconstructed using (i) the 
pitch, yaw and -roll sensor outputs to determine the attitude 
of the spacecraft with respect to sun-Canopus reference'd 
celestial coordinate system (ii) signals from the clock and cone ^ 
angle gimbal actuator potentiometers to determine the orientation 
of the scan platform with respect to the spacecraft -structure 
(iii) fixed offset angles to reference the television to the 
scan platform and (iv) the location of the central line and 
pixel and rotation of the TV line--pixel coordinate system relative 
to the telescope image plane axis. 

Thus errors in each of these transformations contribute to the 
error involved in reconstructing the camera pointing relative to 
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space-fixed coordinates. In addition there is data noise introduced 
into the sensor measurements for (i) and (ii) due to quantization 
of the engineering telemetry. 

(b) Using stars to determine camera pointing: This is the mode 

used for the second and third optical observation methods 
described in Section 1.4, i.e., the satellite-star measurements 
with one camera and the two-camera method for planet-star 
measurements. The accuracy of the star image detemination is 
governed by the errors in the measurement imprecision category 
described earlier. The amount of error in the camera pointing 
as estimated from the star images is a function of this accuracy, 
the number of stars and their relative location — this issue is 
further addressed in Chapter 4. The errors however turn out to 
be of the order of the TV pixel resolution. In addition to 
these errors, for the planet- star two camera method there is 
contribution from alignment errors between the two cameras and 
from spacecraft limit cycle motion during the time elapsed 
between the shuttering of the two cameras (since they are not 
shuttered simultaneously) . 

2) Object Modelling Errors: Errors in the accurate representation 

of the planet or satellite figure and illuminated profile cause 
errors in the location of the object. Thus the true location of 
the object is different from the observed location due to (i) errors 
in the mathematical model of the shape and size of the planet or 
satellite, (ii) the center of mass being offset from the center of 
the geometrical figure, (iii) errors in the illuminated profile 


OF POOR QUALIT'59 



caused by Irregularities on the surface, (iv) limb darkening effect 
due to the falling off in the brightness close to the limb, and 
(v) atmospheric effects on the surface of the planet. Tha last two 
effects may significantly degrade the accuracy of limb measurements 
or even its definition. 

3) Satellite Ephemeris Errors: In the cade of satellite 

observations, errors are caused owing to uncertainties in the ephemeris 
of the satellite relative to the target planet. There is an error 
due to the uncertainty in the planet relative satellite vector at 
an epoch and errors in the propagation of the ephemeris in time. 

The planet mass and gravity field uncertainties contribute to the 
latter. The spin axis of the planet is also a source of error when 
the satellite theory used and the planet gravity field are referenced 
to the planet equator coordinate system. 

Classification of Error Sources 

The random components of all error sources along with the pixel 
resolution can be grouped together as data noise. We regroup the 
other error sources into the following classes for ease in analyzing 
their effects on navigation accuracy. 

(i) TV pointing errors: All errors which produce a global 

shift of all the picture elements fall into this class. 

Thus all the images would be shifted a uniform amount, 
but their relative locations would be unaffected by this 
error. This is shown schematically in Fig; 2-5a. These 
are modelled as a combination of constant angular error 
and a varying time correlated component. 
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(ii) TV geometric distortion: These errors corrupt the 

relative location of all the picture elements and 
therefore of the images in a picture (shown schematically 
in Fig* 2-5b)* Ttiis is modelled as a polynomial of 
the radial distance from the center of the vidicon*- 
•(iii) Image center-finding errors: Errors that contribute to’ 

the center location of a specific image fall into this 
class. These errors are modelled as constant errors 
proportional to the image size (see Fig. 2-5c). 

(iv) Satellite ephemeris errors r Since the position of the 
spacecraft relative to Mars has to be derived using the 
position of the satellites relative to Mars, uncertainties 
in the satellite ephemeris have an important effect. 

These are modelled as errors in the orbital parameters of 
the satell-ite theory chosen to represent the motion of 
the satellite around the planet. 


65 



2.4 Optical Data Information Content Analysis 


In this section we discuss the basic information content of the 
optical data as it relates to navigating the spacecraft in the 
interplanetary approach phase of a mission. We examine the problem 
in two dimensions. Further simplifications are made to seek a 
minimal set of parameters which is sufficient to describe the 
problem. 

Figure 2-6 (a) shows the two-body motion of a spacecraft 
approaching a spherical planet located at coordinates (Xp, yp) • 

The dotted line, making ah angle relative to the x axis, is the 
approach asymptote of the hyperbolic trajectory. The eccentricity 
of the hyperbola is determined by the mass of the planet and the 
asymptotic'veloclty V^; for a massless planet the trajectory would 
be the dotted straight line at a perpendicular distance B from the 
center of the planet. The point of closest approach to the planet 
is labelled E (for "enco unter ") and the spacecraft arrives there 
at time T, the time of flight. 

With the description above the position of the spacecraft at 
any time t is at a distance H from E on the approach asyiq>tote, 
given by 

Jl(t) = T , (2-4-1) 

where 

T = T - t . (2-4-2) 

At the initial time, t == 0, the position is therefore given by 

£(0) = V T . 

CO 
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DIRECTION TO REFERENCE STAR 
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Figure 2-6. Optical Measurements during Planet Approach 



Optical Data Content 


Optical observations j either planet limb or natural satellites, 
essentially measure the angle between the target planet center of mass 
and a reference direction, e*g*, a star direction. Let T\(t) be this 
angle as shown in Figure 2-6 (b). For the purposes of this .analysis 
the error in this angular measurement represents the accumulation of 
all error sources such as center- finding errors, satellite ephemeris 
errors,, biases and camera pointing errors. These errors effectively 
result in a degradation of the observation angle. 

Since we are dealing with planet- relative angular observations 
and a planet-relative description of the spacecraft trajectory, 
the coordinates (Xp, y^) of the planet become irrelevant. The 
primary quantities of interest are B, x, and where the first 
two relate to the two dimensional position and the latter to the 
velocity- 

Let 0(t) be the angle, shown in Figure 2-6 (b), given by 
6(t) = nCt) - tjj . 

Taking variations we obtain 

60 (t) =6n(t> - S\i> . (2-4-3^) 

Now 0(t) is given in terms of the trajectory parameters by 

tane(t) 

Since the observations are taken relatively far from the planet, 

6(t) is small and 

9 I + 0(03) . .. (2-4-4) 

X/- 
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TJeglecting higher order terms we may write 


39 1 SB sa 

-6 B A 


However from Equations (2-4-1) and (2-4-2) 

6V 


Z 


^ 3t 

V T 


V r 


Substituting (2-4-6) and (2-4-3) into (2-4-5) , 


6B 50 




6T 


B 


e V 


6V 


- J(6r) - 6(f.) + — + — 


Substituting from (2-4-4) and rearranging. 


6n = <5iJ) + 


6B 
V T 

CO 


B 

V 2- 


5V - 


V T‘ 


5T 


(2-4-5) 


(2-4-6) 


(2-4-7) 

(2-4-8) 


(2-4-9) 


which is the variation on the measurement equation. 

We observe from this equation that it is not possible to 
separate an error in the direction of the approach asymtpotlc velocity 
StJ;, from a bias error in the measurement angle, Sq. 

Since optical observations involve the measurement of angles 
between the target planet and a reference direction, they suffer 
from the inability to determine the velocity (V^) of the spacecraft 
and -the accurate time of flight. To illustrate "this, consider a 
spacecraft moving on trajectory No. 1 (Fig. 2-6(c)). Let 0^ and 



represent two angular measurements of the direction between the 
reference direction (which we have for simplicity assumed to be along 
V^) and the target planet* The observable equation is 

tan 0^ = ’ ' (2-4-10) 

00 i 

where T. = T - t. and t. is the time of ith observation, 

X XX 

From Eq. (2-4-10) it is seen that the time of flight, T, can 
be determined from two perfect observations of 6. However, only the . 
ratio B/V can be determined from observations of 0 . This is 

GO 

because the observation history for any parallel trajectory with the 
same value of B/V^ (for example, trajectory No. 2) will be identical 
to the true trajectory. This could also have been observed from 
equation (2-4-9) above — i.e., since the coefficients of 6B and 6V^ 
have the same temporal behavior it would not be possible to solve 
for these two parameters separately. These parallel trajectories 
also will have the same time of flight as and will be indistinguishable 
from the true trajectory. Note also that two perfect direction 
observations determine the plane of motion. 

From Eq. (2-4-7) it is seen that even with perfect observations 
the limiting accuracy for B is determined by V^, i.e., 

6B = 6V . (2-4-11) 

y CO ' ' 

CO 

To obtain the time of flight, T, assume that two observations of 
0 are taken; then from (2-4-10), 

tan 0^ = tan 0^ , (2-4-12) 
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or 


(T - t^) tan 6 ^ = (T - t^) tan 62 , 

which yields, 

t„ tan - t. tan 8 . 

T = — . (2-4-13) 

tan - tan 0 ^ 

Even though in theory two perfect observations of 0 uniquely 
determine time of flight, in practice this quantity is rather poorly 
determined by optical data since it is extremely sensitive to 
errors in 9. This can be illustrated by examining an expression 
for the time of -flight uncertainty. Taking the variation of 
Eq. (2-4-12) ,' 

+ T^60^ == 6T202 + ^2'^®2 ' (2-4-14) 

However, from T^ = T - t^, 
using 

STj^ = 6 T 2 = 6 T, 

5 T (02 - ©2^) = Tj^S©^^ - T2602 , 


or 

T,60\ - T„60„ 

6 T = ^ ~ . 

®2 “ ®1 

Assuming independent observations, the standard deviation of 
T is then given by 

(T 2 + T„2)'^ 


^1 ■ "2 
(T S3 

T 0. - 0, 


71 



Substituting for 6^ from 0^. = 


i V t" ’ 
« i 




V T T CT^ + T?S 

V 1 2^ 1 2 ^ 

B(Tj^ - T^) 


1/2 


^0 • 


(2-4-15) 


Srom Eg. (2-4-15) it is seen that tbe uncertainty in time -of 
flight is very sensitive to the uncertainty in pointing angle when 
the spacecraft is far from the target planet and decreases as the. 
spacecraft approaches the planet. ' The Equation emphasizes the 
importance of stars in the data since they minimize the contribution 
of pointing errors to Og. It is seen that -a larger B, which increases 
parallax, minimizes the error. A smaller also gives more 
parallax by decreasing the range at which the observations are taken. 
Finally, for a fixed measurement time, T^, Eg. (2-4-15) is minimized 
as T 2 is taken closer to encounter. 

Satellite Observations 

It is emphasized that the analysis presented here only applies 
far from encounter. As the spacecraft approaches the planet, 
parallax effects in the case of natural satellite observations allow 
one to solve for V^. Also -the time of flight solution becomes less 
sensitive to pointing errors. Furthermore, sufficient data will 
have been taken to estimate the natural satellite’s aphemeris relative 
to the target planet thus reducing effects of this error source. 

In the case of planet limb observations, cannot be accurately 
determined until planetary bending of the approach -trajectory occurs. 
In the case of Deimos, however, parallax effects are discernible 
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long before planetary bending occurs. In addition, the small size 

of Deimos makes image center-f inding errors negligible. These two 

- * 

factors make satellite observations significantly more accurate than 
Mars limb 'observations for approach navigation as we shall see in 
Chapters 4 and 5. 

Combined Radio and Optical Solutions 

Solutions which use a combination of doppler and optical data 
are of particular value during planetary approach since these data 
types complement one another. The primary error sources prior to 
encounter in solutions using only doppler data are target planet 
epheraeris errors, station location errors and non-gravitational 
accelerations. Optical data is insensitive to these errors since 
it directly relates the planet and spacecraft positions. On the 
other hand, as discussed above, optical, data suffers from the 
inability to determine accurately time of flight and velocity of 
the spacecraft, quantities which, for favorable approach geometry, 
are well determined by doppler data. The optical data and radio 
data complement each other in the navigation information that they 
provide. Hence, the combination of radio and optical data yields 
extremely accurate solutions and give a good estimate of encounter 
conditions much earlier than either data type taken separately. 

. For the radio plus optical strategies, two solutions could be 
generated : 

1) a combined radio plus optical solution obtained by processing 
radio data' and optical data simultaneously 
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2) a radio plus optical solution obtained by processing the 


radio and optical data sequentially; i.e. , after the radio 
data processing, using the a-posteriori covariance from the 
radio analysis as the apriori covariance for the optical 
data arc. This process should yield solutions that are 
very insensitive to radio data errors. 

The radio data in each case could consist of 

(i) only doppler 

(ii) only range 

(iii) both doppler and range. 
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CHAPTER 3 


OPTICAL NAVIGATION SYSTEM 

In this chapter we describe the optical navigation system required 
and the details of the methods to perform the navigation task using 
optical data. We will begin with a description of the overall functions 
involved for navigation in general. The sub functions are general but 
are described in the context of optical navigation and the details pre-- 
sented for the subfunctions pertain specifically to the use of optical 
data; details specific to the use of radio data can be found in 
Ref. 27. All subsequent sections of the chapter deal primarily with 
optical data planning and processing, with the exception of Section 3.6, 
which deals with the orbit determination process, including the overall 
d 3 mamic and measurement model description, for any data type. 

Figure 3^1 gives a diagram describing the functions involved in 
the navigation of a spacecraft. These fall into the five major cate- 
gories shown in the figure; 

(i) system modelling, 

(ii) data constitution planning, 

(iii) , measurement extraction and processing^ 

(iv) estimation, and 

(v) maneuver computation. 

The navigation process begins \<f±th a characterization of the sys- 
tem for both the dynamic and the measurement components. The dynamic 
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modelling (Section 3.1) involves a description of the time evolution of 
the position and velocity of the relevant celestial objects and of the 
spacecraft* This coupled with appropriate initial conditions yields a 
flight trajectory. The measurement modelling includes a definition of 
the data type used together with modelling the behavior of the instru- 
ments from which these data types are acquired (Sections 3.2 and 3.5). 

The , next step in the navigation process consists of data sequence 
planning and its evaluation from the standpoint of navigation require- 
ments. For optical measurements the picture sequence design (Sec- 
tion 3.3) consists of (i) determining the opportunities when the desired 
images can be viewed within the spacecraft and environmental constraints 
imposed; (ii) designing individual pictures with appropriate image 
geometry; and (iii) from a sequence of these, determining if adequate 
overall data coverage is obtained or if a new spacecraft trajectory 
design is required. The navigational evaluation of this picture 
sequence is based on requirements and constraints imposed by the over- 
all mission objectives. These Include the orbit determination, trajec- 
tory control and instrument pointing accuracies needed for the 
scientific investigations to be conducted. The generation of ‘these 
predicted accuracies is preceded by the linearization of the dynamic 
and measurement models about nominal conditions, which enables the use 
of linear estimation techniques. 

After an acceptable data sequence design satisfying navigation 
requirements has been obtained, the spacecraft is commanded to acquire 
the data. After execution of these commands and return of the measure- 
ment data to earthy this data must be processed (Section 3.4). for 
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optical data this consists of two streams of data — the data related to 
instrument pointing and fhe video data from the television cameras. 

The pointing measurement reduction consists of processing the data 
from the spacecraft attitude control system sensors and from the scan 
platform gimbal angle sensors. The TV data processing consists of 
defining the image profile and the subsequent determination of the 
image center. 

The use of the measurements for the estimation of parameters con-- 
sists of computing residuals (departures from expected measurement 
values) and then using these for the generation of the solutions, along 
with the statistical accuracies for these. In the calibration mode 
(Section 3.5) only the parameters related to the instrument model are 
treated. In the orbit determination mode (Section 3.6) all model 
parameters may be treated, with the emphasis being on spacecraft 
trajectory estimation and possibly on dynamical model improvement. 

Both functions require a linearization of the respective models (Sec- 
tions 3.5 and 3.7). 

Based on the spacecraft orbit estimate from the orbit determina- 
tion function, a trajectory correction maneuver can be computed. When 
executed, the maneuver alters the spacecraft state, thus yielding a new 
spacecraft trajectory. 

In the .sections to follow we discuss these functions in more 
detail, as indicated above. The chapter concludes in Section 3.8 with 
a brief description of the overall optical navigation system structure. 
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3 . L Spacecraft-Object Vector 


The geometry of obtaining the vector from the spacecraft to the 
object being imaged in the television camera was illustrated in 
Figure 2-2. This section describes the process of obtaining each 
of the component vectors required to establish the spacecraft-object 
vector. These are (a) the spacecraft vector, (b) the target planet 
vector (c) planet centered satellite vector (when the image is a 
natural satellite) and (d) vector direction to a star (for star 
observations) . 

Spacecraft Trajectory 

From Newton's law of universal gravitation the force of attraction 
acting on a particle of mass m^ due to a particle Pj of mass 
m, is given by 




G m.m. 

- 1 J 



(3-1-1) 


where G is the Gravitational Constant and are the position 

vectors of P^, respectively. From Equation (3-1-1), it follows 
that the equation for the motion of a spacecraft relative to a central 
body can be written as 


d^r 




(3-1-2) 


with 


£ = position vector of spacecraft relative to the 
central body 

“S/C (3-1-3) 
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and 


[i = G (M ^ ^c/p) ^ 


o . S/& o 


(3-1-4) 


r^y = position vector and mass of the central body 

^S/G position vector and mass of the spacecraft 
* th 

= k perturbing acceleration 

Perturbing accelerations are defined to be those that lead to 
departure from two-body motion, i.e., accelerations other than the 
second term in Eq. (3-1-2). The principal perturbing accelerations 
on the motion of the spacecraft arise from: 

1) N-body accelerations due to the gravitational effect of 
bodies in. the solar system other, than the- central body of integration; 

2) the departure of the gravitational field of the central 
body from spherical symmetry; 

3) . solar radiation pressure on the spacecraft; and 

4) propulsive maneuver thrusts by the spacecraft. 

The abovementioned are described below: 

1) W-body acceleration: The perturbing acceleration due to 

the ith gravitating body, is composed of two terms — the . 

direct acceleration on the spacecraft and the indirect acceleration 
due to the perturbing acceleration -on the central body. Thus, 
using Eq. (3-1-1), 


^,3 


hi - Ss/cl 


and 


r, - r , 


1=1 


(3-1-5) 


(3-1-6) 
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where 


'= position vector and mass of the ith body and 
N = number of perturbing bodies* 

2) Nonspherical Gravitational Field: The gravitational 

potential, V, at the point jr due to a distributed mass of density 
D is given by, 

dv(£) , (3-1-7) 

where D(^) is the density and dv(§) is the volume element at location 
5 , The acceleration due to V on the spacecraft at _r is then given 
by 

r - VV , 

where V is the vector gradient operator. 

The potential function can.be expressed as (Reference 27) 

J (^f p (si„0) 

n 

(sin0) fc cos m X -h S .sin mX 
n ^ L nm nm J 

(3-1-8) : 

where 

r, 0 , \ = body-centered radius, latitude, and longitude of 
spacecraft ^ 

a^ *= mean equatorial radius of body 
B^ (sin 0) - Legendre polynomial of degree n In sin 0 
B^ (sin 0) = associated Legendre function of first kind 




VCI) = 8 JJ-J 


D(§) 



G , S = numerical coefficients (tesseral and 
n,m n,m ' 

sectorial harmonic coefficients) • 

The associated Legendre function is defined by 

- m/2 ,m 

P” (V) = (1 - V ) ^ P^(V) , (3-1-9) 

dv 

where is the Legendre polynomial of degree n given by Rodrigues' 
formula 


•j 2 

^ I J ^ 

2 n! dy 


(3^1-10) 


The three terms inEq. (3-1-8) correspond to the potential of 

a point mass, zonal harmonics and the tesseral and sectional 

harmonics C and S • The perturbative acceleration can therefore 
n,m n,m ^ 

be derived from terms beyond the first in Eq. (3-1-8). 

3) Solar Radiation Pressure: There is a perturbing acceleration 

on the spacecraft, A£gp^ impact of photons from solar 

radiation. This was discussed in Section 1.3 where Eq. (1-3-1) 

gives the acceleration in the sun- spacecraft direction. In addition 
to this there are two smaller forces orthogonal to this because of 
the asymmetry of the spacecraft configuration. 

4) Propulsive Maneuver Thrusts: Accelerations due to thrusts 

generated by the spacecraft during propulsive maneuvers, 

given by 


^-PM m(t) ~PM'‘- ^ 


(3-1-11) 


where t is the time from the start of the burn, r_,(t) is the thrust 

^ — PM 

direction, m(t) is the spacecraft mass, and F(t) is the thrust 


82 ’ 



magnitude expressed as 



The mass m(t) at time t is 

m(t) - m<'t)dt , 

where ni(t), the mass flow rate, is given by 



A maneuver of very short duration can be regarded as an impulsive 
bum, as opposed to a finite burn* This can be represented by an 
instantaneous change in the velocity of the spacecraft, “^he 

corresponding change in the instantaneous position is 1/2 tj^, where 

t^ is the bum timeo 
Target Planet Vector 

The vector to the target planet is obtained from the precomputed 
position and velocity ephemerides for the celestial bodies in the 
solar system* These planetary ephemerides are obtained by a 
simultaneous numerical integration of the N-body equations of 
motion for the celestial bodies (Reference 28)* 

Planet Centered Satellite Vector 

The motion of the natural satellite in its orbit around a planet 
is obtained based on the analytical ephemeris theory developed by 
H. Struve and described in Reference 29* Wilkin^ s orbital elements 
(Reference 30) are used to define the coordinate system -- see Figure 
3-2, In this theory the orbital plane of the satellite is approximated 
to be inclined at a constant angle to a fixed plane, called the 
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Laplacian plane, upon 'which the ascending node of the satellite orbital 
plane regresses. Short period variations in the orbits are ignored^ 

The angles shown in Figure 3-2 are defined below: 

= longitude of node of fixed Laplacian plane on standard 
equator (1950.0 Earth Equator). 

= inclination of fixed Laplacian plane to standard equator 
“ the argument of the ascending node of the orbital plane on 
the fixed Laplacian plane 

- inclination of the satellite orbital plane to the fixed 
Laplacian plane 

L = the mean longitude of the satellite measured along the 

standard equator, the Laplacian plane and the satellite plane 
P = the longitude of pericenter of the orbit of the satellite, 
measured along the standard equator, the Laplacian plane 
and the satellite orbit plane. 

As mentioned earlier, is held a constant in the theory. The 
angles J^, K^, and P are modelled as linear functions of time 
given by 


= Ng + N^f 


•"a = 




(3-1-12) 


and 




•where the elements (O 7 the values of the angles (Oa at a 


specified epoch and the elements (On ®te their ratesj the time, t, 

K 
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is measured in days past the epoch. 

The mean longitude,’ L, is modelled as having a secular acceleration 
and a long period variation term. This latter term is due to Bom 
(Reference 31); it is the combined effect of the direct solar perturb- 
ation and an interati 9 n perturbation induced by the variation in 
inclination of the satellite orbit relative to the planet equator due 
to solar perturbation. The mean longitude is given by 


- sinCK^-K^)] 


(3-1-13) 


where is the coefficient of the periodic variation term and is 

the node of the planet orbit about the sun measured relative to the 

planet equator. (For Deimos L = 0.27 deg; for Phobos this effect is 

K 

negligible . ) 

At a given time, t, the position vector ^ of the natural satellite 
relative to the planet in the 1950.0 Earth Equator and Equinox Coordinate 
System is given by 


£ = NJKw Wr . 

where the rotation matrices are defined, as 


(3-1-14) 


N = 


cos -sin 0 


sin N. 


cos N. 0 
A 


J = 


0 cos J , -sin J . 

A A 


0 sin cos 
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K *= 


cos “Sin 0 


sin cos 0 




W ^ 


0 


0 cos -sin 


0 sin cos 


cos -sin 0 


sin W* cos W. 0 
A A 


and 


ij 


The rotation angle is defined by 

W = w + V 
A 


where 


w » P - K. - N* (argument of periapsis) 
A A 


(3-1-15) 


(3-1-16) 


and 

2 

V =® M + 2e sin M + 5e sin 2 m/ 4 + . . . (true anomaly) (3-1-17) 

and where 

M. = L - P (mean anomaly) (3-1-18) 

The radius vector ^ is given by 

r I 

r t* < o I (3-1-19) 

o I 

\ / 

where 


r » a 


/ l ..r. e! ■ A 

\ 1 + e cos V j 


(3-1-20) 
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and 

a, e semi-major axis and eccentricity of orbit of the 
satellite around the planet* 

Vector -Direction to a Star 

The apparent locations of stars are obtained from a dictionary 

of stars derived from the Smithsonian Astrophysical Observatory 

star catalog* Other catalogs of stars can also be used* These 

catalogs list the stars in terms of the right- ascension and declination 

in the 1950.0 Earth Equatorial and Equinox coordinate system; 

these are locations corrected for stellar aberration effects at 

the time of observation and are referenced to a heliocentric frame 

of reference. The direction to a star q, in that frame is given by 

I cos 6 cos a I 

( 

cos 8 sin or > (3-1-21) 

sin 6 I 

where a and 6 are the right ascension and declination of the star, 
obtained from the star catalog and corrected’ for proper motion since 
the epoch® 
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3.2 Observation Equation Model 


In this section we develop the observation equation in terms 
of the vector direction to the object, There are two steps 
involved in arriving at the observation in TV pixel- line coordinates 
as defined in Section 2.1: 

1) obtaining the direction vector to the object in the vidicon 
frame of reference; and 

2 ) mapping the vidicon referenced direction vector to pixel and ' 
line coordinates. 

We now discuss both of. these in the following two subsections. 

/\ 

1. Obtaining ^ in Vidicon Coordinates 

In Section 1.5 we developed the coordinate transformations to 
go from the frame of reference through the frames C, X, S^ and P, 
to the frame of reference T; these frames of reference are also 
defined in Section 1.5. 

The television camera system, T, can be directly referenced 
to the Celestial system, C, by its clock angle cone angle 
and rotation angle as shown in Figure 3-3. The transformation 
which takes the ABC system to the television MNL system (as 
opposed to the Platform MNL system), is' then given by 


T> — p p p 

^CT 0X30X201^ 


where 


R. 


CT 


1 


cos 

-sin 

0 


sin or^ O" 

cos or^ 0 

0 1 


(3-2-1) 


89 





Figure 3-3. Celestial/ Television Coordinate System Relationship 
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R, 


CT, 


and 


R. 


■CT. 


This yields 



0 - 

•sin 


1 

c 0 


0 

cos 

COS 

sin 

1 

o 

-sin 

cos 

Y« 0 

0 

0 

1 



cG.cB.cA - 

sG.sA 

cG*cB,sA + 

sG.cA 

-cG. s-B 

^CT 

-s6.cB.cA - 

cG. sA 

-sG.cB. sA 1 

cG.cA 

SG m sB 


sB.cA 


sB, sA 


cB 


(3-2-2) 


where the prefixes "s" and "c" refer to sin(*) and cos(*) respectively 
and where A, B and G refer to and respectively. 


Since can be written 
^GT ^T^CP 


(3-2-3) 


using Eqs, (3-2-2)^ (1-5-14) and (1-5- 15a) through (1-5-18) we 
can determine the television camera pointing angles given: (a) 

the platform pointing angles obtained from the clock and cone angle 
gimbal potentiometers, and (b) the pitch, yaw and roll angles obtained 
from the sun sensor and star tracker signals. It is useful to have 
these angles directly when we wish to solve for the camera pointing, 
Vidicon Coordinates 

As discussed in Section 2 pi, the definition of V, the vidicon 
coordinate system xyz, is determined by the vidicon electronics — 
specifically, the direction of scan during the image charge readout 
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process. This is rotated relative to the T coordinates through an 
angle p about the positive camera ^ axis, as shown in Figure 3-4. The 
corresponding rotation matrix is given by 




cos p sin p 0 

'Sin p cos p 0 

0 0 1 


(3-2-4) 


The vector ^ expressed in inertial coordinates can therefore 


be transformed to vidicon coordinates, _t^ through 


where is given by 

^IV \v- ^CT 


(3-2-5) 


(3-2-6) 


It may be observed that the rotation R could have been absorbed 


along with in. the rotation . However, it is convenient to 


GT 


define the such that it has a small value ('^ 0). Then, if the T 
offset angles \Ji, x>.^ small, the angles of^, and y^ can be 
easily determined from similar angles for the platform, a, P ,. y through 
the approximate equations 
or + x/sin P 

« P + ^ ' (3-2-7) 


and 


Yj Y ■ X cot P + U) 


instead of through equation (3-2-2) , We can make these approximations 
because the rule of vector addition holds for infinitesimal angular 
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rotations . 


2) ImaRe Location in Pixel and Line Coordinates 
Mapping of the vector coordinates through the camera optics 
onto the camera image plane is quite straightforward using the 
collinearity equations of photogrammetry (Ref. 32), and is given 
by 



(3-2-8) 


where f is the focal length of the optical system and 

T 


^ = t t, t 

-V V 3 


(3-2-9) 


1 ''2 ''3_ 

The transformation from the image plane to pixel and line 
coordinates is obtained by a mapping through a scale factor matrix 
K and a translation of the origin, i.e., 

(P) 


K 


a -Cl 


(3-2-10) 


Here (p , X ) are the coordinates of the central pixel and line of the 
o o 

TV target raster; K is given by 


K = 


k k 
X xy 


k k 
yx y 


(3-2-11) 


where the diagonal elements of k are the scale factors from the 
image plane to the pixel- line plane while the off-diagonal terms 
provide a rotation with respect to x-y coordinates and a non- 
orthogonality of the P, L axes. 
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3 .3 Picture Design 


In this section various quantities necessary when designing a 
picture are computed . 

In Section 3.2 we went through the development of the equations 
for the locations of the observed images in the picture. In addition 
to this, during the picture design, the geometry in both inertial 
coordinates xyz and in celestial coordinates ABC is important, due 
respectively to the requirement of having a good star background 
and the requirement that scan platform movement constraints not be 
violated. 

If a particular right ascension, a, and declination, 6, direction 
(in inertial coordinates) is desired, will be given by 
cos a cos 6 
sin a cos 6 
sin 6 

The clock and cone angles, a and p, of the direction can then be 
obtained using 



(3-3-1) 


-c 



(3^3-2) 


atid 

a B tan“^ (t ft ) 

C2 

P = cos"^ (t ) (3-3-3) 

^3 


where 0 s: a < 2tt and 0 ^ p ^ rr . 

Similarly, for specific clock and cone' angles, using 
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I cos O' sin p 


^ = (sin a sin p 

-c 


(3-3-4) 

\ 

cos P 

and using Eq. (3-3-2) we can obtain the right ascension and declination 


a = tan ^ (t ^ 

H n 

8 » sin ^ (t ) 


(3-3-5) 


where 0 ^ a < 2tt and -rr/2 < 6 ^ rr/2 . 

For ease and accuracy of image detection there should be a 
minimum of smear caused by the motion of the image during the time 
the camera shutter is open. To compute the velocity of an image in 
pixel- line coordinates, let _t be the velocity of the object. Then 
t. 


r 

^2( 


-V “ ) ^ 2 ) " ^3 


(3-3-6) 


From Eqs. (3-2-8) and (3-2-10) we can write 


= f[K] 


where t. ^ I tl t„ 
X ' V, 


' VSi 


Differentiating we obtain 


- f[K] • ^ 





*^ 3^2 “ ^ 3^2 


(3-3-7) 


(3-3-8) 


where v^, v^ are the image velocities in the pixel and line directions 
respectively. In this computation we have not considered the image 
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velocity due to the change in in Eq. (3-3-6) caused by the 
attitude control limit cycle motion. 

To determine the illuminated shape of the image, the phase 
angle of the object is required, where the phase angle is defined 
as the angle between the directions from the object to sun and to 
the observer (spacecraft) respectively. Thus the phase angle gives 
a measure of the departure from full illumination of the object 
as viewed by the observer, zero degrees implying full illumination 
and 180° implying zero illumination. The phase angle, \jr, is 
computed by 

lit = cos"^(2 ‘ ^ (3r3-9) 

where 

^ = sun-object unit vector. 

The brightness of an image is characterized by the visual 
magnitude V of the object as viewed from the spacecraft. For stars 
this is available directly fromthe star catalogue; for other 
celestial bodies it is computed by 

V = + 5 loSj^Q ty - 2 • 5 logj^Q 3> (3-3-10) 

where 

y, t = distance of the object, expressed in astronomical 

units (AU) , from the sun and spacecraft respectively, 
« visual magnitude of the object as viewed from the 
sun at a distance of 1 All, ' 
and $ is computed from 

^ = — [sin ilf + (tt - ilr) cos ilij , (3-3-11) 

TT ^ 
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An important consideration during the design of pictures 
containing the natural satellite is the avoidance of the planet 
in the field of view. Exposure times that are suitable for the 
satellite images could cause damage to the camera if the much brighter 
planet were imaged. Even if the planet is not in the field of view 
but close to the edge, problems could arise in detection of the 
images within the field of view due to stray light from the planet. 

For typical interplanetary spacecraft, instrument pointing control 
accuracies are of the order of 0.5 degrees (3 cr) in any direction; 
this then causes a (3 a) constraint to be placed on the picture 
design that the separation angle between the edge of the satellite 
image and the edge of the planet be larger than 1.0 deg — i.e., 

0.5 deg to avoid imaging Mars and 0.5 deg to ensure imaging the 
satellite. 

To compute the separation angle, first the semi-angular diameter 
0 of both the satellite and the planet are computed by 

0 = sin R /i tl , 

0p = sin’^ Rp/ls| (3-3-12) 

where 0 = ^emi- angular diameter, 

R » radius of the object, 

and subscripts n, p refer to the satellite, planet respectively. 

The separation angle is then given by 

0 -n-(0 + 0), (3-3-13) 

sep ' n p' ’ \ / 

where fj is the angular separation between the directions to the 
centers of the satellite and planet and is given by 
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( 3 - 3 - 14 ) 


^ / A A. “1 y-A A. 

U ^ cos ^ ^ ^ ^ (s ♦ t) 

If the value of f^rom Eq, (3-3-13) Is negative, this implies 

that the satellite is either (i) occulted by the planet or (ii) is 

in transit across the planet as viewed from the spacecraft. These 

two cases occur respectively when (i) |^| < |^j or (ii) 


t . 
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3.4 Measurement Processing 

As discussed in the beginning of the chapter, there are two 
categories of measurements required for optical data processing, viz* 

(i) measurements that lead to a determination of the television camera 
pointing direction using the spacecraft telemetry engineering data; and 

(ii) measurements using the television picture data leading to a deter- 
mination of the image location in picture coordinates. These are 
referred to in the following as pointing data and imaging data respec- 
tively, and are discussed below. 

Pointing Data 

The process of obtaining the pointing of the television camera 
from spacecraft sensors was described in Section 1.5. This consists of 
(i) obtaining the inertially referenced spacecraft attitude from the 
sun-sensor and star- tracker data (i.e. the transformation followed 

by (ii) the use of the scan platform clock and cone angle data to obtain 
the spacecraft relative- platform orientation, Rgp 5 ^nd therefrom the 
television camera orientation, Rg^, using 

^ST ^ ^T ^SP (3-4“l) 

These two pieces can then be combined to give 

=*IT - ®ST hs • 

which determines the pointing direction, of the ^television camera in 
inertial coordinates , or equivalently in celestial coordinates (^) 

using 
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The television orientation can be expressed in terms of the 
television clock, cone and rotation angles, and Ytji> defined in 

Section 3* 2* 

Imaging Data 

The raw imaging data consists of video intensity for each pixel in 
a picture, discretized to 2^^ gray levels; n typically ranges from 
seven to nine. As mentioned at the beginning of the chapter, the 
reduction of this raw data to image location information can be divided 
into two functional steps, viz. (i) some method of defining and deter- 
mining the profile of the lit limb, and (ii) a method for estimating 
the center of the viewed object as mapped into picture frame 
coordinates* 

For small images (-1-5 pixels), such as those of stars and of 
small diameter natural satellites , both of these functions may be per- 
formed by display and visual inspection of areas around the expected 
image locations in the form of plots showing video intensity as a 
function of line and pixel number. The detection of dim images can be 
facilitated by contrast enhancement within the picture. 

For large images such as that of ^ planet the processing has to be 

more elaborate to minimize the error in locating the center. The lit 
\ 

IxTDb is defined by a curve joining all points with the same predeter- 
mined level of video intensity. These points are located by reading 
the digital video data along a direction parallel to the P or L axis of 
the vidicon (P, L) coordinate system. A process is then required to 
determine the center of the image from these limb points. 
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If the planet is modelled as a three-dimensional ellipsoid the 
cross-section observed from any direction is given by an ellipse; the 
parameters of the ellipse are a function of the orientation of the 
ellipsoid relative to the line-of-sight ("look") direction, from the 
spacecraft at the time of observation. In the following we give a 
development to determine these parameters. 

The equation of an ellipsoid oriented with its principal axes along 
the coordinate axis is given by 


2 2 2 

X y z 

2 ^2 2 

a b c 

o o o 


(3-4-4) 


where a 


o' o' o 
respectively; or 


c are the semi-major, semi-mean and semi-minor axes 


$ = X Ax = 1 


(3-4-5) 


where 


X A 


X 


L^o. 


(3-4-6) 


and 


A == 


l/a' 


1/b: 


lie 


(3-4-7) 


Let the vector x be represented by x„ in Ellipsoidal body-fixed 

' Tj 


coordinates, E, where the z axis is aligned with the planet pole and the 
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X axis is at longitude X. relative to the Prime Meridian reference, M, 
for the planet. Thus 

(3-4-8) 


— E ^ME 


where 




cos X 
—sin X 
. 0 


sin X 
cos X 
0 


0 

0 

1 


(3-4-9) 


Let the planet be rotating at angular velocity w about its pole 
and the planet pole direction be given by right ascension a and 
declination 6 in inertial coordinates, I; then 

~ —I (3-4-10) 

where, letting c(*) and s(«) represent cos(*) and sin(*)» 


IM 


c(u)t) s(u>t) 0' 

-s(cot) c(oot) 0 

0 0 1 


s6 

0 

c6 


0 

1 

0 


-c6' 

0 

s6 


' ca 

SQ 

0" 

-sa 

cff 

0 

. 0 

0 

1_ 


(3-4- 

-11) 


The rotation to vldicon coordinates, V, is given by Rj-^. obtained 
from Section 3.2; then we can write 

— E ~ ^ (3-4-12) 

where x,j, is the vector from the center to the surface of the ellipsoid. 


expressed in vidicon coordinates, and where R__ is given by 

Vri 

^ 

Substituting this into Eq. (3-4-5), . 


T 

Sv 


§ = x^, B ^ = 1 


(3-4-13) 


(3-4-14) 
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,where 


B - Rev* Rve (3-4-15) 

Now, the cross-section of the ellipsoid seen will be an ellipse 
defined by those limb points where the tangent is parallel to the look 
direction _2, i.e. the gradient 9$/93C will be normal to the vector 
Therefore 


9f 


= 0 


9 3^ -V 

where ^ is _£ expressed in V coordinates, given by 

'o' 

iv ■ 


(3-4-16) 


(3-4-17) 


Differentiating (3-4-14) and using B = B , Eq. (3-4-16) yields 

'o' 

^ = 0 (3-4-18) 


T 

2^2 


'lO 

1 


Using this in Eq, (3-4-14) will yield the ellipse, the projection 
of which in the plane is the desired limb profile. We would 

like to define a coordinate system U for which the reference plane is 
parallel to this ellipse; defining u such that 


or 


u H {u- 


= C u 






(3-4-19) 


(3-4-203 
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where 


c ^ b“^ 


-3_ 

and where A is given by 

0 
o 


A-l = 


(3-4-21) 

(3-4-22) 


0 


0 


(3-4-23) 


In obtaining Eq. (3-4—22) we used the fact that the R*s are orthonormal 
matrices. Using Eq. (3-4-19) in (3-4-18), 

= 0 


T 

u •< 


or 


0 

IJ 


“3 = 0 


(3-4-24) 


Substituting Eq. (3-4-20) into (3-4-14) and using (3-4-21) , 
T T 

$ = u C BC u = 1 


T -1 

or, sxnce C = C = B , 


1 = $ = u C u 


■ “2 1 “3] 


i . <^2 

1 

pT I p 

L 2 ! 03 J 


= [Uj_ U2] 




Uo 


U., 




U. 


(3-4-25) 


(3-4-26) 


(3-4-27) 
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From Eq. (3-4-26) and (3-4-20), 



- 

r "1 
"i 

J ' 


to 

to 

1 

1 



(3-4-28) 


or 











L 2j 


L 2j 


Substituting (3-4-29) into (3-4-27) , 


(3-4-29) 


1 -1 

L 1 2J 






= 1 


(3-4-30) 


Thus, if is obtained as the upper left hand 2x2 partition of 
C which is obtained using Eq. (3-4-22), then Eq. (3-4-30) gives the 
projected limb profile in vidicon coordinates, translated such that the 
center of the ellipse is at the origin. 

Let this ellipse be represented as in Figure 3-5 with the x-y 
axes parallel to the principal axes, where a, b are the semi-major, 
semi-minor axes of the limb profile; Pj^, are the coordinates, of the 
center of the ellipse; and i|j is the orientation angle of the ellipse as 
shown. The equation of this ellipse is 


2 2 


(3-4-31) 


or 


cl 2,2 2.2 _ 

(1 r-j X + y - a + c = 0 

a ' 


(3-4-32) 
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t 

L (DIRECTION OF INCREASING 
SCAN LINE NUMBER) 


Figure 3-5. Planet Limb -fitting Process 
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where c is the distance of the focus from the center given by 

A1/2 


c * (a^ - b2) 

The observable equation is defined to be (Ref. 190 
2'^ 


<Jr = 1 - 


2.2 J- ^ 2 , i, 

X +y - a + c +§ 


(3-4-33)“ 


(3-4-34) 


where t is the observation noise and where (x, y) are given in terms of 
the limb points (p, f) by 

cos -sin il;" 


P - Pi- 


f - 


(3-4-35) 


[_-’sin 4* ■ . "COS tpj 

A minimum variance estimation algorithm- is used to identify the param- 
eters a^ c, i}^, p^ and 

To determine the effect of errors in the limb points, we take 
the variation of Eqs. (3-4-34) and (3-4-35), 

(3-4-36) ■ 


7 c^\ 1 


(6x1 

64>= 2 1 - ^ X y 



U, aj 

- 



= -2 


I 

-|1 - X cos + y sin 4^ J |1 j I ^ 4^+7 cos 4> 

a / I \. a- 


6p 

bH 


(3-4-37) 


We will use this expression in Chapter 4. 


06«P0^- 
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3^5 lustrumeutaliloii Error Models 


'll! this section we presenl: the mathematical models for the 
Instrument related errors in oi)tical navigation measurement processing. 
Referring to the classification of’ optical data errors in Section 2.3 
the TV pointing errors and TV distortion errors stem primarily from 
spacecraft related errors and camera related errors respectively. 

These are discussed below. 

1. TV Pointing Errors 

The apriori TV pointing is determined using the attitude 
control system pitch, yaw, and roll sensor signals and the cioch 
and cone gimbai angles for the scan platform. In addition to the 
error in the signals themselves, there are contributions to it from 
mechanical misalignments in the spacecraft instrumentation. We 
discuss models for these here. 

Effect of Mechanical Misalignment Errors 

When there are instrument misalignment errors to be modelled, 
the observable equation has to have these errors factored in. In 
general' the misalignment error at any step can have components in 
oach of the three coordinate axes; consider an error rotation about 
the % coordinate ^is in an xyz system. This will be given by 


which for small values of €3 reduces to 


cos 


sin 

®3 

0 

-sin 

to 

cos 

®3 

0 - 

0 


0 


1 
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1 G3 0 

■63 1 0 


(3-5-1) 


= (I + E 3 ) , 


where I is the identity matrix and 

r 




0 63 0 

•63 0 0 


000 

Using similar error rotations, and about the other two axes. 


and writing 

R = I + E. ■ 

6 . 1 
1 

the general error rotation transformation can be written 
= (I + E 3 ) (X + E 3 ) (I + 

or, neglecting higher order terms, 

R^ = I + (E^ + E 2 + E 3 ) 


(3-5-2) 


i-5-3) 


(3-5-4) 


= I + E 


where 


E » 


0 


^3 

0 


-e. 


-e. 


0 


(3-5-5) 


We note here that if R^ is a rotation through an angle 9 ^ (not 
necessarily a small angle) about the ith axis, then E. commutes 


with R. , i.e. , 
1 


R. E. — E.R., 
11 11 


i = 1 , 2, 3 


(3-5-6) 
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Using the noietion in Section i.5, the observation vector in 
vidicon coordinates without the effect of errors is written as 

~V “ ®TV ^ P- ^SP \s ®X X ®IX -I 

3 2 1 or O' 0 0 

Each of the R's after is a one axis rotation except for ^ 

which is a three axes rotation where the angles are small. 

TXo 

is an idealized rotation of coordinates. Thus with the effect 
of alignment errors the coordinates ^ will be, in general, 

4, ■= V (I + Rj) (1 f Ep^) (I + Ep^) Epp^ (I + Ep^) 

*p (I + E„) Ejp Cl + E ) Rjj (I + E^) Ejj p t 
O' a 0 0 

(3-5-8) 

where each of the E's is a three axis error rotation matrix as 
represented in general by Eq. (3-5-5). 

If new we make the assumption that .each of t]), y., and uo are 
small (when this assumption is not true all the terms must be 
carried through) then each of the Rp„ is a small angle rotation 
which can be represented by Eq. (3-5-2); thus the terms from (I ^ E^) 
through (I + Ep ) in Eq. (3-5-8) can be represented as (I + E^) 
and we can write as in Ref. , 

- Rpv Bpp (I + V P + V V. « + V 


\s ^X^ ®XqX ^Xq~I 


(3-5-9) 


Neglecting higher order terms this givws 




(3-5-10) 
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where 



and where J represents' P, o;, S, K and Q^. are respectively the 
products of the U’s preceding, following the corresponding + Ej.) 
term in Sq. (3-5-9) . 

At this point we note that, any antisymmetric matrix U 
operating -on a vector V can be represented, as the cross product of 
two vectors. In particular, if 

U = 

T 

then have 



Ujv = {3-5-11) 

and the error terms in Eq. (3-5-10) can be written 

* A 4’^ 

“ -cPj%) X <'a4^ 

■ A%> 

or, 

wA " 
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(3-5-13) 


where is the rotation vector angle corresponding to E^, ■ 

^ = r t t t 

and , ' . 


T = 


r 0 

V 


t 

■ 

'\i 

03 

0 

01 

02 

~^V 

01 

0 


• (3-5-14) 


Using this result in Eq. (3-5-9) we obtain 

-V ~ P%r ®P P^SP -S ®P P^'SP (3-5-15) 

0 O' (X or or O' 

where ^ is given in Eq. (1-5-18). Using Eq. (1-5-14) and 

replacing cos and sin by c and s. 


^SP " ^ f^SP 

O' O' 


and using Eq. (1-5-17) and combining terms. 


o»pC?p 

“fpCpp 

-Sp. 


oy_ 

0 

-sc^p 

p 

CQfp^p 

S^pSPp. 



(3-5-16) 


’’xp \p®xs 




: sCo'x ■ “p) 

• "p’ 


-C(cifjj. - Q?p 

sia^ - ofp) 


^P 

a 

-cp 


(3-5-17) 

The errors arise from the attitude control system pitch, yaw, 
and roil angle determinations. These angles, are determined, using 
the voltage signals from the sun sensor and star tracker. Modelling 
this determination simply as 
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(3-5-13) 


9. = Tc.v. 

X XX 


where are angles whicli are a function of the voltages, and 
scale factors, the error is given "by 


69i - 8k^l/^ 


6k. 

k.dy. d- — ^ 0. 

X X ^ 


or 


69. — e . + X .0 . 

X X XX 


(3-5-19) 

where are sensor angle null offsets (biases) and -are percentage 
errors in the scale factors. The e„ can therefore be written as 



i ( 

\ ^ 

1 \ 


> + < 

)y| 

i } 


1 

A r 

f 1 




y y/ 
X 9 

r r , 


(3-5-20) 


where the subscripts denote the pitch, yaw, and roll components 
respectively , 

There is a redundancy of error parameters in Eq. (3-5-15) due 
to Eq. (3-5-6). The four ^ are transmitted successively through 
three single angle rotations and E^g about axes numbers 

Of Of 

2, 3 and 3 respectively. There are therefore three redundant parameters 
the second and third components of ^ will have the same effect as 
the second of ^ and the third of ^g, while this latter will have 
the same effect as the bias portion of the third component of 
also. This is also evident upon examination of Eqs. (1-5-18), 

(3-5-15) and (3-5-17). 


114 



If 



then Eq. (3-5-15) can be written 


(3-5-21) 


(3-5-22) 


(3-5-23) 


(3-5-24) 
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The B e in Eq. (3-5-24) maps the effect of all the error rotations 
6 “ 

to platform coordinates. If we let 


b 

m 



= B e 

e “ 


(3-5-27) 


then similar to Eqs. (3-2-7) we can write 

Q, s= Of _ h /sin p . (3-5-28) 

r m 

Pj - p + (3-5-29) 

Yp = Y \ P + b^ (3-5-30) 

2* TV Geometric Distortion Model 

As discussed in Section 2.3 distortion effects in the television 
arise from (a) the camera optics during the imaging process and (b) 
the vidicon electronics during the photoelectric charge readout 
scan process. Corrections to correct for these effects need to be 
added to the vidicon (x, y) coordinates expected for an image (see 
Ref. 26). 


Optical Distortion : The optical distortion can be represented 

as the sum of four components ^ as follows: 

(i) a null offset of the optical principal point- from the 
target raster center given by 




aYo]- 


(3-5-31) 


(ii) a symmetric radial distortion component given by 


r n 


i2 = 


1 

2i-2 

2^ ^2i-2 ^o 

L i=l 


% 


(3-5-32) 
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where ^ is the image location relative to the optical 
principal pointy obtained by ijsing Eq. (3-2-8) and (3-5-24); 


^0 = i^l 


and 


^2±—2 ~ radial optical distortion coefficients j 

(lii) asymmetric radial and tangential distortion effects given 

— o 

by 


l3- 


n. 


2i 

2_ j '^Zi-l ’^o 


i=l 


-sin 0 


cos 0 


(3-5-33) 


where 


. and 


= asymmetric optical distortion coefficients 


0^ = orientation angle of axis of maximum tangential 


distortion; 

(iv) a distortion caused by lens misalignment, given by 


i4 = J 


jx y 
o o 


-yn 


X y 

O'' o 



'f/ 


1 

■ 

h- 

J 

1 ^ 

% 


(3-5-34) 


where are the corresponding coefficients. 

Relative to target raster coordinates the location of the image 
after optical distortion effects but before electromagnetic distortion 
effects is given by 

4 


X = X 




(3-5-35) 


1=1 
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Electromagnetic Distortion; The electromagnetic distortion can be 
represented as the stun of two components as follows: 

(i) a symmetric radial distortion component, given by 


% = 


n. 


E 2i 

’^2i-2 ^E 


2i-2 


i=l 




( 3 - 5 - 36 ) 


where 




'E 


r = x_ - 4 


4 = electromagnetic distortion null point 


^E -E 


and 


V 2£_2 ~ symmetric electromagnetic distortion 
coefficients; 

(ii) a s}rmmetric tangential distortion component, given by 


( 3 - 5 - 37 ) 


The location of the image after the distortion corrections is 
given by 

2 





“^E 

^2 = 

\ ^ 2i-l 

''2i-l ^E 

* 


i=l 


."E . 


X = X + 




i=l 

4 


= X + 


i=l i=l 


( 3 - 5 - 38 ) 
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The error models described above can be used for calibrating the 
TV camera optical and electronic systems through the use of inflight data 
(Ref. 33). 
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3^6 Orbit Determination Models and Processing 

In this section we describe the process of spacecraft orbit 
determination using optical data* The basic data processing is con- 
ducted using the well-developed theory of the discrete Kalman filter 
(Kef* 34) for a linear dynamical system* However , to use the Kalman 
filter equations, it is necessary for the non-linear equations to be 
linearized about a nominal spacecraft trajectory. The spacecraft state 
vector can be augmented by ether model parameters to be simultaneously 
es timated. 

We develop the models for the processing here; the actual optical 
data equation linearization is described in the following section. In 
addition to obtaining an estimate and statistics for the state vector, 
a sensitivity analysis can be conducted to determine the effect of 
unestimated error parameters and of incorrect modelling of tKe apriori 
statistics. To facilitate interpretation and evaluation of different 
estimates it may be desirable to map them along with their statistics 
to -a more convenient time and coordinate system. 

The orbit determination process then consists of the following 
functions : 

(i) Generating a nominal trajectory based on assumed initial 
conditions and nominal model parameter values, by integra- 
tion of the spacecraft dynamical equations; 

(ii) Integration of the spacecraft variational equations to 

obtain the variation in spacecraft state for variations in 
nominal initial conditions and parameter values. This is 
required for the linearization about the spacecraft 
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trajectory and is conveniently done simultaneously with 
step (i); 

(iii) Computation of the data partial derivatives, required for 
the linearization about the nominal measurement equation, 
using the result from step (ii) ; 

Civ) Computation of the data residual — 'deviation, 'of the 
measurement from the nominal predicted value of the 
observable; 

(v) Filtering the data residuaJs to obtain the optimal -estimate 
and statistics, under the assumptions, for the spacecraft 
initial conditions and model ''parameters; 

(vi) Mapping the resulting spacecraft state estimate and 

statistics to the desired time and coordinate system .for 

I 

proper evaluation (see Appendix A) . 

In the following we 

(i) describe the dynamic and measurement models along with the 
treatment of stochastic variables; 

(ii) present the Kalman filter equations with a consider option. 

Dynamic Model 

Let the six dimensional spacecraft state vector be represented by 


H(t) 4 4 V 
lf(t)J 


( 3 - 6 - 1 ) 


and the equations of motion for the 'spacecraft dynamics (discussed in 
Section 3.1) be put into the form 

u(t) = f[u(t), v(t), t] (3-6-2) 
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where consists of the set of constant parameters affecting the 

dynamics of Some of the parameters could be modelled as stochastic 

dynamic parameters — we discuss that later. 

Given some initial conditions on u and the vector v , where 

— —o’ 



the nonlinear equations of motion (3-6-2) can be integrated to yield a 
nominal trajectory for the spacecraft. This can be symbolically written 

‘] ’ 

where the subscript N denotes "nominal,” This trajectory is required 
to compute the nominal values of the observables which we discuss in 
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the following subsection. Also required is a linearization about this 
nominal and a computation of the effect of small departures from the 
nominal , 

Taking variations, the state deviation from the nominal trajectory, 

Eq. (3-6-2), is given to first order by 

v(t), t] 9f^[u<t), v(t), t] 

= 3u(t) 9v(t) 

(3-6-6) 

Where the variation 6(*) indicates departure of the function from 
nominal, i.e. 

60(t) - 0(t) - 


and where 


ai raf^ afi 
9v Lap^ BpJ 


The variation in v(t) can be written 


9v(y , ±) 

^ il, 

— o 


which, using Eqs. (3-6-5), can be expressed as 


6v(t) = 

'9P j 9£ ‘ 

9p [ 9a 

“‘-o — o 


► 


:7T'7] 


^6a ^ 


(3-6-7) 


(3-6-8) 


(3-6-9) 


Here 9p/9p^ and 9£./9^ are assumed to be available as precomputed 
quantities. Note that and hence v^, are not restricted to the same 
coordinates as p^, v, and could in general have different descriptions 
from £>, V respectively. 
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Since by the chain rule 

a^(u(t), w , t) 8f^(u(t), v(t), o t) 


3v 


3v 


3v(t) 

then using Eq;. .(3^6-8) in (3-6-6) we can wri-te- 

81(u(t), V , f) 3f(u(t), V , t) 


(3-6-10) 


6u(t) 

N* 

6v 
— o 



3u(t) 

0 


— c 

0 





6u(t) 


6v 


—0 


(3-6-11) 

where the lower partitions in the matrix are zero because the are 

constants • 

The solution of this can be written in the form 

6u(t)l 


6v 

“-0 


-u(t, t^) 

vct, t^)' 


6u(t ) 
— 0 

0 

I 


6v 


- 


^o 


(3-6-12) 


where, again, th.e zero and identity matrices in the lower partitions 

arise because the v ’s are constants. The state transition matrix must 
— o 


obey the equation (Ref. 35) 






d 

U 

V 


dt 

0 

I 



8f 3f . 


9_u 9^ 

0 ' 0 


U V 


0 


(3-6-13) 


That is, U and V must obey 

dU(t, t ) 9jf(u, y, t) 


dt 


3u 


U(t, t^) 


and 


dV(t, t^) 9f^(u, V, t) 


dt 


9u 


V(t, t^) + 


(3-6-14) 


t) 9v(v , t) 


3v 


,3v 

“O 


(3-6-15) 
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with tha initial conditions 


U(t , t ) = I 
o’ o 

(3-6-16) 

. '('o- ‘o> - 0 

Thus integration of the spacecraft variational equations (3-6-14) 
through (3-6-15) yields the solution Eq, (3-6-12), 

We note from Eq. (3-6-12) that the transition matrices U(t» 

V(t, t^) are the partial derivatives of the spacecraft state at time t 
with respect to spacecraft Initial conditions and other constant 
parameters, i.e* 


and 


U(t, = 


3u(t) 


(3-6-17) 


8u(t) 

. (3-6-18)- 

“O 

We shall use them in the following development for the lineariza- 
tion of the measurement equations. 

Measurement Model 

The general measurement equation could be written as a non-linear 
function of the vectors u(t), v(t) , w(t) given by 

_z(t) = li^^u(t^)> Vq 5 tj, v(a^, t), w(t), tj + n(t) 


(3-6-19) 


where _n(t) is the data noise. Here ii(t:^) are the spacecraft initial 
conditions and ^ are dynamic parameters affecting the spacecraft 
trajectory. We note that the parameters affect the data through the 
spacecraft state _u(t) and also directly through v. For instance these 


ommu pxge is 
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could be planetary ephemerls parameters which affect the gravitational 
influence on the spacecraft dynamics through the planet position vector; 
the latter however may also directly affect the data since the planet 

m 

image can be the observable. 

The w(t) are physical quantities that affect the data but are 
dynamically uncoupled from the spacecraft dynamics* In general w will 
consist of variables _q(t) and constant parameters ^ 


-a(t> 

w(t) = • 

• b 


in general obeying the dynamical equations 


i(t) 

b(t) 


[^(^(t), b, v(t), t' 
0 




} 


( 3 - 6 - 20 ) 


( 3 - 6 - 21 ) 


Let the solution of these equations be 


q(t) 



b 


b 

— o J 


( 3 - 6 - 22 ) 


where is a set of constant parameters. 

The ^(t) could be, for example, the natural satellite position 
vector — which for small satellite masses would not perturb the space- 
craft yet would be themselves perturbed by the planetary y_ parameters 
or image size proportional center finding errors, which are a function 
of time through their dependence on range to the target. The ^ could 
be instriament error parameters. In addition, some of the parameters 
could be stochastic — we discuss those later. 

Based on the nominal spacecraft trajectory ii(t)^ and nominal time 
functionals v(t)^, w(t)^ we can construct the nominal measurement 
vector 
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z(t)^ = h 




Given the actual measurement vector _z(t) , we can construct the data 
residual vector 

6^(t) ^ z(t) - C3“6“2^ 

given by the model 

6^(t) - 6h(t) + n(t) (3^6-25 

Using Eq* (3-6-22) in (3-6-19) we could write 

_z(‘t) = y V , w , t] + n(t) (3-6-26 


Defining 


= i 


m A m(t ) 
— o — ^ o' 


we can wrxte 


where 


6m^ = K ^ bin 
“t t , t ^ 
o 


U(t, t ) V(t, t ) 


9v(t) 


(3-6-30) 


9w(t) 9w(t) 
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the transition 



The vector m in the previous subsections was developed as consist- 
ing of uncertain but constant initial conditions and bias parameters* 

In addition to these it is often necessary to model some of the process 
noise effects on the spacecraft dynamics or the errors on the measure- 
ments, from a' variety of causes, as random time varying phenomena. This 
is done through a set of stochastic variables ^(t)* These could be time- 
varying representations of the bias parameters Included in any of the 
vectors u , v , w or could be due to unmodelled parameters or effects. 
Commonly used quantities, for instance, are random nongravitational 
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accelerations on the spacecraft, primarily useful for radio data, and ' 
camera pointing errors for optical data. 

These stochastic parameters can be treated by augmenting the 
state vector to include these quantities. Thus defining 


K(t) A i 


(3-6-34) 


L^(t) J 


the vector x(t)' becomes the quantity to be estimated. For the estima- 
tion of 2 c(t) we need to describe its evolution by determining the state 
transition matrix $(t, t^), such that 


ro-' 


x(t) = ^(t, x(t^) +• 


(3-6-35) 


A large variety of random processes can be modelled by the linear 
stochastic differential equation representing exponentially correlated 
process noise — i.e, the Langevin equation 


d(t) = - i s(t) + |(t) . ■ (3-6-36) 

where t is the correlation time constant and ^(t) is a Gaussian purely 
random (white) zero mean process, i.e. 


E(e(t)) = 0 


(3-6-37) 


with the correlation function 


E(e<t) |(t)) = Q(t) 6(t -t) . (3-6-38) 

This can be integrated and converted to discrete form (Ref. 36) 

®i+l = *^i®i ^i (3-6-39) 


where 


= exp (-At^/r) 


At. 


1 


*^i+l 


t. 

1 


(3-6-40) 

(3-6-41) 
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and 


{ - +i) 

From Eq, (3-6-39) and using Eq. ('3-6-42),- 

= ■I'i '^fi) - '•'i) 


(3-6-42) 


(3-6-43) 


which shows that in the steady s.tate 

This description can be used to represent purely random noise with 
T — 0, randomly varying drif-ts using an appropriately non-^zero finite 
value of T, and constant parameter for r oo ; this indicates the flexi- 
bility of this stochastic sequence. 

The -transition matrix for a vector ^ of these random ^processes is 
given by 


s.,T = [S.,- ,] s. + 
““X+l i+l,!-* —1 “X 


(3-6-45) 




s 


i+l,i 




(3-6-46) 


L 

where j represents the j th component of the vector _s . 

The transition matrix ® of Eq. (3-6-35) for the vector x^(t) is 

obtained by deterministically mapping the effect of into the 

"initial condition" deviation vector 6m (t.,,) using the stochastically 

— o x+1 ^ 

mapped deviation, of 6iS;(t^^j^). We can express 
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(3-6-47) 


^ ^ (t.J ^ ^ — TTT^ 

— o x+1 —O X 9s(t.) -X 

— X 

vhere the second term is obtained as follows: We can obtain the matrix 

M (t.., , t,) 
s x+1’ x' 




(3-6-48) 


as a subset of the matrix _£ in Eq, (3-6-32). Then combining 

these with S.,, . from Eq. (3-6-46) we can construct the transition 

x+1 ,1 ^ 


relationship 


6m(t^^j^) 


a< 1=1+1 ) 


M M 


6m (t^) 


1+1,1 


(3-6-49) 


The left hand side of Eq. (3-6-49) and the inverse of the transition 
matrix M in (3-6-30) can be used to obtain the mapped effect of the _s^ 
at time mapped deterministically to time t^. Indicating the func- 

tional dependence on we write 

6m I (s.) ” ^ 6m I (s,, ) 

-oL-^-i 1 1+1.0 -L.^1 1 


“lii.o [“ 


1+1.1 


6m(t^) 


= MTii „ M s. 

1 + 1,0 s... . — 1 
1+1 , 1 


(3-6-50) 


which yields the required transition S , 


S ^ M~^(t. , t )• M (t. . , , t.) 
m 1 + 1 ’ o' 1 + 1 ’ i' 


(3-6-51) 
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The state transition matrix ^ of Eq* (3*-6-35) for x is then 
given by 


4,(t.^i. t.) = 


can in general be 


u 


S = 
m 


w 


m 


(3-6-52) 


i+l,i 


(3-6-53) 


where 


^e = 


9-So<"±+l> 

as(t^) 


e = u, V, w 


(3-6-54) 


Considering two cases in particular: 

(i) Non- gravitational accelerations: In this case both 

and S are zero: 
w 

(ii) Camera pointing errors: The entire matrix is zero* 

We now present the equations for processing the ^measurement 

t 

residuals, where we will need the partial derivative matrix of 
Eq. (3-6-33) and the state transition matrix ^ of Eq* (3-6-52).^ 
Sequential Filtering of the Data 

The Kalman filter equations can be derived in various ways and 
are readily availatrle from several sources (e.g. Refs. 34 thru 37). 

Here we just present these equations; we then present the equations to 
consider the sensitivity of the estimate to unestimated parameters. 

In this subsection we drop the underscoring of vectors — all lower case 
quantities are vectors unless other\<rise clear from the context. 
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Given the meastiremeni: model 


■z. = X* n. 

1 IX X 


(3-6-55) 


with the' initial state a zero mean ^aussian random vector, 
E[x(0)] - 0 

and positive semidefinite assiimed covariance matrix 


(3-6-56) 


E[x(0) x(t))^] - X(0) 


(3-6-57) 


and n^ a zero mean guassian white sequence 

E[n.3 = 0 


(3-6-58) 


with covariance 


fifn nTl = r! 6.. 

L 3- jJ 3 ij 


(3-6-59) 


where is a positive semidefinite matrix, the optimal estimate of x 
is given by 


X. - X. + K.(z. - H, X.) . 

X X X^ X XX 


(3-6-60) 


The gain is computed using 


JC. = X. H? Fh. X. hT + rM 
a ii[_xix aj 


-1 


(3-6-61) 


where x^, X^ are the parameter estimate and covariance mapped from the 
previous measurement -epoch through the equations 


^i+1 " ^ 


“ * rj 


(3-6-35) 


and 


^i+1 " ^i + 1, i) + 




(3-6-62) 
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Here 


m 


^(i + 1, i) = 


(3-6-52) 


J i+lji 


and 


r[e,ej] - q. 


13 


(3-6-63) 


The covariance of x^, assuming is optimal, is then given by 


X. = X, - K.H.X. 

1 X XXX 


Consider Covariance Option 


(3-6-64) 


When the system model or apriori statistics assumed are inaccurate 
the filter becomes sub-optimal yet the formal covariance, without 
knowledge of the inaccuracies, x^70uld tend to be overly optimistic. The 
generalized consider option (Ref. 4) involves considering the effects 
of ignoring certain parameters from the estimated state vector and of 
assuming incorrect apriori statistics. 

Let the observable equation be of the form 

z. H.x. + G.y, + n. ^ (3-6-65) 

where the parameters are the consider parameters which are deemed to 
affect the data but are not formally estimated. Let the actual apriori 
covariance be given by 

-1 r 

P 


P(0) = E 


— 






X 





0 






r T 

Tl 





y 

= 



0 

oj 



^O 

L. 



— 

- 





X 






( 3 - 6 - 66 ) 


and the actual covariance of the noise n be 


Efq .qTl = S.-6 . . 

L 1 3 J 3 ^3 


(3-6-67) 
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With the assumptions (3-6-55) , (3-6-57) and (3-6-59) the filter 
gain in C3-6-61) is subop timal and the computed covariance will be 
overoptimistic. Rewriting (3-6-65) and (3-6-60) as 

'x" 


- [H G] 


yj 


+ n. 
1 


( 3 - 6 - 68 ) 


and 


A 

X. 

X 


■ 1 

X. 

1 

' + 

'k.“ 
. 1 

A 




0 



^i 



(z. - H.x.) 
1 11 


(,3-6-69) 


we see that for the complete model the data partial derivative matrix 
and the (subop timal) filter gain are, respectively 

k [H G] (3-6-70) 

and 


A 


K. 

1 


0 


(3-6-71) 


For a filter of this structure, for any filter gain, the update consider 
covariance equation is 

= [^ “ ^<3rl P^r^ - (3-6-72) 

where ^ is the identity matrix; that is 


P. - 

X 


(I - K.H.) ! -K.G. 

^ XX XX 


i 

I 

t ' 
j 
1 


P. 

X 


(I - K.H.) 
11 

-K.G. 

0 

\ I _ 


K, 

1 


r[kJ I o] 


(3-6-73) 
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• The consider -covariance mapping equation is 



( 3 - 6 - 74 ) 


where 

§ ^ §(i + 1, i). 
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3.7 Optical Data Measurement Equation Linearization 


The description of the measurement model x^as given in Section 3.6. 

In this section the linearization of the optical data measurement equa- 
tion is performed; that is, the partial derivative matrix H of Equation 
(3-6-33) is obtained for the optical data. 

We begin with a general form of the optical data equation* and 
proceed to linearize it, developing the partials in terms of the 
direct and indirect terms ; for this formulation the terms are con- 
veniently grouped, in that ,the direct terms all stem from instrument 
related parameters while the indirect terms stem from celestial geom-' 
etry related parameters and are dependent on the celestial object imaged. 
The celestial geometry related components of the latter terms are dealt 
with first, for each of star, planet and satellite images. All space- 
craft and planet related dynamic parameters are obtained as the solu- 
tions from the integration of the respective spacecraft and planetary 
variational equations in Section 3.6. Satellite related parameter 
variation effects arc developed for the satellite theory described in 
Section 3.1. 

Finally the variation of the data equation is taken to yield the 
partials for the vectors upon which the data is directly functionally 
dependent, thus completing the linearization. 

Optical Data Partial Derivatives 

From the discussion in Sections' 3.2 and 3.5 for the optical data, 
the observed pixel and line image locations on the vidicon are func- 
tionally given by 

z == r, b) (3-7-1) 
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where b is a vector of instrument parameters; the direct functional 
dependence on _r arises through the transformation in Eq. 3-2-5. 

If t, T are the spacecraft observation time and light-time respectively, 
and j^(t - r ) is the position vector to the object being viewed at time • 
(t - t), then _t Is a function of ^(t - t), j;(t) and r(t) as discussed 
in section 2.2-. 


Representing all solve-for parameters by jn, we can write 
z = z[t, r(t, m)', m] { 


where 


^ m), r(t, m), x(t - ®)] 


Differentiating Sq, (3-7-2) w.r.t, m, we have 

9m 9t 9m ^ 9r Bm I 9m ( 

— — — I Jt — I -Jt-, r 


(3-7-2) 


(2-7-3) 


(3-7-4) 


where the subscripts indicate the quantities held constant in the 


partial differentiation, and 
9t Bt 9r(t, m) 


d_t 9jt(t, ra) 


8m 9r(t, m) 8m 8f(c, m) 8m 


8x(t - 


fr9x(t - T, m)"| 1 

'ttJ) ■ [ - h' - a)— „ r 


(3-7-5) 


The term 3T/9m must be obtained from the solution of the light time 
equation, Eq. (2-2-3) where w^'are _^Ct - t) respectively* 

Here hoxirever, we can afford to neglect the small effects due to 
stellar aberration and light- time, for the purpose of computing the 
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partial derivatives, A considerable simplification of the expression 
for _t and 9_t/3m results; we obtain 

t m) "^(t, m) (3-7-6)' 

3^ on) 9^(t, m) 

= ^ 7 (3-7-7) 

0JI 3m 3m 

Using the equations 

j: = £ + £ (3-7-8) 

and 

y = £ + £ (2-2-2) 

in (3-7-6), the planet vector drops out, to give 

t = £(t, m) - £(t, m) (3-7-9) 

and 


9£ S£(t, m) 0£(t, m) 

9m 9m 9m 


(3-7-10) 


Combining this with Eq, (3-7-4) and using Eq* (3-7-8) we obtain 


the complete equation 


Bz ra 

\zl "0£ 

8s- 

~d^~ 


”9_p 9£“ "9£^ 

am a 

£ _9m 

amj 

ar_ 

t 

_0m ^ 9m ^ ^ 9m 


(3-7-11) 


The terms 9£/9m and 9^/9ni are obtained from Eqs. (3-6-9) and 
(3-6-13) for all the dynamic parameters m that perturb £ and £ respec- 
tively* The latter, 9£/9m, is obtained from the solution to the space- 
craft trajectory variational equations. 

When the planet is the center of integration in the heliocentric 
phase these yield the entire term 9r;/9m. We develop the terms 9£/9£, 

9_z/ 9£ and 9£/9m later, after discussion of the 9£/9m terms for each of 
the observable image types. 
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Observable Image Types ,and 9_t/8m 

Star Observations: For star observations jt is the heliocentric 

direction to the star and can be taken to be constant since 

(i) the accuracy to which the star locations are known 

(-"0*1 arcsec) is more than an order of magnitude better 
than typical spacecraft television camera resolutions and 
(ii) the effect of 6_r is negligible compared with the distance 
to the star; therefore 


Qt 

— = 0 , (3-7-12) 

om 

for stars, and the first term in Eq. (3-7-11) drops out. 
Planet Observations: In this case _t is given by 

- s , (3-7-13) 

the vector ^ and the partial derivative matrix 3^/8m are identically 
zero, i.e. 


a^ 

9m 9m 


(3-7-14) 


which is obtained using Eq. (3-6-13). 

Satellite Observations: For observations of the natural satellite 

_t and a_t/am are given by Eqs. (3-7-9) and (3-7-10). We now develop the 
term needed in Eq. (3-7-10). The parameter vector m here consists 
of the satellite orbital elements for the Struve-Wilkins-Born Satellite 


theory described in Section 3.1, 

At a given time, t, measured in days past epoch, the position vec- 
tor ^ of the natural satellite of a planet in the 1950.0 Earth mean 
Equator and Equinox coordinate system can be written in the notation of 
Section 3.1 as 
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5. = N(I +-6N) J(I + 6J) K(I + 6K) ffi(I + 6ffi) W(I + 6w) r 

(3-7-15) 

where ^ is the planet relative satellite radius vector in the orbit 
plane, I is the (3 x 3) identity matrix, and the small angle error 
rotation matrices -are defined as 



OBIQlNAIi PAGE IS 

OS 


141 



This equation can be rewritten, keeping only first order terms, as 
£ = NJKffiWr + N(6N)JKffiWr_ + NJ(6J)K5:Wr + NJK(6K)ffiWr 

+ NJKffi(6K) Wr + NJKtqW(6W)r (3r-7-16) 

which can be rewritten, following the method of Section 3.5, as 

£ = + QN6N + QNJ6J + QNJK6K + QNJKffiSffi + QT6W (3-7-17) 


where 


= NJKffiWr = T£ = 


with 


and 


T 4 NJKKW 


Q i 


0 

~q. 


-q. 


0 

-q: 


0 


12 'll 

The error vectors are defined as: 
0 


6N = -c 


0 


6 J = ^ 


5N. 

A 

J 
r 

61 


6J, 


6ffi = 


6W = 


0 

0 

5W, 


(3-7-18) 


(3-7-19) 


(3-7-20) 



6K = 4 


0 


9 
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Writing the last five terms of Eq. (3-7-17) as a 
taking variations 

s 

6 ^. = 65 ^ + E 6 E^ 

i=l 

The 6E^ are obtained by taking variations on Eqs 
(3-1-18). From Eqs. (3-1-12), 

“a - 

“a ' ”2 + 

'^‘'a = ^ V 

6 P = 6P, + 6P„t 

The term 6,W^ is obtained as follows: combining Eqs. 

(3-1-18) and taking variations 

= L - + (V - M) 

6W^ = 6L - 6K^ - 6N^ + 6(v - M) 

From Eq. (3-1-17), 

6(v - M) = (2e cos M + 2.5 e^ cos 2M) 6M 
+ (2 sin M + 2.5 e sin 2M) 6e 

and from Eq. (3-1-18) 

6M = 6L - 6P 

The variation on Eq. (3-1-13) yields 6L, 

6L = 6L^ + 6L„t + 5L„t^ + sin(K. - K„) 

Z N M KAO 

+ ■='’=»A - ‘‘0) “a 


summation and 
(3-7-21) 

(3-1-12) through 

(3-7-22) 

(3-1-15) through 

(3-7-23) 
(3-7-24) 

(3-7-25) 
(3-7-26) 

(3-7-27) 
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The term 6^ of Eq. (3-7-21) is obtained using Eq. (3-7-18), 


6q = T 


(3-7-28) 


where from Eq^ (3-1—20) we can obtain 

2 


6r 


= (.. j- “ e ] 

\1 4- e cos -v/ 


6a - 


a£2e + -(e + 1) cos v] 
2 

(1 + e cos v) 


6e 


, r e sin v . 

+ — - 6v 


(3-7-29) 


a(l - e ) 

Collecting and combining the appropriate variations in Eqs. 
(3-7-21) through (3-7-29) we .can write. 



(3-7-30) 


where 


^ *-2 ^2 ‘'z “z -^zhiW \ \ hi hi 

(3-7-31) 

and is the required partial derivative matrix. 

Variation on the Optical Data Observable Equation 

To obtain the 9_z/9(*) terms of Eq. (3-T-H) the variation is taken - 
of the optical state observable equation, developed in Section 3.2. 

Of the three remaining terms of Eq. (3-7-11) the term arises 

from the variation of the transformation due to variations in the 
spacecraft trajectory, since the ABC coordinate system is defined with 
respect to the spacecraft-sun direction. Thus clearly the effect is an 
angle change which is inversely proportional to the spacecraft-sun 
range. For example, a spacecraft position error of 100 km at 
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1 astronomical unit from, the sun would have aa effect of the order of 
(100 km/1 AU) radians »0»7 x 10 ^ radians; this effect, is therefore 
negligible for typical spacecraft position errors, -relative to the 
other error sources. 

The terms (9^/9_t) and thh direct instrument related terms, 
(9,z^/9m), are evaluated in Ref, 38 and will not he developed here. 
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3. 8 Optical Navigation System Structure 

In this section we give a brief description of the structure of 
the optical navigation system. A schematic flow diagram for the optical 
navigation system is sho^m in Figure 3-6. 

A parameter initialization program is used to create a file of 
basic optical parameters to be used throughout the system. The SAO 
star catalog is used to generate an appropriate subset of stars and the 
integration of the planetary and spacecraft dynamical equations yields 
the planetary ephemeris and a spacecraft trajectory; these are input to 
the Optical Geometry Program to perform the picture sequence design 
function, with the Plotting Program used for the necessary illustration 
of picture sequence geometry in the celestial, inertial and television 
camera coordinate systems. Once a picture sequence design is finalized 
spacecraft commands are transmitted for its execution. 

After picture requests are implemented, the data is processed to 
extract pointing and imaging optical measurements. The Optical 
Geometry Program generates a picture sequence file for use in the 
Optical Observables and Partials Program which computes the residuals 
and perforins the linearization aromd the nominal trajectory. For this 
purpose in addition to the planetary ephemeris and spacecraft trajectory^ 
the planetary ephemeris partials and spacecraft trajectory variations 
are required; these are generated by integrating the respective varia- 
tional equations. The data partials and residuals generated are used 
in the Estimation and Mapping Program to perform the Calibration or 
Orbit Determination functions. The data residuals both before and after 
the fit can be displayed by the Plotting Program, which can also display 


146 




4 >* 

•-vj 


. 5 1 

o ^ 
io ^ 











the spacecraft trajectory estimate history mapped to the B-plane. In 
addition to the spacecraft trajectory estimate the filtering process 
produces updates for parameter values and picture pointing estimates. 
If necessary the entire process can be repeated with these new nominal 
conditions for iteration to convergence. When the trajectory has 
satisfactorily converged it is used to compute the maneuver commands 
required to correct the spacecraft course. 
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CHAPTER 4 


OPTICAL NAVIGATION APPLICATION TO MARINER IX DATA 

In this chapter the first two of the three optical navigation 
methods outlined in Section 1.4 are evaluated using actual flight 
data obtained from the Mariner IX Mission to Mars. Tne objective 
of this mission was to place the space probe into orbit around Mars 
so as to carry out scientific investigations related to the planet. 

The data used for the evaluation here consisted of selections 
from the Science TV pictures taken during the approach phase to Mars. 
These included pictures of Mars taken with short exposure times (10-20 
msec) and of Phobos and Deimos, the two natural satellites of Mars, 
with the maximum exposure times possible (6.144 sec). Several stars 
were imaged, along with Phobos and Deimos, in the latter pictures, 
enabling the satellite- star method for optical navigation to be 
evaluated. However the exposure times for the Kars pictures was 
too short for the stars to be imaged. The TV -pointing for these 
pictures was therefore obtained by using the telemetered engineering 
data. There was no data available for the two- camera method since 
the Mariner spacecraft had only one narrow- angle TV camera. 

In the following discussion we first describe the data in detail 
in Section 4.1 followed by a summary of instrument calibration results 
in Section 4.2. Sections 4.3 and 4.4 respectively deal with the 
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analyses of the satellite- star method and the planet limb method 
(without stars). The results show that the satellite- star method 
yields greatly improved performance over “radio- only" navigation. 

Also, as expected j the planet- limb method is not found to be as good 
as* the satellite-star method. However the results indicate that it 
does contribute to an improvement in accuracy relative to “radio- only" 
data; it would thus be a viable choice in situations where natural 
satellites may not be available. 
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4.1 Data Sources and Content 


In this section we describe the data obtained from Mariner XX 
and used for the optical navigation analyses performed. The sources 
of data were the digitally transmitted TV pictures along with supporting 
engineering data* 

Onboard Measurement Systems 

Instruments onboard the Mariner spacecraft which were utilized 
to obtain the data included the narrow and wide angle science TV 
cameras and the attitude control sensors. TV cameras were mounted 
on a scan platform provided with two degrees of freedom. The 1.1-deg 
X 1.4-deg and ll-deg x 14-deg £ieId-of--view vidicons with 9.6 mm x 
12.5 mm selenium targets were electronically scanned in 700 lines with 
832 picture elements (pixels) per line. Continuous video intensity 
level was sampled and digitized to 9 bits (512 levels) prior to 
transmission to Earth. Each picture element , therefore ^ would be 
defined by its pixel number (1 to 832) , its line number (1 to 700) 
and its intensity (0 to 511). 

The attitude control subsystem provided spacecraft stabilization 
and orientation. Its celestial sensors, i.e.. Sun sensors and star 
(Canopus) sensor, produced pitch, yaw and roll position signals. 

The attitude control signals, together with the scan platform gimbal 
position data, were transmitted to the ground station through the 
engineering telemetry channels. 

Data Types and Content 

The first few pictures during ttie approach phase to Mars came 
from two sequences of planet pictures known as Mars Calibration 
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Sequences I and II. These had been taken to calibrate the TV cameras 
using Mars images and to safeguard against possibilities of damage 
later from excessive exposure; they are different from the scan 
platform calibration pictures which were taken about two months before 
encounter. 

During the 3 day approach period prior to the insertion of 
Mariner 9 into orbit about Mars, three Preorbital Science picture 
sequences (POS I, II, and III) were taken (Fig. 4-1). Each POS 
sequence covered a 24- hour period with 31 pictures being recorded 
aboard Mariner 9 and then transmitted to Earth during a 3- hour period. 

Planet Limb Data: From the two Mars calibration sequences which 

took place 5 and 4 days before Mars encounter (abbreviated as E-5d 
and E-4d, respectively), 22 narrow angle TV pictures which contained 
Mars lit limb images were taken. Wide angle TV pictures of Mars 
during the Mars calibration sequences were not processed because of 
poor angular resolution. During POS (preorbit Science) I, II, and 
III sequences which started from E-3.1d, E-1.9d and E-0.7d, 
respectively, 36 narrow angle TV pictures and 4 wide angle TV 
pictures were selected for lit limb data processing, where POS-I 
and II pictures would have been available before Mars Orbit Insertion 
(MOI) . Towards the end of POS-II and throughout POS-III sequences, 
the Mars images taken by the narrow angle TV were so large that only 
a small portion of the lit limb was visible; hence such frames 
were not processed. 

Satellite- Star Data: The positions of Deimos as viewed from 

Mariner 9 against the star background are shown in Figure 4-2, 
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Figure 4-2. Approach Deimos/Star Geometry 


154 




The first Deimos picture frame was bad, leaving only 12' Deimos 
pictures available before MOI. The POS III phase yielded five 
additional Deimos pictures. 
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4.2 Instrument Calibration. Summary 

Two sequences of pictures, viewing dense star fields, had been 
taken within two months before Mars encounter. These two sequences 
of 31 pictures each, which were aimed at star clusters, single stars. 
Mars and Saturn, provided a good set of data for calibration of the 
scan platform subsystem, the attitude control sensors and of the 
geometric distortion in the TV cameras. A summary of the results 
is given here. 

Scan Platform Calibration 

Calibration of the scan platform subsystem was necessary to 
establish the true scan platform orientation (hence TV pointing 
direction) relative to the spacecraft. The calibration of the scan 
platform subsystem and attitude control sensors was performed by 
referencing to the orientation of the narrow angle TV imaged stars 
whose locations were known apriori. The scan platform errors were 
characterized by systematic errors such as gimbal mounting errors, 
gimbal axis misalignments nnd TV instrument mounting offsets placed 
on the scan platform. Attitude control sensor signal values were 
needed to define the spacecraft attitude relative to the celestial 
references. Systematic errors of the telemetered attitude control 
sensor signals were characterized by sensor null offsets and scale 
factor errors. Mechanical misalignments of the wide angle TV were 
calibrated with respect to the narrow angle TV. 

The star-referenced inflight calibration .(Reference 39) determined 
the acan platform subsystem and attitude control errors to a total 
accuracy of 0.005° (1 a) about the TV line of sight (LOS) and 0.020° 
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in cone and cross- cone axes for the narrow angle TV camera and 0.015° 

(1 a) and 0.17° (1 a), respectively, for the wide angle TV camera. 

This then determined the camera pointing error levels for the planet 
limb data type where stars are not available for accurate pointing. 
Television Calibration 

Calibration of the TV camera geometric distortion was necessary 
to establish undistorted image direction relative to the TV coordinate 
system. This was performed by utilizing both reseau marks etched 
on the vidicon target and the star fields from the two calibration 
picture sequences . Since, ground calibration of the TV cameras 
(Reference 39) revealed practically no optical distortion, optical 
distortion parameters were not solved for. The star and reseau 
images were fit to the analytic geometric distortion model described 
in Chapter 3 . 

The geometric distortion,, was calibrated (Reference 4(J) to 
the accuracy of 0.5 (1 a) pixels and 0.7 (1 cr) pixels in the line 
and p^ixel directions, respectively for the narrow and wide angle TV 
cameras. Figure 4-3 illustrates the changes of observed reseau 
locations in narrow angle TV pictures determined by the calibrations. 
These are reasoned to have been caused by the absence of Earth’s 
magnetic field, and possibly a slight displacement of TV beam deflection 
coils. 
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4.3 Analysis of Satellite-Star Data 


In this section we present results of the analysis of the satellite- 
star data from Mariner IX. The analysis consists of an evaluation of 
the data for orbit determination, followed by its sensitivity to 
number of stars in the field of view, TV distortion, planet ephemeris 
error, and the amount of radio tracking data. 

An approach picture is shown in Figure 4-4 which has the images 
of Deimos, ten stars of magnitude 3.9 to 9.2, and the 7x9 reseau 
grid. Stray light from Mars is seen in the lower left portion of 
the picture. The picture was enhanced to bring out the dim images. 

Figure 4-5 is a computer drawn version of the picture which is used 
to distinguish the star pattern from noise or vidicon blemishes. 

Orbit Determination Using Satellite-Star Data 

The optical data was processed in three iterations (Reference 23). 

In the first iteration, star images were processed to remove TV 
pointing errors. Pointing errors were generally large enough to 
introduce nonlinearities into the data processing because of the 
nonlinearity of the TV distortion model. In the second iteration, 
star and satellite images were processed. Residual TV pointing 
errors from linear corrections in the first iteration were removed. 

Also, Deimos ephemerides and the spacecraft trajectory were estimated 
to reduce Deimos data residuals within a linear region. A trajectory 
based on radio data only was used in the second iteration; however, 
a loose apriori uncertainty (thousands of km) was used. The third 
iteration produced combined spacecraft trajectory and Deimos 
ephemeris estimates using Deimos data and radio Doppler data. After 
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Figure 4-4. 
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Figure 4-5. Predicted Approach Deimos Picture 
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a Deimos ephemeris update from POS I data, the second and third 
iterations could be combined. 

The optical residuals (difference between observed and computed 
image locations) for Deimos (D) and stars (o) are shown in Figure 4-6. 

The residuals in Figure 4-6a were obtained using apriori Deimos 
ephemerides, a trajectory based on radio data to encounter minus 19 
hours (E-19 hours), and TV pointing based on reduced spacecraft 
telemetry data. The clusters of star residuals reflect the global 
offset of the images due to TV pointing errors. It is noted that the 
TV pointing errors become large at the end of POS II and throughout 
POS III in the pixel direction but were smaller and more random in 
the line direction. 

The residuals in Figure 4- 6b were generated using the same 
conditions as those in Figure 4-6a with the exception the TV 
pointing errors have been removed by using the star images. Note 
that the star residuals which are only sensitive to TV pointing errors 
are zero mean. With TV pointing errors removed, the Deimos residuals 
reflect Deimos ephemeris errors (20-hour period) and spacecraft 
trajectory errors. The Deimos ephemeris errors are evident in the . 

30- hour periodic cycle seen in POS II and POS III residuals'. The 
spacecraft trajectory error is seen as a slope in the line residuals. 

The residuals in Figure 4-6c were generated after solving for 
the spacecraft trajectory and Deimos ephemerides. A 100 km Deimos 
ephemeris error and 80 km spacecraft trajectory error were determined. 

The residuals are seen to be random with zero mean and a standard 
deviation of less than 0.5 pixels (3 arc sec). 

iORIQINAL, PAOS m 
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Figure 4-6. Deirnos Observation Residuals 
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Figure 4-7 shows the B-platie trajectory estimates which were 
generated* "The radio only solution and its 1 a error ellipse were 
based on data to E-13 hours. The radio plus optical iterated solution 
and 1 a error ellipse were based on radio data to E-19 hours and 
optical data from POS I and POS II. The current best estimate is 
denoted as GBE» 

Sensitivity to Stars 

The orbit determination (OD) accuracy sensitivity to the number 
of stars per picture was investigated. The following three cases 
were studied: A) no stars/picture ; B) one star /picture; and C) an 

average of five stars/picture. All three cases had apriori TV 
pointing information from reduced spacecraft telemetry data. Also, 
the nominal trajectory was based on radio data only from E-30 days 
to E-16 hours. OD accuracies for the three cases are shown in 
Figure 4-8 and the associated trajectory estimates are shown in 
Figures 4-9 through 4ril- 

Figures 4-10 and 4-11 show that the first picture in POS I for 
cases B and C drives the trajectory estimate to within 15 km of the 
GBEo The time behavior of cases B and C are very similar, with the 
trajectory estimates to within a 2 km agreement at the end of POB III 
data. The expected accuracies of cases B and C (Fig, 4-8) are the 
same . 

TV pointing derived from spacecraft telemetry is an order of 
magnitude less accurate than pointing derived from star images. This 
degradation is reflected in both the expected accuracy and the actual 
trajectory estimate of case A (Fig, 4-9) as compared to cases B and G. 
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Figure 4-9> Trajectory Estimate Using No Stars 
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The case A trajectory estitaate is well behaved when compared to its 
expected accuracy until the last picture in POS II ♦ Then the large 
non- random TV pointing errors (modeled as random errors) in the 
remaining pictures drives the trajectory estimate to a 3 a error. 

The full accuracy potential of the optical data can be obtained with 
only one star/picture . These star sensitivity results can be explained 
by examining the TV pointing errors • For a given picture, all sources of 
pointing errors can be modeled as three independent rotations about the 
axes of an orthogonal coordinate system (e.g., TV pointing has three 
degrees of rotational freedom) . One star image (a pixel and line obser- 
vation) yields two of the three degrees of rotational freedom. The third 
degree of freedom is obtained from a second star or from apriori TV 
pointing which has an accuracy of a few tenths of a degree. Apriori TV 
pointing to this accuracy can be obtained from either reduced spacecraft 
telemetry data (0.015 deg- la) or from the desired pointing (0.15 deg- la) • 
Additional stars, however, do not improve the Deimos- inertial reference 
information in a picture. Deimos image location measurement errors con- 
trol this accuracy and are not affected by star observations. 

Even though only one star/picture is needed, it is desirable 
to have many stars/picture from a reliability standpoint. Many 
stars /picture give independent checks on the TV pointing and also 
indicate the accuracy of the TV distortion model. Any discrepancy 
between image location residuals within a given picture would flag 
it for further evaluation. 

Sensitivity to TV Distortion 

A comparison was made of navigation performance as a function 
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of the data type used to estimate TV distortion and also the order 
of polynomial used to model distortion* Data used for distortion 
calibration included: 1) only reseau images from ground pictures; 

2) only reseau images from flight pictures; and 3) only star images 
from flight pictures. Distortion polynomials of 1st (linear) and 
3rd order were compared to results from the nominal 6th order 
polynomial . 

In comparing calibration data, it was found that all three types 
gave equivalent trajectory estimation results. The TV distortion 
did not change from prelaunch to the end of the mission. Only a 
linear shift and rotation of the reseau grid of a few pixels was 
measured, which was easily absorbed in the TV pointing error model. 
Approximately 200 star images from Plieades pictures and optical 
navigation pictures were used to produce equivalent results as from 
reseau data. Stars are a more desirable data type for distortion 
calibration because they also enable the calibration of optical 
distortion, in addition to electromagnetic distortion and are more 
easily detected than reseaux. 

A tradeoff for increased optical data linearity at the expense 
of reduced model accuracy was evaluated. It was found that the 
linear and 3rd order distortion polynomial gave equivalent trajectory 
estimates to the nominal 6th order model. The increased linearity 
was accompanied with a slightly noisier trajectory estimate behavior 
which, however, was well within ghe predicted accuracy. Therefore, 
it is concluded that a 1st order (linear) model would suffice if 
time constraints would not allow iteration of the optical data. 

ORIQINAL PAGE B 
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If time is available , which will generally be the case, the 3rd 
order model gives the full accuracy potential of the optical data 
with iteration of the data. 

Sensitivity to Mars Ephemeris 

One major source of error in the use of Doppler data for navigation 
estimates is planetary ephemeris errors. This results from the fact 
that the data is taken by stations on Earth and must be related to 
the target planet by using assumed station locations and a planetary 
ephemeris. However, from onboard optical data the spacecraft state 
is directly related to the target planet. 

To demonstrate the independence of optical navigation estimates 
from the planetary ephemeris errors, a solution was made with a 
Mars ephemeris error of about 500 km. The results of processing the 
optical data with this ephemeris error are shown in Figure 4-12, 
giving the B-plane solution history. The origin of the plot is at 
the current best estimate (GBE) , The figure shows that the first 
pass through the data moves the estimate from an apriori more than 
500 km away to within 10 km of the GBE* The final iteration moves 
the estimate to within 2 km of the GBE. 

Sensitivity to Radio Tracking Data 

To evaluate the strength of optical data, a trajectory solution 
was made without the aid of any other tracking data. From Section 2.4 
it is expected that POS I and II data alone would an accurate B • R 
and B • T solution but limited pictures and observed parallax would 
degrade the time of flight accuracy. However, POS III data containing 
both Deimos parallax and trajectory bending would yield a complete 
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Figure 4-12. 


Sensitivity to Mars Ephemeris Error 
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trajectory determination from only the optical data. 

A nominal trajectory was generated from Atlas/Centaur injection 
conditions. These injection conditions gave a 25,000 km aim point 
bias at Mars for planetary quarantine. The use of this trajectory 
did not allow the optical data to "knpw,” apriori, that a midcourse 
maneuver had been performed five days after launch. The midcourse 
maneuver changed the actual trajectory aim point by 15,000 km and 
arrival time by about a day (250,000 km). In other words, this 
apriori trajectory gave the optical data only a vague indication 
that the spacecraft was going in the vicinity of Mars. 

Initially, only POS I and II data was iterated. These results 
could have been available in real time since the POS III data was 
played back after insertion* Five complete iterations of the POS I 
and II data were needed to obtain a converged solution because of 
the nonlinearity due to the 250,000 km apriori trajectory error. 

After a converged solution was obtained for the POS I and II data, 
an additional solution was made which included the POS III data. 

This final solution allowed the full potential of the approach optical 
data to be evaluated. 

The B- plane trajectory estimates are shown in Figures 4-13 and 
4-14 at the end of a complete iteration of POS I and II data. It 
can be seen that an accurate estimate of B • T and B • R can be 
obtained using only POS I and II data as expected. The B- plane 
estimate after 5 iterations was within ten km and ten seconds of 
the current best estimate. Adding the POS III data brought the B- 
plane estimates from optical data only to within five km and three 
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Figure 4-13. Trajectory Solution from Injection 
Using Only Optical Data 
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Figure 4-14. Flight Time Solution Using Only Optical Data 
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seconds of the current best estimate. 


The time of flight estimate and expected uncertainty from the 
final solution are shown in Figure 4-15. It is seen that the uncer- 
tainty does not go below a few seconds until 10 hours from encounter. 
This level of accuracy would be available about a day before encounter 
from radiometric data. It is concluded, therefore, that the optical 
data by itself only can yield an accurate trajectory estimate, using 
data within 10 hours from Mars encounter. By combining optical and 
radio data, an estimate of comparable accuracy can be obtained much 
earlier. 
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Figure 4-15. 


Flight Time Estimate Using Only Optical Data 
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4,4 Analysis of Planet*>Limb Data 

In this section we present results of the analysis of the planet 
limb data from Mariner IX. This consists of the processing of the 
raw limb data to locate the center of the planet image, followed 
by an evaluation of the spacecraft orbit determination performance 
using this data, and its sensitivity to the amount of radio tracking 
data. 

Limb Data Processing 

A typical picture is shown in Figure 4-16. As described in 
Section 3.4^ the limb-finding algorithm searched the digital video 
data TV scan line by scan line. That is, the limb search started 
from appoint off the planet .image and continued onto the limb. The 
first of three adjacent pixels, all having video intensity levels 
exceeding a predetermined threshold level, was selected as the 
location of the lit limb on that line. The requirement for three ' 
adjacent pixels eliminated the detection of false limb point caused 
by one- or two-pixel bit-error noise (Fig. 4-17). 

Typical lit limb structure imaged by the narrow angle camera 
(Fig.* 4-17) indicates a relatively sharp increase in brightness near 
the equator and a gradual increase at the pole regions. The 
transition region between the dark background and the apparent limb 
ranged from 15 km to 25 km wide near the equator, and 25 km to 30 km 
wide at the pole regions. 

I / To. model the imaged lit limb. Mars was assumed to be an ellipsoid; 

4 'I ’ 

as discussed in Section 3»4, with this assumption the image of M^s can 


176 





Figure 4-16. Typical Mars Lit Limb TV Picture 
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Figure 4-17. Mars Limb Image Structures 


be reasonably well modeled by an ellipse, with the center of the 
ellipse being coincident with 1:he center of the ellipsoidal planet 
bodyo 

The radius estimates obtained from the lit limb data (Fig. 4-18) 

are approximately 100 hm greater than the mean, equatorial radius 

(3393,4 km) determined by the S-band occultation experiment (Reference 41). 

Similarly, the identified shape parameter of the planet (Fig. 4^19) 

is smaller than the expected value (e = 0.145) by about 3 percent. 

That is, the Mars images were somewhat larger and more circular than 

expected , These results also indicate that the lit litnb of the 

imaged planet does not represent the Martian surface but light 

scattered from the Martian atmosphere: They may also have been 

caused by the limb darkening effect at the polar regions and lightening 

near the south polar cap, marked albedo change near the lit limb, 

atmospheric activities, i,e., the dust storm which prevailed in the 

Martian atmosphere throughout the approach phase, and the spectral 

response characteristics of the TV camera subsystem. 

The limb fit residuals, which are a measure of the observation 

data noise after a successful limb parameter estimation process, were 

evaluated for each scan line containing a limb measurement (Fig. 4-20) 

2 

with the units of (pixel) , The observation data noise had been 

caused by quantization in picking discrete limb points, uncalibrated 

geometric distortion of the TV, residual limb image from the previous 

picture (triggering a false limb data) , and limb model approximation 

error. To interpret the limb residual statistics with nearly zero- 

2 

mean and a variance of several tens of (pixel) , the residual random 
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Figure 4-18. Observed Mars Mean Equatorial Radius 
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process is mapped onto an equivalent random process along the limb 
image in the TV coordinate system^ From Eq« (3-4-37) we obtain 

Var (60) = Var (6p) + Var (6 jS) , (4-4-1) 

p ^ 

where 


2 , 
c - 4 
P 


1 ) X cos ^ ^ y sin 


(4-4-2) 


c2 = 4 


1 c , 

1 2 ) ^ ^ d* y cos ^ 

a 


(4-4-3) 


If it is assumed .that the observation noise affects equally in line 
and pixel directions, Eq. 4-4-1 can be rewritten as 


Var (6p) = Var (6jJ) = , 


(4-4-4) 


or for typical values for the known parameter relations, e.g*, 


^ = 0.145 , 

a 


and 

2^2 2 
X y p- a 

Eq. (4-4-4) is approximated by 

Var^ (6p) - Var^ (6jJ) , (4-4-5) 

Eq, (4-4-5) typically took a value of 0*3 pixels, reasonably well 
coinciding with the limb quantization error statistics (a uniform 
distribution over a pixel) . 

Orbit Determination Using Planet Limb Data 

The TV data through POS II was processed using as apriori a 
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trajectory based on radio tracking data processed through E--13 
hours. The picture data from Kars calibration sequences I and II and 
POS sequences I and II were used for the analysis. Of these pictures 
the Mars calibration II pictures (at E-4 days) were dropped due to 
vary large data residuals caused by large pointing errors. 

Figure 4-21 shows the B- plane solution history of the trajectory 
estimate obtained by processing the optical data* The parameters 
estimated were the spacecraft Cartesian state at epoch, TV pointing 
biases and errors, and image center finding errors. The origin in 
the figure is at the current best estimate (CBE) . The apriori 
value shown for the optical data arc was the aposteriori estimate 
from a short radio data arc (5 days) through E-13 hours. The estimate 
stabilized only towards the last few pictures. Two iterations were 
performed, indicated by the dashed and full lines, respectively, 
with the last two points on the full line representing the estimates 
obtained from adding POS III data. The error ellipse for the final 
estimate at the end of POS II data had a semimajor axis of 70 km. 
However, the estimate was within 25 km of the current best estimate 
(CBE) . 

Figure 4-22 (a) shows the raw data residuals^ in the pixel and 
line directions, for the POS data« Figure 4-22 (b) shows the residuals 
after the solved for pointing biases and trajectoiry error have been 
removed o The magnitude of the pointing biases estimated agreed with 
scan calibration results o The center-finding errors, modeled here 
as being proportional to the image size, appeared to be very small 
and independent of phase angle o The pointing errors in the raw data 
included systematic as well as random components o 
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Figure 4-21. Trajectory .Fs.timate Using Fit Fixnb Data 
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Figure 4-22. POS Lit Limb Data Residuals 
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The raw data residuals for the Mars calibration sequence X 
pictures are shown in Figure 4-23, The expanded time scale facilitates 
the identification of attitude control limit cycle motion of the 
spacecraft with a period of about 1.2 hours, Thi^ residual limit 
cycle motion may be due to scale factor errors in the attitude 
control model. These were, however, not separately included in the 
estimation but were absorbed in the pointing errors. 

The orbit determination (OD) accuracies (Fig, 4-* 24) indicate 
a slow reduction of the SMAA and SMI A with time up to about E-50 
hours. Rapid reduction in uncertainty occurs thereafter. At the end 
of POS II the SMAA is about 70 km while at the end of POS III it 
is about 30 km. The time of flight uncertainty shows little 
improvement, as would be anticipated from our discussion of the 
approach trajectory. 

Sensitivity to Radio Tracking 

As indicated before, optical data suffer a serious limitation 
in that an accurate determination of the time of flight and is 
not obtained. However, for the Mariner IX trajectory geometry these 
are precisely the parameters that can be accurately determined from 
radio data» A combination, then, of optical and radio data, to give 
better estimates of encounter conditions much earlier than either 
data separately could, would be ideal. To study the reliance of 
optical limb data solutions on apriori values obtained from radio 
data arcs, two different analyses were carried out. In the first, 
a few hours of radio tracking was assumed beyond the first midcourse 


187 








maneuver at launch (L) -h 5 days, and no more radio data thereafter. 

The spacecraft state estimate and covariance at that time were input 
as apriori estimates to the optical data arc. The mapped B-plane 
uncertainties were very large at this point owing to large maneuver 
execution uncertainties being mapped several months forward. The 
analysis was conducted using only the TV pictures through POS II 
for the first pass. 

The B-plane solution is shown in Figure 4-25. The apriori 
estimate was more than 2500 km from the GBE in the B-plane. The 
first pass moved the estimate at the end of POS II to within 450 
km of the GBE. The next pass did not appreciably move the POS II 
solution, but the POS II solution, but the POS III data moved the 
estimate to within 30 km. 

The second analysis was performed with no radio data at all. 

The injection conditions, for purposes of planetary quarantine (PQ) , 
had a 25,000 km aim-point bias and 22 hours lag in time of flight 
at Mars. The analysis made here for this optical-only solution was 
not cognizant of the mid- course maneuver performed at L -F 5 days 
which removed the PQ bias. The B-plane trajectory solution (Fig. 4-26) 
stabilized at about 7000 km from the CBE, with negligible improvement 
in time of flight. 

The large final error involved here indicates that for optimum 
results using optical limb data, some radio data is necessary. Even 
small radio data arcs when used to give apriori estimates and covariances, 
improve the effectiveness of optical limb data enormously (Fig. 4-25). 
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CHAPTER 5 


OPTICAL NAVIGATION APPLICATION TO VIKING 

The Viking Mission consists of two unmanned spacecraft. Viking 1 
and Viking 2, launched in August 1975, which will be inserted into Mars 
orbit in June and August 1976 respectively and will softland instru- 
mented packages on the surface of Mars* The primary purpose of the 
Mission is to further our understanding of the origin and evolution of 
the solar system and of life in the solar system* In this chapter we 
evaluate the application of optical navigation to the Viking mission 
during the approach phase to the planet. 

In section 5.1 we discuss optical navigation related mission design 
aspects, and considerations for selection of data types to be used. 
Section 5.2 evaluates the navigation accuracies obtainable using a 
combination of radio and optical data; this is followed by "optical 
only" accuracies in section 5.3 where we also examine the effect of 
various different data arcs in combination. In section 5.4 a method 
is developed to evaluate consistency between any two solutions and is 
shown through a simulation application to be very effective in the 
detection of data anomalies. We conclude the chapter in section 5.5 by 
the description of a parametric probability analysis technique devel- 
oped to evaluate navigation performance as a function -of system 
reliabilities. 
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5.1 Optical Navigation Considerations in Mission Design 

In this section we briefly describe the mission to place in 
perspective the navigation critical times. This is followed by a dis- 
cussion of the optical navigation data types selected for each phase of 
the mission and the rationale for doing so. We close the section with 
the possible impact of some of the optical navigation considerations 
on trajectory selection and mission planning. 

Typical Interplanetary Mission Profile 

Most space missions go through a similar pattern- of events. 

Shortly after launch the spacecraft is placed into a parking orbit 
around the earth from which it is then injected into the desired 
trajectory away from the earth for the remainder of the mission. At 
the current state of the art the major portion of the space flight 
occurs primarily in "free fall" — i.e. under the gravitational influ- 
ence of bodies in the solar system. This free fall is interrupted at 
selected times to alter the trajectory through spacecraft propulsion 
maneuvers; these are the navigation critical times. 

For interplanetary trajectories, midcourse maneuvers are required 
shortly after injection into the heliocentric interplanetary orbit and 
shortly before arrival at the target planet. Typically a minimum of 
one and a maximum of two midcourse maneuvers are designed at each of 
these two .phases -• the near earth phase and the planet approach phase. 

The near-earth maneuvers are needed to correct for an aiming bias 
imposed at injection due to the planetary quarantine (PQ) constraint 
requirement. In addition to removing any remaining aimpoint PQ bias, 
the planet-approach maneuvers are needed to correct any trajectory 
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errors prior to reaching the planet* .The last approach maneuver 
designed determines the accuracy of the spacecraft delivery at planet 
encounter and provides the first opportunity to control this delivery; 
it is referred to as the "control" point. 

To insert the spacecraft into a planet centered elliptical orbit 
a maneuver is essential near the periapsis of the approach hyperbolic 
trajectory of the spacecraft relative to the planet. At some time 
prior to this, referred to as the "knowledge" point, the best orbit 
estimate^ known is used to calculate the commandable quantities of the 
maneuver; the precision of this determines the accuracy of the orbit 
insertion. The knowledge point and the control point are the times 
when navigation accuracies are most critical for the planet approach 
phase of a space mission. 

For Viking an approach maneuver could be executed from 30 days 
before Mars Orbit Insertion (MOI) at the earliest to 10 days before 
MOI at the latest. 

Television Cameras and Data Type Selection 

The Viking spacecraft orb iter (Fig. 1-4) contains two identical 
narrow angle television cameras the boresight of which is offset in the 
cross-cone direction by 1..38 degrees. The field of view of each 
camera is 1.51 degrees, by 1.69 degrees with a* scan of 1056 lines by 
1204 pixels /line yielding a resolution of 25 micro radians. The focal 
length is 475 mm and the exposure time varies up to a maximum of 
2.66 seconds. The two cameras, can be shuttered at a minimum time 
spacing of 4.48 seconds. 
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The investigation in Reference 22 demonstrated that the satellite- 
star method ^method (ii) of section 1.4 ^yields the best results^ 
followed by the planet-star two camera approach, at both the control 
and the knowledge phases. Howevere as we discussed in section 3.3 a 
constraint that must be met when designing pictures is to ensure a 
separation angle of twice the accuracy of pointing obtainable i.e. a 
separation angle of 1.0 degrees for Viking. Figure 5-1 shows the plot 
of the separation angle between the Mars limb and Deimos for the 
Viking 1 trajectory. The 15 hour sinusoidal variation is observed 
because the spacecraft approach is inclined relative to the orbit of 
Deimos and the period of Deimos is approximately 30 hours* 

As can be seen from Fig. 5-1 the earliest that we can safely take 
Deimos star pictures based on the pointing constraint is about 
MOI - 104 hours. The semi-major axis of Deimos is 23500 km while that 
of Phobos is 9400 km; the availability of Phobos-star pictures therefore 
occurs considerably later. This explains why Deimos is the better 
satellite to concentrate on for approach optical navigation purposes at 
the knowledge phase. Moreover since the control point is at MOI - 
10 days we must rely on the planet-star two camera method for the con- 
trol optical navigation. These planet-star picture pairs are referred 
to as diads. These can be extended to star-planet-star triads; obtain- 
ing the planet picture frame pointing direction by interpolating between 
the pointing obtained from the star picture frames we can thereby 
eliminate the effect of attitude control system limit cycle motion 
during the time elapsed between camera shutters. 
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Trajectory Selection Considerations 

For radio navigation an important geometrical consideration during 
the Mission planning process of selecting a set preferred trajectories 
is the low declination problem. As we discussed in section 1.3 this 
could cause severe degradation in radio navigation accuracies; Viking 2 
radio only navigation accuracies are worse* than those for Viking 1 for 
primarily this reason. 

For optical navigation there is a counterpart consideration during 
the mission planning process. This is in the availability of stars in 
the field of view that are bright enough to be imaged by the vidicon. 
Thus, for example, if the sensitivity of the camera was such that it 
could view no objects dimmer than, say 7.5 visual magnitude, then there 
must be stars of this magnitude or less available in the picture 
sequence planned. This can be a constraint when designing the control 
point picture sequence. 

During the planet approach the spacecraft follows approximately a 
hyperbolic trajectory relative to the planet. The solar gravitational 
influence on' the spacecraft and on the planet (neglecting its mass) is 
about the same since both are approximately in the same orientation and 
location relative to the sun. Therefore the inertially referenced 
spacecraft hyperbolic asymptote changes very little during planetary 
approach. The T.V. look direction to Mars differs from this because 
of (i) bending of the spacecraft hyperbolic trajectory and' (ii) the 
parallax angle between the spacecraft velocity and look directions. 
Since these effects do not become appreciable until very close 
approach to the planet>. the inertial look direction does not change ^ 
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appreciably (see Figures 5-2a and 5- 2b) during the control phase. 
Therefore the star background variety is quite limited; if the star 
availability for particular trajectories turns out to be sparse, this 
could pose a severe problem for control optical navigation for those 
trajectories. 

This problem does not occur for knowledge phase picture sequence 
design since the satellite motion around the planet provides adequate 
look direction variability and therefore star background coverage (see 
Figure 5-3)* However, another problem could arise here for approach 
geometries where the spacecraft trajectory is in a plane parallel to 
the satellite orbit plane. In that event viewing of the satellite 
would be prohibited for a major portion of the satellite orbit, due 
to the presence of the bright planet, possibly leading to (i) limited 
regions of star background variability and (ii) limited mean anomaly 
coverage which may hinder solving for the satellite orbit. 
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(a) Mars Right Ascension 



-40 -35 -30 -25 -20 -15 -10 -5 0 

TIME FROM ENCOUNTER, days 


(b) Mars Declination 

Figure 5-2. Spacecraft Based Mars Inertial Geometry at Control 
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5.2 Radio Plus Optical Navigation Accuracy Analysis 

This section describes a study made to evaluate the orbit 
determination (OD) accuracies obtainable during the interplanetary 
approach phase o£ the mission using both radio and optical data. The 
evaluation is conducted using a strategy yielding accuracies deemed 
representative o£ the CD performance expected. Several combinations of 
data arcs are analyzed to enable the assessment of the sensitivity of 
overall navigation performance to the timing of the approach midcourse 
maneuver. 

The following subsection gives an explanation of the various 
approach maneuver scenarios to be investigated. This is followed by 
subsections describing the assumed composition of the data and the 
filter configuration employed. Finally the OD accuracies are given, 
mapped to the B-plane at Mars encounter (E) . 

Maneuver Scenarios 

The study consisted of evaluating the OD accuracy before each of a 
series candidate approach~midcourse maneuver epochs and before the Mars 
orbit insertion (MOI) maneuver. 

In conjunction with the scheduled MOI maneuver, four different 
approach midcourse maneuver strategies were conceptualized, namely: 

(i) only one approach maneuver at the control point — 10 days 
before encounter. 

(li) only one approach maneuver at 20 days before encounter 
(iii) only one approach maneuver at 30 days before encounter 
(iv) two approach maneuvers at E-30 days and at E-10 days, 
respectively. 
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These maneuver scenarios are shown schematically below. Of these the 
first is the nominal plan. 



MOI-30 DAYS 
1 

MOI-20 DAYS 
1 

MOI-10 DAYS 

MOI 

(i) 

1 

' 1 ' 

1 

A 

1 

A 

(ii) 


A 


A 

(iii) 

A 



A 

(iv) 

A 


A 

A 


Maneuver Scenarios 


Data Constitution 

For the radio data two-way doppler data was assumed to be available 
from 40 days before encounter (E-40 days) , and simulated as such in a 
multistation tracking strategy. The three tracking stations assumed 
for the simulation are those at Goldstone, CA. (DSS 14), Madrid, Spain 
(DSS 61), and Canberra, Australia (DSS 42). Simultaneous data was not 
assumed during the overlapping view periods from a pair of stations. 

The station switching pattern was based' on elevation angle and specifi- 
cation of a preferential order of station coverage desired. The order 
specified x^as DSS 42, DSS 14, and DSS 61, x^ith all data eliminated below 
an elevation of 15® . 

For the simulation, one range-rate point was sampled every hour, 
with the data weight adjusted to simulate continuous data. In addition 
to the range-rate data, one two-way range point was sampled for each 
pass of each tracking station. The range point was always sampled 
when the spacecraft was at the zenith for that pass. 

For the optical data a 3 day arc of Mars- Stars data was assumed 
available prior ta the mid course maneuver. This xrould consist of a 
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pair or triad of frames using the two-camera approach wherein Mars 
would he viewed in frames from one camera and a star field in those 
from the other. The pictures are near-simultaneous (a few seconds 
apart) with the scan platform held fixed between shutter times. 

An arc of data extending from M-4 days to M-1 day (where M is the 
time of the approach midcourse maneuver) was assumed available at a rate 
of one set every six hours, regardless of when the maneuver occurs. 

Thus for an E-30 day maneuver the arc extends from E— 34 days to 
E-31 days and for an E-10 day maneuver from E-14 days to E-11 days. 

For the knowledge phase, the data consists of pictures of Deimos 
against a star background within the same frame. This is possible here 
since the visual magnitude of Deimos is sufficiently high (i.e, the 
brightness is sufficiently low) to enable both Deimos and some stars to 
be viewed with the same exposure setting without saturating the vidicon. 
A data arc extending from E-3 days to E-18 hours was assumed with one 
picture every three hours. 

Analysis Sfrategy and Filter Configuration 

Covariance analyses for the different maneuver scenarios were per- 
formed using combined radio and optical data. The achievable accuracies 
were established for each of the control and knowledge phases. The 
strategy employed for the combination of the radio and optical data was 
to use the a— posteriori covariance from a radio analysis as the apriori 
covariance for the optical data arc. One of the incidental facets of 
this procedure is that the effect of adding Increasing quantities of 
optical data to a fixed amount of radio data can be displayed and 
examined . 
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The only parameters estimated in the baseline study were the 
spacecraft state. Other parameters were considered as error sources 
but not explicitly estimated by the filter. For the radio analysis 
these included the equivalent station location errors, attitude control 
accelerations and solar radiation pressure reflectivity coefficients, 
Brouwer and Clemence Set III planetary ephemeris parameters (see 
Ref. 42) for the Earth--Moon barycenter and Mars, and the mass of the 
moon and of Mars. 

For the control point optical data arc the parameters considered 
as error sources were the measurement and camera pointing biases and 
the image proportional center-finding errors. For the knowledge point 
processing the optical data related errors considered were the l^ilkin^s 
parameters for the Deimos ephemeris errors and again the camera pointing 
biases. For these pictures the camera pointing biases assumed reflect 
the accumulation of error due to- uncalibrated electromagnetic and 
optical distortion, local distortion owing to the Deimos image read-out 
process and asymmetric image blooming -effects in the vidicon. These 
biases in the control data arc havje, in addition, the errors due to 
the offsets between the two cameras and the motion of the spacecraft 
during the time between the shutter times of the two cameras. The 
satellite ephemeris parameters would normally be estimated in-flight 
along with the spacecraft state as one of the solutions generated. In 
the initial phases of the trajectory estimation, in fact, a greater 
degree of reliance may be placed on this solution than on a solution 
based on estimating the spacecraft state only. However it is reasonable 
to expect that the satellite ephemeris parameters thus estimated could 
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be in error up to a level compatible with the apriori uncertainties 
for those parameters assumed in the analysis here. Thus the results 
presented here would be only slightly conservative. 

The error levels assumed for the various parameters considered 
are consistent with the baseline error model given in Table 5-1. The 
apriori uncertainty at the solution epoch was assumed to be 100,000 km 
(1 cr) for position and 1 km/sec (1 cr) for velocity in each direction. 
Accuracy Results 

The control point, for purposes of this accuracy analysis was 
defined as the time at which the approach midcourse occurs, and radio 
data was assumed up to that time. The knowledge point on the other 
hand was defined as MOI-12 hours; only radio data up to E-2 days was 
assumed to be factored into the knowledge point accuracy assessments. 

With these assumptions the maneuver strategies outlined above 
required a covariance analysis sequence with data arcs as indicated 
below. In each case the last sequence listed for each maneuver was 
for the knowledge case where the optical data was Deimos-Stars, while 
the others were for control accuracy and the optical data was 
Mars-Stars . 




Radio 

Data Arc 

Optical Data Arc 

(i) 

a) 

E-40 days 

to 

E-10 days 

E-14 days to E-11 days 


b) 

E-10 days 

to 

E-2 days 

E-3 days to E-18 hours 

(ii) 

a) 

E-40 days 

to 

E.-20 days 

E-24 days to E-21 days 



E-20 days 

to 

E-2 days 

E-3 days to E-18 hours 

(iii) 

a) 

E-40 days 

to 

E-30 days 

E-34 days to E-31 days 


b) 

E-30 days 

to 

E-Z days 

E-3 days to E-18 hours 
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TABLE 5-1 


Baseline Error Model 


Parameter" 

One-Sigma Error 

Equivalent Station Location Errors 


Spin axis distance (all stns.) 

1*5 m 

Height off equator (all stns.) 

15.0 m 

Longitude (all stns.,) 

3.0 m 

(Correlation between all longitudes = 0*9) 


Attitude Control 


Acceleration (each component) 

1.2 X 10“^^ Wsec^ 

Solar Radiation 


Reflectivity coefficient (each component) 

0.05 

Planetary Ephemeris 


Earth-Moon Barycenter 

20 km 

Mars 

30 km 

Satellite Ephemeris (Deimos) 



0“.l 


0“.l 


1°.0 


0“.0003 

I 

0°.02 


0°.l 


0°.0001 

• 

5'*.0 


O^.OO^, 
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TABLE 5-1 


Baseline Error Model (Contd) 


Parameter 

One^Si'gnia Error 

Satellite Ephemeris (Deimos) (Contd) 
a 

0.1 km 

e 

0.001 

Planetary Masses 
GM Mars 

3 2 

0.3 km /sec 

GM Earth 

0.4 kni^/sec^ 

GM Moon 

0.02 km^/sec^ 

Camera Pointing Errors 


Biases (each component) 

0.003 deg 

Random (each component) 

0.02 deg 

Center-Finding Errors 


Image proportional c.f. error 


Cone direction 

2% of image diameter 

Cross-cone direction 

1% of image diameter 

Data Noise 


Radio : 


Range rate (2-way) 

0.015 Hz (= 1 mm/sec) 

Range (2-way) 

100 nsec (= 15 m) 

Optical: 


Star image location (each direction) 

1.0 pixels 

Satellite image location (each direction) 

1.5 pixels 

Mars image location (each direction) 

2.0 pixels 
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Radio Data Arc 


Optical Data Arc 


(iv) a) E“40 days to E-30 days E-34 days to E-31 days 

b) E-30 days to E-10 days E-14 days to E--11 days 

c) E-10 days to E-2 days E-3 days to E-18 hours 

Keeping the initial epoch fixed at E-40 days for the radio data 

automatically builds in a degree of conservativeness desired for the 
earlier midcourse maneuvers. Note that (iv)a and (iv)c are identical 
to (iii)a and (i)b respectively. Thus only seven analysis sequences 
were required to cover all the maneuver scenarios. An additional 
sequence has been added to examine the sensitivity of (i)b to adding 
radio data up to E-18 hours instead of only E-2 days. 

Figures 5-4 and 5-5 summarize the effect on OD of varying manewer 
times. Figure 5-4 shows the control accuracy obtainable as a function 
of days before encounter for the maneuver epoch. The radio only curves 
are also shown for comparison. The curves clearly indicate that the 
maneuver should be performed as late as feasible, from the point of 
view of OD accuracy, particularly so as to benefit from the power of 
the optical data. 

Figure 5-5 shows the effect on knowledge accuracy of moving the 
maneuver time. Also shown are the values for "optical only" OD (to be 
described in the following Section — Section 5.3)» which are unaffected 
by maneuver timing. The x shows the effect of adding more radio data 
beyond E-2 days. It is seen that fo'r Viking 1 there is an appreciable 
effect but not for Viking 2, 

Figure 5-6 shows the variation in Radio plus optical accuracy for 
the knowledge point as a function of days from encounter. It displays 
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b-plane covariance semi-major axis, km 



E - 30 E - 20 E - 10 

EPOCH OF M/C MANEUVER, DAYS BEFORE ENCOUNTER 

Figure 5-4. Effect of Maneuver Timing on Control 
•Radio + Optical Accuracy 
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Figure 5-6. Knowledge Accuracy as Function of Radio and 
Optical Data Taken Concurrently 
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the accuracy (SMAA) as a function of adding radio and optical data 
concurrently • In contrast, all- of the follox^jing plots will show the 
ef-fect of adding increasing quantities of optical data after all the 
radio data has been processed first. 

f*igures 5-7 and 5-8 give the control and knowledge results 
respectively, for the Viking 2 trajectory^ of the strategies and solu- 
tion arcs outlined above. For each .case the plot shows the B-plane 
SMAA and SMIA. The control accuracy results are given first for the 
cases (iii)a, (ii)a, (i)a and (iv)b respectively; next are given the 
results for the knowledge case (i)b, and of the special case to study 
the sensitivity of (i)b to adding radio data up to E-18 hours instead 
of up to only E-2 days. These are followed by results for the knowledge 
cases (ii)b, (iii)b respectively. 

In the Viking 1 E-30 day Control case the SMAA shows degradation 
when optical data is added. This is caused by the bias error assumed 
in the optical data which has an r66 perturbation in the B-plane 
(where r is the range and 60 is the pointing bias) , it is therefore 
smaller for the E-20 day and E-iO day cases where the range is smaller. 
This explains why the effect is not observed for those cases. In the 
case of Viking 2 the computed covariance^ due to data noise, without 
adding the effect due to the sensitivity to error parameters, is much 
higher then for Viking 1 but the bias perturbation is at the same level 
as for Viking 1 (since it is governed by the range). The latter is 
therefore not as noticeable in Viking 2. 

Figure 5-8 shows that for the knowledge case the accuracies 
undergo marked fluctuations with a period of Deimos’ orbit around Mars. 
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The variations are due to errors in the Deimos ephemeris* The fact that 
the consider statistics show this behavior should not be alarming. 

These fluctuations would disappear if the satellite ephemeris were 
formally estimated. However the latter results would be too optimistic. 
The quotation of the consider statistics is an admission that there 
may be unmodelled effects in the data which would limit the ability to 
effectively solve for the satellite ephemeris to accuracy levels below 
those we are using for the analysis. 

The formal results, for control and knowledge, of estimating all 
parameters (instead of "considering” them) along with the spacecraft 
state are shown in Figures 5-9 and 5-10 respectively for the Viking 2 
trajectory. The results in those plots should be regarded as the most 
optimistic expectation of the radio plus optical OD performance. 
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(c) E - 40 d to E - 10 d 


(d) E - 30 d to E - 10 d 


Figure 5-9 


. Radio + Optical Control Accuracies; 
Solving for all Parameters 
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5*3 Accuracy Analysis Using Optical Data Only 

In this section we describe the results of a study to evaluate the 
accuracy analysis obtainable using optical data only^ with no radio 
data at all* In addition to providing significant information regard- 
ing the strengths of optical data alone, having "optical only" solutions 
could serve as a check for radio data; we shall, discuss this consistency 
analysis aspect in more detail in the next section. We shall examine 
the accuracies obtained by various different combinations of optical 
data arcs. 

Nominal Data Arcs 

The maneuver scenarios examined and the optical data constitution 
assumed here are the same as outlined in section 5.2, as is also the 
filter configuration described there for the optical data arc. The 
apriori uncertainty (1 o') for the spacecraft state was assumed at 
10000 km for position and 1 m/sec for velocity in each direction; this 
is not as large as for the radio plus optical analysis, which was close 
to infinite, but is sufficiently loose to enable the navigation evalua- 
tion of optical data by itself; the figures are certainly reasonable if 
preceded by even a small amount of radio tracking. 

The results, in terms of the B-Plane SMAA, of this analysis are 
shown as part of Figures 5-11. and 5-12 for the control and knowledge 
phases respectively. As indicated in Figure 5-11 the curves labelled 
0(1), 0(11) and O(III) give the results for each M-4 day to M-1 
day optical data arc processed by itself. The corresponding, curve 
in Figure 5-12 for the knowledge data arc is the one entitled "nominal 
sequence." The knowledge phase optical only accuracy is only marginally 
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SMAA/ km 



34 32 30 28 26 24 22 20 18 16 14 12 10 

TIME FROM ENCOUNTER, days 


Figure 5-11- Control OD Accuracy Using Optical Data Only 
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worse than the radio plus optical accuracy in section 5,2. The results 
clearly indicate that as expected from our discussion in section 2.4, 
the planet limb data is relatively weak as compared with the satellite 
star data. This is despite the improved pointing accuracies from the 
processing of the star images in the adjoining camera. Of course it 
must be borne in mind that for the planet limb data there is also the 
additional effect here due to the much larger range at which those 
pictures are taken. 

Extended Data Arcs for Control 

Figures 5-11 and 5-12 also show the results of an additional study 
made to evaluate the effect of extending the data arcs. For the con- 
trol case the three data arcs, entitled (I), (II) and (III) in the fig- 
ure, were processed in various combinations as follows: 

(a) I + II + III 

(b) I + III 

(c) II + III 

Also shown in the figure is the effect of Including all three arcs for 
the radio plus optical processing, entitled R + 0(1 + II + III), 

It is immediately apparent from the curves in Figure 5-11 that; 

(i) the accuracies of processing the single data arcs improve 
dramatically for 0(1), 0(11), O(III) respectively — i.e® a 
considerable improvement in navigation performance would be 
obtained by delaying the maneuver epoch as much as possible, 
(ii) the accuracies at any maneuver epoch would improve con- 
siderably with the addition of some more pictures earlier 
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thus extending the total span of the dataware ” e.g. 
comparing O(III) and 0 (II + III) or 0(H) and 0(1 + II). 

(iii) the effect above increases monotonically with the length 

of the total span but the improvement obtained is marginal 
after a point (comparing for instance 0(11 + III) and 
0(1 -b III)). 

(iv) given a total data span determined by the first and last 
pictures, the presence of additional data within the span 
does not contribute at all (in terms of accuracy at the 
final point) , as evidenced by a comparison of curves 
0(1 + III) and 0(1 +.11 + III). 

(v) the radio plus optical accuracy improves considerably from 
the 167 km SMAA of the nominal R + O(III) processing (Fig- 
ure 5-7c) to 100 km using the entire optical span here. 

(vi) the radio plus optical accuracy is only slightly better than 
the optical only accuracy at E-10 days, given a long enough 
span of optical data; this is clear from curves 0(1 + III) 
and R + 0(1 + II + HI). 

The indication from these observations is that although as 
expected from the discussion in Section 2.4 we do need velocity informa- 
tion, this is available by taking a long span of data — i.e. it is not 
always essential to obtain velocity information from radio data. More- 
over the strategy for the data acquisition should be to take a long 
span of dafa for the velocity information and to take the last few 
pictures as late as possible before the. control maneuver epoch for accu- 
rate position information. 
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Extended Data Arc for Knowledge 


Examining the effect of extending the data arc for the knowledge 
case we observe from Fig. 5-12 that 

(i) the final accuracy obtained is the same for the extended 
E-104 hours to E-18 hour arc as it is for the shorter 
E-70 hour to E-18 hour arc. 

(ii) this limiting accuracy is obtained much earlier in the case 
of the extended arc; this consideration could often be 
quite important for estimate reliability considerations in 
the decision for the final orbit insertion maneuver. 

Satellite star pictures would not normally be available much 
before the extended arc in Fig, 5-12 begins. However the planet-star 
two camera data could be taken before that during the knowledge phase. 
A study made to examine the effect of this on the knowledge point 
accuracy yielded results similar in character to the change from the 
shorter to the extended arc of only satellite^star pictures of Fig- 
ure 5-12 — i.e. limiting accuracies were the same but were achieved 
earlier in time. 
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5*4 Radio and Optical Consistency Analysis 


A primary concern for any application of estimation theory is that 
of validating the data — that is, ensuring that the data coming from 
the real world is consistent with the models assumed* We show in this 
section how optical data, in addition to its use with radio data in 
improving navigation accuracy, can contribute in the alleviation of 
this concern by corroborating the validity of the radio data* This 
would thus yield enhanced performance for the total navigation system* 
The method described also serves to establish consistency between the 
radio and the optical data. 

As discussed in Chapter Two the optical and radio data are comple- 
mentary in their information content and their major error sources are 
different. The "best-estimate” of the spacecraft trajectory would 
therefore normally be a radio plus optical solution. However before 
the selection of the best estimate is made the optical only and radio 
only solutions should be checked for consistency. This is necessary so 
as to obtain confidence that there are no significant unexplained 
anomalies remaining and to establish overall consistency between the 
radio and optical data sets. If by this time consistency has not been 
established within the radio only and optical only solution sets 
respectively, the optical data solutions can also aid in the detection 
of inconsistencies within the radio data solutions, as discussed below. 

In the following we first give a general procedure for conducting 
a consistency analysis between any two different estimates of the 
same physical quantity, evaluated in terms of a defined figure of merit* 
Then we develop an analysis for obtaining the covariance of the estimate 
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difference, required to compute the figure of merit, when the two 
estimates are obtained using any two different data types, or their 
combination. This is followed by a discussion of the use of the set of 
figures of merit to give indications of possible problem areas when the 
two observation sets consist of radio and optical data, and the dis- 
crepancy resolution philosophy, if an inconsistency were suspected. 
Finally we give a particular application of this, using simulated 
radio and optical data, which demonstrates the effectiveness of the 
process. 

Consistency Analysis Process 

The following procedure should be conducted for each pair of solu- 
tions to be examined for consistency: 

1) Form the vector difference between the two estimates of the 
system state vector x. 

Ax = _x^ - X 2 (5-4-1) 


2) Calculate the covariance of this difference, A. In general 
this will be a function of the physical model assumed for 
the system and the particular filter or filters used to 
obtain the estimates and ^ 2 * 

3) Calculate a figure of de-merit, f, defined by 


f 4 Axj 


1/2 


(5-4-2) 

The value of this figure of de-merit indicates the extent of 
the discrepancy between the pair of estimates x^ and 
is, qualitatively, an indication of the difference between 
the estimates, normalized with respect to the uncertainty, 
and with all the dimensions collapsed into one scalar measure. 
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Discrepancy Statistics for Different Observation Types 


We now develop the covariance A to be used in Eq* (5-4-2) for the 
case when the solutions x^, x^ arise from the use of two different data 
types. As in section 3*6, let the observation vector z depend on the 
estimated spacecraft state vector x and the unestimated ^consider" 
parameter vector y through the equation 

z'= Ax + By + n (5-4-3) 

where 


E(n) = 0 

E[n(t^) n(t^)] = R6.J 

and where we. have made the notational replacement of A and B for G and 
H respectively, to conform with conventional batch least squares filter 
notation. 

Given any two data types, denoted here by (r) and (o) , the vector 
y can in general be 


fy. 


y4 


( 5 - 4 - 4 ) 


L^b, 

where y^, y^ and y^ are respectively comprised of parameters to which 
(i) only data type (r) , (ii) only data type (o) and (iii) both data 
types, are sensitive. The corresponding matrix B can then be written 


» - I»rl®oPbl 

Let the observations be similarly grouped with subscripts 1, 
respectively; the corresponding equations will then be 

Zj. = AjX + B^y + n^ 


( 5 - 4 - 5 ) 
2 , 3 

( 5 - 4 - 6 ) 
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Z3 = A3X + B^y + 1I3 


(5-4-8) 


E[n. (t.)n. (t, )] = R.6., , i = 1, 2, 3 

1 3 1 k X jk 


rir rfal' 


(5-4-9) 


(5-4-10) 


®2 - ‘\2i'o2i“b2] - t'>|'>02Pb2l 


(5-4-11) 


the zero matrices arising because of the definition of z, , B and of 

1 o 

z^, Since subscript 3 denotes both data equations together, we will 


(5-4-12) 


(5-4-13) 


and, since n^^, are independent white noise, 


(5-4-14) 


The least squares estimate, x, with x estimated and y considered 
from equation (5-4-3) and assuming no apriori on x Is 


A .T „-l . .T „-l 
x=AR A AR z 


T -1 

= XA-^ R 2 
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where X is the ’’computed" covariance of x, 

X ^ E I [x - xjy = 0] [x - xjy = 0]^} 


= [aT r“1 a] 


-1 


0 


Substituting z from Eq^ (5“-4~3) into Eq* (5-4-15) we obtain 


X - X = S y + XA^ R ^ n 
xy 


(i 


where is the sensitivity of x to errors in y, given by 
T -1 

S ^ XA R B 
xy 


(i 


Using Eqs. (5-4-16) to (5-4-18) we can obtain the expressions 
P ^ E[(x - x) (x - x)^] 


S P S + X 
xy y xy 


C 


^ E[(^ - x) y^] 


= S P 
xy y 


C 


or 


S = P ■ P ^ 
xy xy y 


(; 


where 


p = 


p p 

X xy 


xy 


Since y is not estimated is constant, given by 
P = P (t ) 

y y o" 


(. 


1 - 4 - 16 ) 


-4-17) 


1-4-18) 


1-4-19) 


- 4 - 20 ) 


1-4-21) 


1 - 4 - 22 ) 
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For each data set and matrices A^, the equations above 

will hold yielding the corresponding estimate with the estimated 

covariance X. , consider covariance 
x’ 


and sensxtivxty S which will consist of 
x.y 
1-^ 


s is is- 

x.y I x.y I x.y, 
i-'rl 1-^0 1 i-^b_ 


We note, using Eq. (5-4-10) and (5-4-11), that 


(5-4-23) 


(5-4-24) 


The quantity that we need to evaluate relative consistency between 

any two estimates x., x. is the covariance of their difference. A, 

3 - 3 

given by 


= E[(x^ - Xj) (x^ - 


= e|^|(x^ - x) - (xj - x)| |(x^ - x) - (x^ - x)|^j 


= P -P -P +P 
X. x.x. x.x. X. 
3- 3. J J 1 J 


(5-4-25) 


where R is the cross variance between the estimates x. and x., 
x.x. X . 3 

3-3 


■ CvJ 


. (5-4-26) 
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which, using Eq. (5^4-17), gives 


= E r/s y + X.aT R.^ n.\ (s y + X.aT R.^ nA 1 

[_\x^y ^ ^ ^ 7 \ j j j jy J 

= S P + X.aT rT^ E(n.nT) rT^ A.X. 

x.y y x.y ixx \ij/ J JJ 


Since n. are white noise, 
X 


■("l"2) ■ “ 


(5-4-27) 


(5-4-28) 


i "2)] 

= [Vil I 

- h [?ii i ^12] 


1 = 1 , 2 


(5-4-29) 


Therefore the second term of Eq. (5-4-27) is zero for (i 3, j 3); 


for j = 3 it is given by 


[2nd Termlj = 3] - X^aJ R^^ j 6 ^ 2 ] 


R“^ 0 1 [a^ 


0 R. 


-1 A 


= x.A^ (r:^ 

X X \ 1 


Aj 6.^ + 


' 4 *12) 


= X^ , i = 1, 2 


Therefore Eq. (5-4-27) with (5-4-20) yields 


(5-4-30) 


P = S PS 
X 1 X 2 x^y y X27 


-1 T 

= P P P 
^lY y x^y 


(5-4-31) 
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and 


1*3 


where is given by 


= s 

T 

P S 


y ^3y 

= p 

p-1 pT 

x.y 

y X 3 : 


i = 1, 2 


X3 - (xf + " 


(5-4-32) 


(5-4-33) 


Substituting these expressions into Eq. (5-4-25) yields A..; thus 

A = P - P “ P + P (5-4-33) 

12 x^ x^y y x.y x y y x y x 


and 

A . . = P - P P"^ P^ - 2X_ - P P~^ P^ + P 
'•i3 x^ x^y y x^y 3 x^y y x^y 

(5-4-34) 

In particular, using Eqs. (5-4-23) and (5-4-24) in 

Eq* (5-4-31) and assuming a block diagonal form for 


? = S OS 

^1^2 L ^l^r ^l^bj 


0 

0 

_J 


0 




r 



OP 0 


s 

^o 


Vo 

OOP 


s 

L 


L vy 


= S PS 

Vb ^b ^2^b 


(5-4-35) 


which reduces to zero if y^ is a null vector, i.e. if the error sources 
for the two data types are completely disjoint. In this case 
reduces to 

(5-4-36) 


Ai2 ' ^ ^3 


(Xb - 0) 
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Radio-Optical Discrepancy Detection, and Resolution 


The consistency analysis process described above should be 
conducted to examine consistency of the optical only solution with each 
of the following radio solutions; 

(i) doppler only 

(ii) range only 

(iii) doppler plus range. 

The difference between the optical estimate ^ and each of the radio 
estimates j = (i) » Cii)j (iii)j should be formed as in Eq. 5-4-1. 
Since the dominant error sources affecting the radio and optical data 
are quite different, the covariance of this difference, is given 
from Eq. 5-4-36 by 



+ A 


r . 

3 


where A is the covariance of the reference optical solution x and A 

o — o r . 

^ J 

are the covariances corresponding to the different radio solutions 
described above. The figures of de-merit are then computed according 
to Eq. (5-4“2). The magnitude of these will indicate the extent of the 
discrepancy between the optical solution and the different radio solu- 
tions; this could be interpreted to be an indication of consistency for 
small values of f,, or inconsistency for large values of f.. 

3 

Eight (=2 ) different situations are possible from the results of 
these three cases, if the solutions are regarded as being either con- 
sistent or not consistent with the reference optical solution. These 
are tabulated in Table 5-2, where a check (i\j ) or a cross (X) implies 
consistency or inconsistency, respectively. 
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TABLE 5-2 


Consistency Analysis 


Situation No. 

opt. vs Dopp. 

Opt. vs Range 

Opt. vs (Dopp + Eng) 

1 

X 

X 

X 

2 

X 

X 


3 

X 

^/ 

X 

4 

X 

n/ 

n/ 

5 


X 

X 

6 

^l 

X 

n/ 

7 

^l 

^/ 

X 

8 


'v/ 

n/ 

Note: V implies consistency 



X implies Inconsistency 
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The probable interpretations of the different situations are, 
respectively; 

1) Either (a) anomaly in optical data or (b) anomaly (or 
anomalies) affecting both .doppler and range, 

2) Anomalies in both doppler and range, having effects of com- 
parable magnitude but opposite direction, 

3) Anomaly in doppler with possible range effect, 

-4) Anomaly in doppler, 

5) Anomaly in range with possible doppler effect, 

6) Anomaly in range, 

7) Combination strategy for doppler and range sxispect, or 

8) No anomalies — the desired goal. 

As is evident, the doppler plus range- case (case iii) is useful 
primarily in providing supporting information for knowledge acquired 
from cases (i) and (ii) . However, in addition to answering whether or 
not an inconsistency may exist between radio and optical data, the 
values of the figures of de-merit for all three cases can give indica- 
tions of possible problem areas. These values would therefore be of 
significant assistance in analysis. 

The process for the resolution of a discrepancy would necessarily 
be an adaptive one. The following factors ' should be examined: 

1) Solutions obtained by deleting the suspect data type and 
combining the rest, 

2) Solutions obtained by solving for- the parameters pertaining 
to - the dominant- error sources of the suspect data type,. 
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3 ) l^istory of the suspect solutions to see if it shows 
convergent or divergent behavior, 

4) Comparison of residuals of each data type against trajectories 
based on the other data types so as to detect any structure 
within the residuals. 

Application to Simulated Data 

We now present an example for the application of the analysis 
techniques described above 5 using simulated radio and optical data; 
this example illustrates the efficacy of the consistency analysis 
techniques in tracking down problems. 

The first four items in Table 5-3 list the B-plane estimates and 
covariances of representative solutions obtained using 

1 ) optical data only 

2 ) doppler data only 

3 ) range data only 

4) doppler plus range data. 

The corresponding solutions~are shown graphically in Figure 5-13. In 
all the radio data solutions only the spacecraft state was solved for. 
The optical data solution was obtained solving for both state and 
satellite ephemeris. The effect of satellite ephemeris error is 
generally apparent from the optical data residuals (see for instance 
Fig, 4-6b of Chapter 4) therefore this strategy is reasonable. 

The figures of de-meri-t for cases (2), (3)^ and (4) were computed 
using the optical only solution, ( 1 ), as the reference for comparison. 
These are listed in the last column in Table 5-3. It is immediately 
apparent that there exists 
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TABLE 5-3 


Application Example: Radio-Optical Consistency Analysis 


Solution Set 

Estimate j 

X. 

Covariance, A . 

Figure of 
De-Merit, 

fi 

i 

No. ! Description 

1 1 Optical only 

1 
1 

1 

1 

1 

£»R j B-T 

I 

6432 I 6779 
1 
1 
1 

ri91 62" 

L 62 264_ 


— 

2 t Doppler only 

I 

1 

6583 1 7901 
I 
r 
1 


45690 -70397" 

-70397 163423_ 


5.77 

1 

3 1 Range only 

1 
1 
1 

1 

6221 I 7751 

1 

I 

1 

1 


"27242 10974" 

.10974 4577_ 


52.77 

4 1 Doppler + Range 

[ 

1 

1 

6537 1 7921 

1 

1 

1 


"9548 3843" 

.3843 1584. 


65.48 

5 . Range only 

1 (after xange 
1 fix) 

1 

1 

6600 1 6829 
1 

1 

1 


"26235 10544" 

.10544 4388. 


1.36 

6 1 Doppler + Range 

1 (after range 
• fix) 

i 

6877 1 6975 

1 

1 


"9652 3873" 

i 

.3873 I69lj 


4.61 
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(±) an unacceptable inconsistency of the optical data solution 


with the doppler data solution, and 

(ii) a gross inconsistency of the optical only solution with both 
the range only and the doppler plus range solutions. 

The indication from this is that there is an inconsistency between 
doppler data and, range data, probably due to some gross problem with 
the latter; in addition there may be another anomaly affecting the 
radio data (both doppler and range) or the optical data. 

The mathematical representation of the range observable p is 
(Ref. 27) 

p = (At)j^^ F, modulo M 

where 

(At) „ = round-trip time of the signal in seconds of station time 
F == conversion factor from seconds of station time to the 
units of the range observable 

M = modulo number* The largest integer multiple of M which 
is less than (At)^^ F is removed from this quantity, 
leaving the observable p, which is less than M* This 
operation on a number n is referred to as "modding" 
n by M, 

The problem with the data in the example was that the range was 
"out of mod" due to a large error in the apriori trajectory used; this 
caused solutions containing range to have the statistically large dis- 
crepancy with solutions using doppler and optical data* Reintegrating 
a trajectory from an optical plus doppler solution and refitting the 
range data brought the range back within the correct "mod"* 
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The range only and doppler plus range solutions, using the 
reintegrated trajectory, are cases (5) and (6) respectively, in 
Table 5-^3 and in Figure 5-13 ♦ The figures of de-merit, f, for these 
and for the doppler only solution (case 2) show that there still exists 
an inconsistency at an unacceptable level between the optical and the 
radio data solutions* Since the values of f are quite different for 
doppler and range the indication is that there may be an error in the 
radio data which affects doppler data more than the range data ^ 
station location errors. This indication was then confirmed when the 
station location errors were subsequently solved for. 
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5*5 Navigation System Fuel Costs Analysis 

In this section we develop a reliability analysis of the combined 
radio plus optical navigation system from the standpoint of fuel costs. 
This consists of a parametric probability analysis to study navigation 
performance sensitivity of the overall radio + optical system to optical 
subsystem reliabilities. The result of this is a method which can be 
an aid in identifying possible promising areas for enhancing overall 
radio and optical orbit determination reliability. 

The criterion for the measurement of navigation performance here 
is chosen to be the amount of propulsion fuel savings achieved. The 
fuel is measured in units of the absolute value of incremental velocity, 
AV, that can be imparted to the spacecraft. This is normally a major 
concern for a mission since the spacecraft has limited fuel capacity 
and therefore limited capability to make velocity changes. For an ideal 
mission which goes perfectly according to plan with no errors, there 
would be a certain nominal amount of AV needed for required propulsive 
maneuvers. The difference between the nominal and the actual AV 
required in the mission in the presence of statistical dispersions is 
termed There is a total budget provided for these statistical 

variations. The total expenditure of fuel in excess of the nominal must 
stay within the budget, and is a measure of navigation performance 
for a given phase of the mission; is therefore a measure of the 

predicted navigation performance. 

The relative likelihood of failure of any one television camera, 
or of both cameras, can be assessed as a function of spacecraft compo- 
nent reliabilities. Failure of one camera would eliminate the use of 
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the planet-star two-camera data type required at the control point; 
failure of both cameras would eliminate optical navigation entirely* 

The lack of a suitable star background could also eliminate optical 
navigation at :the control point, as was discussed in section 5.1* 

There are therefore four possible optical navigation (ON) system 
operating modes as follows: optical navigation 

(i) operative at both the control and the knowledge point, 

(ii) inoperative — i*e* radio-only at both control and 
knowledge, 

(iii) operative only at knowledge i.e* only radio available 
at control; this can occur (a) if only one camera is 
operational or (b) due to an inadequate star background; 
and 

(iv) operative only at control* 

Let "C" and refer to, respectively, the events of having optical 
data at .control and knowledge; then denote the probability of the four 
cases by, respectively, P(C, K), P(-C, -K), P(-C, K) and P(C, ~K) where 
means "no optical data at;" we then have 

P(C, K) + P(-C, -K) + P(^C, K) + P(C, -K) - 1 (5-5-1) 

However, 

P(C, -K) 0 (5-5-2) 

because this could occur only in the extremely unlikely event that in 
the interval between control and knowledge, both cameras fall* This 
leaves only three optical navigation operating modes and we therefore 
have 

P(C, K) + -K) + P(~C, K) - 0 (5-5-3) 


240 


OF POOB. QUAKEM. 



The total probability density of can be expressed as 

f(AV) = f(AV/C, K) P(C, K) + f(AV/~C, K) PC-C, K) 

+ f(AV/~C, ~K) P(~C, ~K) (5-5-4) 

Based on control and knowledge radio only and radio plus. optical 
covariances, a Monte Carlo maneuver analysis gives AV„„.„ histograms 
for each of the three ON operating modes; these correspond to the three 
conditional probability densities in Eq* (5-5-4) • 

Since the three probabilities in Eq. (5-5-4) are related by 
Eq. (5-5-3), there are only two independent parameters on which the 
distribution of depends. The two parameters chosen for analysis 

here are P(OPT) and P(C/OPT), i.e., the probability of having any 
optical data, and the probability of having optical at control given 
that one has optical at all. The first parameter, which implies the 
existence of any optical data at all, is rather an obvious choice, 

r 

since we are interested in determining the efficacy of the optical 
navigation system in its contribution to the overall reliability of the 
orbit determination, as measured by reduction. The second 

parameter, P (C/OPT), can be used to assess the influence of obtaining 
optical data at the control-point. 

In terms of these two parameters, noting that "OPT** = "K", we 
can write 

P(C, K) = P(C/OPT) • P(OPT) (5-5-5) 

and 

P(~C, ~K) = P(~C/~K) P(~K) 

= [1 - P(C/~K)I [1 - P(K)] 
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which, due to Eq. (5-5-2)j gives 

P(~C, -K) = 1 - P(OPT) (5-5-6) 

Similarly (or using Eq. 5-5-3) 

P(~C, K) - P(OPT) • (1 - P (C/OPT)) (5-5-7) 

Using Eqs. (5-5-5) through (5-5-7) in (5-5-4) and computing 
based on 


"^^STAT 

f(AV) = 0.99 


(5-5-8) 


we can therefore parameterize in terms of P(OPT) and P(C/0PT). 

The results of this parameterization are displayed in Fig- 
ure 5-14 (a) and (b) for a case where the radio OD was considerably 
degraded owing to very large (10 km) station longitude uncertainties 
assumed. Figure 5-14(a) demonstrates the importance of obtaining 
optical data in order to reduce AVg^^^. In the complete absence of 
optical data (P(OPT) = 0) AV^ ^^ ^ is 265 m/sec; (typical allocation 
budgets are -125 m/sec). — On the other hand, with P(OPT) = 1, the 
acquisition of optical data at control reduces AVg^^^ to 26 m/sec. 

Even with optical data available at only the knowledge point 
■97 m/sec. Also, we note that increasing the reliability of successful 
control OD- is ef fective (in the sense of significantly reducing AVg^^^) 
only with a fairly high reliability of the overall optical system — i.e. 
P(OPT) -greater than 90%. 

The sensitivity of TO P(C/0PT) for constant P(OPT) is 

illustrated in Figure 5-14(b); for instance, for a single camera sys- 
tem reliability of P(OPT) = 94% (if at least one camera system is 
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^VsTAT ^ 

P (CIOPT) AS PARAMETER 



P (CIOPT) = ] 









operating, we are assured of optical data at knowledge), AVg 
from 127 m/sec for P(C/OPT) - 0 to 93 n/sec for P(C/OPT) = 1. 


varies 


The marked reduction in for high values of P(C/OPT) indi- 

cates that it would be beneficial to maintain this probability as high 
as possible* In fact, this relates to the star availability discussion 
in section 5.1. We’ have 

P(C/OPT) = P(2 camera ON syst^m/l camera ON system) 

X P (stars) 

= P(2 TV/at least 1 TV) x P (stars) 

. P(2 TV) X P (stars) 

where P (stars) is the probability of viewing stars on the vidicon at 
control. Since this probability may be zero in some extreme cases, it 
would be desirable to plan the flight trajectory so as to ensure a 
suitable star background at control. 
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CHAPTER 6 


CONCLUSIONS 

The central conclusion of the dissertation is that Optical 
Navigation can be a very effective means of navigating an interplanetary 
spacecraft during its approach phase to the planet, particularly with 
the combined use of radio and optical data. More specifically, in the 
following we discuss the principal conclusions of the investigations 
reported hare. 

1) Of the three observation methods examined the satellite-star 
method yields the best accuracies followed by the planet-star 
two-camera method* The former is therefore the preferred 
method for the knowledge phase; for the control phase the 
latter must be used due to camera pointing limitations. 

Optical navigation using the planet limb without stars is the 
weakest method, but does contribute to an improvement in 
accuracy relative to "radio-only” navigation; it would be a 
viable choice in situations where the other two methods are 
not available. 

The primary error sources for the three methods are, respec- 
tively, satellite ephemeris errors for the satellite-star 
method, camera offset biases and planet center-finding errors 
for the two-camera method, and TV pointing errors for the 
planet-limb one- camera method. 
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2) Light- time effects are Important and must be considered, but 
stellar aberration effects can be neglected* This is because 
(i) in the data types using stars, only star relative 
information is sought and differential stellar aberration 
effects are very small and (ii) in the planet limb one- 
camera method, camera pointing errors' are an order of magni- 
tude larger than the stellar aberration effect. 

Relativistic effects can be neglected for the current state 
of the art. 

3) The optical navigation performance is not heavily dependent 
on the number of stars in the field of view when this is 
greater than one. However, calibration of the television and 
of the spacecraft instrumentation, prior to acquisition of 
optical data, is of significant aid in improving data 
accuracy . 

4) The optical data and the radio data are complementary in 
their inherent information content for spacecraft navigation. 
The optical data types are relatively weak in their ability 
to determine velocity and time of flight information; radio 
data determines these same quantities quite precisely. On 
the other hand radio data is weak in the determination of the 
out-of-plane (ecliptic) component of position, where the 
optical data can contribure significantly. 

5) The major error sources of the optical and radio data are 
quite different. For the latter the primary errors stem from 
uncertainties in the tracking station locations, planetary 
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ephemerides and non-gravitational accelerations. Optical data 
is insensitive to these errors. On the other hand optical 
data is affected by uncertainties in the satellite 
ephemeris, TV pointing and camera characteristics, to all of 
which radio data is insensitive. 

6) The accuracies of the optical navigation processes, using the 
three optical observation methods, differ markedly in their 
sensitivity to the amount of radio data used. The final 
accuracy using the satellite's tar method is almost unaffected 
by changes in the amount of radio data; the plane t-s tar two- 
camera method is significantly affected; and the planet limb 
one-camera method is practically useless without some radio 
data. However, the dependence of the two-camera method on 
radio data reduces considerably if the data arc is lengthened, 
and almost disappears for sufficiently long optical data 
spans. In these situations the optimal strategy is to acquire 
data at the beginning and at the end of the data span, with no 
data in between. 

In all cases a combination of optical and radio data yields 
the best results and much earlier than either separately. 

7) From the orbit determination point of view, the maneuver 
timing should be delayed as much as feasible, so as to benefit 
from the strength of the optical data, which increases with 
decreasing range from the planet. 

8) The consistency analysis technique should be applied to the 
* radio and optical estimates. This can help in establishing 
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the validity of each data type or aid in the detec±ion and 
resolution of inconsistencies* 

9) The sensitivity study of navigation fuel costs versus 

optical subsystem reliability shows that enhancement of con- 
trol orbit determination reliability is effective only for a 
high overall optical system reliability (greater than 90%). 
Since this latter reliability should t 3 T)ically be quite high, 
attention to control optical OD, particularly in the area of 
star availability, would be fruitful. 

Looking towards the future of space navigation, the next 'major 
step will be in the area of Autonomous Navigation (see Ref. 43). Since 
this must be a completely self-reliant system with no earth-based 
dependence, optical navigation will play a major role in its develop- 
ment. In addition some means to supplement optical angular data with 
a data type providing range information may be needed in the future. 
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APPENDIX A 


B-PLANE DEFINITION 


The target parameters used for accuracy analysis are the _B plane 
parameters (Ref. 44) defined in Figure A-’l. The uncertainty in these 
parameters is mapped into the 1 cr aim plane dispersion ellipse. 

The following definitions for the B-plane refer to Figure A-1. 

Parameter Definition 

^ The vector from the center of the target body 

directed perpendicular to the incoming asymptote 
of the target centered approach hyperbola. 

S_ Unir vector from the center of the target body 5 

in the direction of the incoming asymptote of the 
target centered approach hyperbola. 

^ Unit vector from the center of target hody, 

defined by the intersection of the plane normal to 
the incoming asymptote with a reference plane 
(usually at the true ecliptic of date) centered 


R 


B • R 
B • T 
SMAA 


at the target body. 

Unit vector perpendicular to ^ in the T plane 
(defined as the B-piane), positive toward the -z 
direction of the reference plane. 

The component of B along R. 

The -component of B along 

The semi-major axis of the one-sigma uncertainty 
ellipse in the R, J plane. 
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Figure A-1. B-Plane Definitions 
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Parameter 


Definition 


SMI A The semi-minor axis of the one-sigma uncertainty 

ellipse. 

6 The orientation angle of SMAA measured positive 

clockwise from T to R. 

The uncertainty of time of closest approach’ to 
the target. 
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